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PREFACE 

The Digital Computer Laboratory of the University of 
Illinois has a three-fold program which is concerned with re- 
search in the field of digital computers, with maintenance for 
University use of a high-speed digital computer, and with teaching 
of design and use of digital computers. The preparation of pro- 
grams for the Illiac is the responsibility of the person who wishes 
his problem solved. But since each prospective user cannot always 
attend the courses on programming offered by the University, it 
has seemed desirable to help make it possible for him to learn 
the elements of Illiac use by himself. 

The result has been this book, first issued in March 195^, 
which makes use of two years of experience in operating the Illiac 
for University research and several years of experience in the 
teaching of programming. Many of the chapters were used in pre- 
liminary form as supplementary material in Mathematics 385/ the 
basic course on programming, and the book is now used as standard 
course material. 

Like most of the work of the Digital Computer Laboratory, 
this is the result of a group effort by many different people. The 
chapters were written by S. Gill, R. E. Meagher, D. E. Muller, 
J. P. Nash, J. E, Robertson, T. Shapin and D.J. Wheeler. Other 
members of the Laboratory staff read the material and offered 
valuable criticism. The typewritten copy was prepared by Helen B. 
Clark and Muriel P. Fetzner; the drawings were made by George 
Ehrlich and Charles Breed. 
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CHAPTER 1 



INTRODUCTION 



In 19^8 the University of Illinois began to look into the 
possibility of acquiring an automatic digital computer. When attempts 
to buy one or to have one built failed, the University decided to 
build one of its own, and the Digital Computer Laboratory was organized 
in February 19^9- Shortly thereafter an agreement was made to build 
a second computer for the Army to be used at the Aberdeen Proving 

Ground . 

It was decided to build the two machines following the de- 
sign of a computer then being designed and built at the Institute for 
Advanced Study. While many modifications were later made, both 
machines' built by the University owe much to the early designs ob- 
tained from Princeton. 

The two machines were built more or less together so that 
advantage could be taken of the savings obtainable by paralleling 
the design and construction work, but work on the Army machine, called 
the ORDVAC, was emphasized, and it was finished first. It passed its 
acceptance tests in February 1952 at Aberdeen and it has been in use 
there since that time. 

Work on the University computer (later named the Illiac) was 
completed in September 1952 and the computer was first made available 
for University use when classes began on Monday, September 22, 1952. 
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The Illiac is an automatic electronic digital computer » 
It is digital because it handles numbers as sets of digits 
which have discrete values, rather than as scale readings or measure- 
ments, which are continuously variable o Apparatus for handling digits 
is more complicated than that for handling continuous quantities, but 
it is capable of giving unlimited accuracy by using suitable numbers 
of digits o 

The Illiac is electronic „ In the last ten years electronic 
circuits for storing, transmitting, adding, subtracting, multiplying 
and dividing numbers in digital form at extremely high speeds have been 
devised. The actual addition of two numbers in the Illiac takes only 
about 75 microseconds. 

Such speed is useless unless the machine can be made to go 
ahead on its own with may thousands of operations, without human inter- 
vention. The Illiac is therefore automatic, in the sense that it can 
be given orders telling it how to proceed, and will then act on these 
orders automatically. 

In common with many other computers of a similar type, the 
Illiac contains the following five essential features" 



(1) An arithmetic unito 

(2) A memory or store, 

(3) Devices for the input and output of 
information (e.go numbers) to and 
from the machine <■ 

(k) Means for the transfer of information 

between the various parts of the machine 
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(5) Means for the automatic control of 
the whole machine. 



The arithmetic unit carries out the individual arithmetical operations 
that make up every computation; it can be thought of as the electronic 
equivalent of a desk calculating machine. It is described in detail 
in Chapter 2. 

The memory is needed because, in any lengthy calculation, 
numbers produced at early stages of the calculation are frequently 
required to be used at later stages; they must therefore be recorded 
or "remembered". The memory is capable of recording 102^ numbers. 
These can be recorded (ioe. transferred to the memory from the arith- 
metic unit) individually, as directed by the computer's control de- 
vice, and recalled again individually in a similar way. The memory 
may be thought of as 102^ little boxes or locations, each accomodating 
one number, and labelled with the numbers through 1023. The label 
of a location is called its address . A number in the memory is 
identified by the address of the location containing it. 

Information enters and leaves the Illiac coded in the form 
of a pattern of holes in punched paper tape; there is a tape reader 
for input, and an automatic punch for output. There is also a tele- 
typewriter which can be used to provide output from the machine 
directly in printed form. Several machines are available for pre- 
paring punched tape, copying it, comparing it, and producing printed 
versions of the Information on it. 



1 

This is the high-speed memory. The magnetic drum furnishes additional 

storage. 
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The problem of controlling the whole computer has been 
solved by stipulating that every Individual operation that occurs within 
the machine must be one of a certain set of specified permissible 
operations^ and that no two such operations can occur simultaneously. 
Thus the design problem was reduced to that of engineering the various 
permissible operations and arranging for them to be executed in any 
desired sequence. It is up to the user of the Illiac to specify the 
sequence of operations or program, which the Illiac must execute to 
carry out his calculation o 

Each permissible operation can be specified in a concise 
coded form called an order o The correspondence between the set of 
permissible operations and the set of orders which specify them is 
called the order code of the Illiac » It is given in detail in Chapter 
S* A coded problem is called a program or routine . 

The machine's control unit has the task of accepting orders 
one by one, and of causing the machine to carry out the operations 
specified according to the order code. If each order were taken by 
the control unit directly from a punched tape, then to make full use 
of the speed of the rest of the machine the tape would have to pass 
through the tape reader at about 200 miles per hour. Instead^ the 
orders are recorded in the memory along with the numbers, so that the 
control unit merely has to take its orders from the memory, which it 
can do electronically at high speed o This is made possible by coding 
each order to look like a number ^ To be more precise, orders are 
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stored in pairs ;, one pair to a memory location. The information 
contained in one memory location is often called a word, meaning 
either a number or an order pairo Of course, the more orders there 
are in the memory, the less room there is for numbers. Both orders 
and numbers are fed into the machine initially on punched tape. 
Normally orders are obeyed by the control unit in the 
sequence in which they are stored in the memory, eog,: 

Left-hand order in location 6, 
Right-hand order in location 6^ 
Left-hand order in location "J , 
Right-hand order in location 7^ 
Left-hand order in location 8, etc. 

Sometimes, however, this sequence is broken and the control unit 
starts over at some new position in the memory; this is called a 
transfer of control o There are special orders which cause this. 
There is also special provision for making a transfer of control 
depend on the value of some number obtained by the machine during 
the calculation. Thus the machine can be made to "choose" one of 
two or more alternative courses of action according to the way 
things happen to work out. 

If control is transferred a few locations back in the memory, 
the machine will repeat the operations specified by the Intervening 
orders. It is possible to cause this repetition to occur any number 
of times, leading to a cyclic behavior of the machine. Practically 
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every calculation which the machine performs contains several such 
cycles, often one inside another. In this way it often happens that 
the same order gets carried out many thousands of times, so that a 
few orders suffice to keep the machine busy for several minutes,. If 
each order in the memory were to be carried out once only, the I Iliac 
would get through them all in a quarter of a second (even if the memory 
contained nothing but orders )o In practice, calculations vary in 
duration from a minute to a few hours o 

The occurrence of cycles is one of the things that compli* 
cates the programming of a calculation^ Another is the fact that^ 
since orders are stored in the memory in the same form as numbers 
they can be operated on and altered during the course of a calculation 
(at the behest of other orders) just as if they were numbers. All this 
makes possible some most interesting calculations; it can also make 
programming difficult. 

Fortunately a coder can of ten, as described in Chapter kj 
make use of bits of programming done by other people. Thus a typical 
program consists of a number of groups of orders, some written by the 
coder, others already available. The latter will be available in 
punched tape form, and can be copied mechanically onto the program 
tape along with the new orders. Tape preparation is described in 
Chapter 9, 

When the whole tape for a particular program has been prepared 
it can be placed in the tape reader of the Illiac, The Illiac reads 
the tape, forms the orders and numbers punched on it and stores them 
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in the memory. When the program is in the memory, the machine begins 
to execute the orders, continuing until it comes to some particular 
order which causes it to stop. If the programming is correct, this 
is the end of the calculation. If there is a mistake in the programming 
various things may happen; remedies are discussed in Chapter 8. 

Somewhere in the program will he some orders which cause 
the machine to punch some output tape. This carries the results of 
the calculation. The program may also contain orders causing the 
machine to read more input tape, carrying data for the calculation. 

The reading of most of the program tape is accomplished by 
the Illiac executing a particular set of orders called the Decimal 
Order Input (See Chapter 5) which is always punched at the beginning 
of every program tape and hence read into the machine before the rest 
of the tape. The Decimal Order Input not only assembles the program 
inside the machine; it also makes certain modifications and conver- 
sions, so that the way in which orders are represented when punched 
is somewhat different from their final form in the memory. The 
object is to make programming easier. It is important to remember 
that the written form of an order and the form which it assumes in 
the memory are not the same thing. The relationship between the two 
is determined by the Decimal Order Input. 

Remaining chapters of this manual are devoted to: The 
arrangement of calculations so that all the numbers encountered 
are the right size (Chapter 6); ways of prograimning certain types of 
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simple tasks (Chapter 7)1 how to estimate the duration of a calculation 
(Chapter lO); an example of a typical program (Chapter il); and how the 
program library is organized (Chapter 13)0 Finally^ concise descriptions 
of the principal contents of this collection are giveno 
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CHAPTER 2 



THE ARITHMETIC OF THE ILLIAC 



The construction of a computer involves a coraprorr'sc 
between engineering economies on the one hand and ease and flexi- 
bility of use on the other. As a result, the details of the 
operations of arithmetic are often dictated by engineering design 
considerations. The following paragraphs describe the peculiarities 
of the representation of numbers and of the operations of arithmetic 
in the Illiac . 

2.1 REPRESENTATION OF NUMBERS- The simplest device for 
storage of numerical information is an electronic element having two 
stable states in which one state is called "zero" and the other state 
"one". One such element is capable of holding one binary digit. 
Registers composed of forty such elements are provided in the Illia.j 
for representation of forty binary digit numbers . The Illiac is 
described as a binary parallel digital computer having a precision 
o ^ forty binary digits — roughly equivalent to twelve decima l 

digits . 

The Illiac is a fixed-point computer^ the location of 

the fixed binary point between the first and second digits is in- 

./ — — ■ — — ' — — — — ~~ 

dicated in Figiire 2.1. 



2-1 



/location of binary point 

3^ 



/ ^ 



■v^ 



sign digit 



59 non-sign digits 



Figiire 2.1 
The Forty Digits of a Word 



The lef -toos-^i git of a number is the ^Ign digit. If the sign 
^ digit is zero, the numb e r is posi tive or zero; if the sign digit is 
^one, the number is negative. For example, + 7/8 is represented in 



the Illiac as 



0. 



11100 . . . 00 = + 1/2 + 1/i^ + 1/8 = +Y/8 



^ Negative Numbers > Negative numbers are represented in 
the Illiac as complements with respect to 2. The process of comple- 
mentation is carried out by forming the digitwise complement (re- 
placing ones by zeros and zeros by ones ) and then adding a unit in 
the least significant (thirty-ninth) non-sign digit. As an example, 
- 7/8 is formed from + 7/8 by complementation as follows: 



0.11100 
1.00011 

+ 0.00000 



00 

11 

01 



1.00100 



00 



+ 7/8 

ones replaced by zeros and zeros by ones 



addition of 2 



-39 



machine representation of - 7/8 



It is sometimes essential to distinguish between the 
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machine representation of a number and its arithmetic value. Suppose 
we have a number whose arithmetic value is x and whose machine re- 
presentation has the sign digit x.^ and non-sign digits x. (i = 1,2,... ,39) 
If X is positive^ i.e., if x^ = 0, then 

39 i 
X = Z 2 X . (2.1) 

i=l ^ 

If X is negative (x = l), then the relationship is 



39 

2""' X..I , 
i=l ^ \ t 



u^ u.^^^^^ ^^ -r 



X = -1 + E 2'^ x..r , , ) (2.2) 



Equations 2.1 and 2.2 may be combined in the one equation 

39 . ^, 

X = -x^ + Z 2 x^. " ' (2.5) 



i=l 

If equations 2.1 through 2.5 are to hold, x must be restricted to 
the range -1 < x < 1. The programmer must keep in mind that the 
Illiac arithmetic unit is designed for numbers in the range -1 < x < 1 
When an operation is performed which yields a number outside this 
range, an unwanted number within the range results. This effect is 
called overflow and, although certain operations are unaffected by 
orerflow, computational errors usually result. 

2.2 FUNDAMENTAL OPEE^TIONS OF THE ILLIAC ARITHMETIC UNIT. 
The following five operations are fundamental in the arithmetic imit 
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of the Illiac: 

1. ccmrplementation, 

2 . addition, 

^' left shifty or multiplica tion by two, 
^' right shift, or division b y two, 
5* clearing to zero. 

Complementation is executed in the Illiac in the manner 
described in the previous section by a circuit called a complement 
gate. The complement gate is actuated by a signal from the control 
and is capable of supplying either the number unchanged, or its 
digitwise complement. Addition is carried out by the adder, which 
is capable of forming the sum of two ko binary digit addends . 

A left shj^t is a displacement of the binary digits one 
digital position to the left and corresponds to a multiplication by 
two. <There will be an overflow if the number to be shifted is 

+ 1/2 or If it has an absolute va lue greater than one -half . The 

■ "^ — — ■ — ' .^ 

right shift is a displacement of the binary digits one digital 
position to the right and corresponds to a division by two. For 
example, 

0.0111 shifted left is 0.1110 5 7/16 x 2 = 7/8 
^0.1110 shifted right is 0.0111; 7/842= 7/16 
l.pOlO shifted right is 1.1001; -7/8^ 2 = -7/16 
1.0100 shifted left is 0.1000; ^3/h x 2 = I/2 because 
of overflow. 
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It should be noted that the sign digit is propagated when the right 
shift is executed. The Illiac arithmetic unit has two shifting 
registers, each capable of executing both left and right shifts. 

Clearing to zero involves setting all digits of a number 
to zero; the corresponding arithmetic value is zero, 

2„3 THE ILLIAC ARITHMETIC UNIT o The structure of the 
Illiac arithmetic unit is shown in Figure 2o2, The arithmetic unit 
is composed of two shifting registers, the accumulator A and the 
quotient register Q, and one non-shifting register, the number 
register R"^. Also required are the complement gate and the adder. 
The A and Q registers are the only registers to which the programmer 
has direct access; the number register R"^ is used to hold temporarily 
the numbers brought from the memory for arithmetic operations. It 
is essential that the programmer be familiar with the roles played 
by the A and Q registers in the operations of arithmetic; many pro- 
gramming errors arise from placing operands in, or removing results 
from, the wrong register. The functions of the registers during the 
operations of arithmetic are described in the sections which follow. 
2»^ ADDITIONS (order Type L) , Before an addition in- 
,, struction begins, the__augen^ lies in the accumulator register A. 



I L ti^ ^ -^ 

'^^ ' During execution of the addition instruction, the addend is trans- 

3 

ferred from a specified memory location to the number register R . 

The digits of R^ are then sensed through the complement gate unchanged, 
so that the addend forms one of the inputs to the adder. The augend 
in A is the second adder input. The adder forms the sum which is 
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transferred to the accumulator k, replacing the augend. The quotient 
.^register Q is undisturbed by the addition Instruction. 
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Figure 2.2 
The Illlac Arithmetic Unit 



Two variants of addition are "hold add" and "clear add." 
! Phe "hold add" instruction leav es the augend in A undisturbed 
until the sum is for med by the adder. The "clear add" instruction 
clears A to zero inlj|:lally^ thus setting the augend to zero. The 
llclear add" instruction is thus a transfer order which moves a 



number from 'a sp ecified memory location to the accumulator A. 
2.5 SUBTRACTIONS (Order Type L). Subtraction in 
the llliac arithmetic unit is performed by a.dding the compleme nt 
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of the sulDtrahend to the minuend . Before the subtract instruction 
begins, the min uend lies in the accumulator k , and is used as one 
of the two adder inputs. The subtrahend is brought from a specified 
memory location to R^; its complement is formed by the complement 
gate and is used as the second adder input o The adder thus forms 
the difference by forming the sum of the minuend and the complement 
of the subtrahend. The adder output is then transferred to the 
accumulator A, replacing the minuend. 

Either the "hold subtract" instruction or the "clear 
subtract" instruction can be used by the programmer. For the former - 
instruction, the result of a previous operation is left in the accumu- 
lator A as minuend; for the "clear subtract," the mi nuend in A is 
set to zero, so that A contains the negative of the number in a 
specified memory location vhen the operation has ended. 



2»6 ABSOLUTE VALUE ADDITION AITO SUBTRACTIOIT ( Order Type iJ . 

It is possible to form the absolute value of the addend (or sub- 
trahend) after it has been transferred from a specified memory loca- 
tion to R^ and before it is added to the augend (or subtracted from 
the minuend) in the accumulator. Ordinarily the setting of t he com- 
plement gate depends upon wh ether the instruction is an addition o r 

3 

a subtraction; for forming the absolute value of the numbe r in R , 

it is necessary to sense the sign digit of R"^ as vei l. For example, 
if the addition of the absolute value of a negative addend is requir ed, 
the sign digit of the addend in R-^ is sensed indicating that a 
complementation is necessary. The add instruction ordinarily does 
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^Ql-i:£gi4ig£..c omplementatlon; the net eff ect of the sensing of the 
_a dd^ Instruction and the R-^ sign digit is that the complement gate is 
/get to form the complement of t he addend. 

2.7 INCREMENT ADD ORDERS (Order Type F). It was noted 
in section 2.1 that the operation of complementation is performed 
by first forming the digitwise complement of the number held in R^ 
and by then adding a unit in the least significant digital position 
of the adder. ^For the ordinary addition and subtraction orders 
^^Jord er type L), the least significant digit insertion occurs only 
^when the complement gate is set to form the digitwise complement. 
In the IlliaCj^ special increment add orders (typ e F) are provided. 
F or these orders ^__thej ^ations hip between the setting of the com- 
pl ^ment gate and the inse rtion oT the least significant digit is 

reversed. Thus,, the "clear increment .add" instruction s ets the 

augend initially in A to zero and adds the addend from R^ to 2"^^ 
^and places the result in A. Conversely, th e "clear increment sub- 
tract" instruction sets the minuend in A to zero, replaces ones 
of_ the subtrahend by zeros and replaces ze ros of the subtrahend 
by ones, and tr ansfers the digitwise complement of the subtrahend 
thus formed from the adder to A . Detailed descriptions of further 
orders of this type are given on page j^-P8 . 

2.8 ADD FROM Q (Order Type S) AND INCREMENT ADD FROM Q. 
(Order Type K) INSTRUCTIONS For the add from Q and increment add 
from Q instructions, the addend or subtrahend is tranf erred to the 
number register R^ from the Q register, rather than from a 
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or-nfied memory location, Oj:;heivlse, the K and S order types 
are the same as the corresponding F and L or der types describe d 
in s ections 2.^ through 2.7. 

2,9 THE SHIFT ORDERS (Order Types 0. l). Since the 
right and left shift operations are fundamental in the Illlac 
arithmetic unit, specific shift instructions are provided for 
the programmer. All digits in both A and Q, except the sign digit 
q of Q> are shifted. A left shift (order type O) of one digital 
position replaces 

^0' ^1' ^2'' ' " ' ^31' ^38' ^39 ^^ ^ " 
and 

^0' ^1' ^2 • • • ^ ^^37' ^38' ^39 ^^ ^ 
by 

a^, ^2' ^-3' • * * ^ ^38' ^39' ^1 ^^ ^ 
and 

qo' ^2' ^3' ' ' - ' ^38' %' ° ^^ ^' 
The right shift (order type l) of one digital position replaces 



a. 



^, a^, a^, . . « , a^^, a^g, a^^ in A 



and 
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^0' ^1' ^' * * • •* ^37' ^38' ^39 ^^ ^ 

by 

^Qf Sq, a^, . . . , a^g, a^^, a^g in A 

and 

^0' ^39/ ^1^ • • • > <l36' ^37' ^38 ^^ ^* 

From the behavior of the left and right shift we see 
that we may consider the shifting to take place in a single register 
AQ of 79 digits consisting of A followed by Q with q^ deleted. 

Thus the left shift replaces 

^0' ^1' ^2' ' ' ' > ^37' ^38' ^39' ^1' ^' • • • / ^3Y^ ^38' ^39 
by 

^1' ^2' ^y ' ' ' > ^38^ ^39' ^1' ^' ^y ' * • f ^-^Qf ^39' ^ 

While the right shift replaces 

^Qf ^i^f ^2' * ' ' > ^37^ ^"38' ^39' ^1' ^' ' ' • » ^3Y> ^38' ^39 

tjy 

Sq, a^, a^, . . . , a^g, a^^, a^gj a^^, q^, • • • > 935, q3^, q38 

The number- n of shifts can be specified in the address 
digits of the shift order by the programmer. A left shift of n 
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digital positions replaces AQ by 2^ AQ. Similarly, a right shift of 
n digital positions replaces AQ by 2~ AQ5 n lies in the range 

1 < n < 65. 

2.10 MULTIPLICATION (Order Type 7). Initially, the 
mult iplier Ij^g-^^j^-^Jae-nunti en-t;, register Q. At the beginning of the 
multiplication instruction, the muO^tlplicand is transferred from a 
location specified by the address of the multiply order into the 
number register R , where it remains throughout the multiplicati on. 
The multiplication then consists of a sequence of additions and shifts. 
More precisely, a multiplier digit in q is sensed. There are two 
cases. (a) If q^Q = 1, the multiplicand is added to the partial 
product in A. A right shift in AQ follows which halves the number 
in A, moves a new multiplier digit into q^Q and transfers 
a product digit from a to q^. (b ) If q^^ = 0, only the right 

shift occurs, which transfers a-,„ to q, and transfers the next most 
' 39 1 

significant digit of the multiplier into q^Q. After thirty-ni ne 
rigfct shifts have occurred, the sign digit of the multiplier is s ens ed 
If the multiplier is positive, the multiplication is c omplete, with 
a double precision product (sign digit and seventy-eight non-sign 
digits) in AQ. For a negative multipli er x, the process of thirty- 
nine shifts and conditional additions constitutes a multiplication 
of the multiplicand y by 1 + x (Equation 2-5) to form a product 
y(l + x) = y + xy. In this case the Illiac automatically subtr acts 
the multiplicand y to produce the correct product xy . 
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A nvunber of variations of the basic multiplication pro- 
cedure described in the previous paragraph are possible. The initial 
contents of A may be either 0, I/2, or some quantity previously 
calculated. If we designate t he initial contents of A by a, the 
final dquble precision product in AQ is xy + 2"^"a, where ^ and y 
are multiplier and multiplicand , respectively. For a/= 0. the 
i nstruction specified is "clear mu ltiply" and the result in AQ is 
t he exact 78 digit signed product xy . For a = I/2, the instruction 
specified is "roundoff multiply" and: the result in A is a rounded 
/39 digit signed product. If a is arbitrary, the instruction 




specified is "hold multiply", and^ the result in AQ is xy + 2"-^^ai 
t tp quantity a is thus added to the least significant part of the 
jproduct in Q . I n. all multiplication instructions th e sign digit 
q^^ of the quotient register is set to zero. ^^^^ 

For each of the three types of multiplication instructions 
described in the previous paragraphs, four additional variations 
can be specified by the programmer. Either N(n), -N(n), |N(n)|, 
or - |N(n)| can be used as the multiplicand, where N(n) is the 
number transferred from memory location n to the number register 
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2.11 DIVISION (Order Type 6). Initially, tl)e doubl e 
.precision dividend lies in A Q. The divisor is transferred a^ the 
beg4.nning of the division instruction from a memory location 
specified by the address of the divide order to the number register 



'^/ 
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B. , For .positive divisor and dividend, the process is analogous 
t p elementary long divis ion. The d ivisor is subtracted f rom a 
partial remainder in A and the sign of the differ ence (in the adder) 

is sensed, -^f the difference is negative, is inserted in q^^ ' 

^ — ^ __ -2- — r-- ^y- 

as quotient digit and AQ i s doubled to form a new p artia l remaind er. 
I f the difference is positive, 1 is inserted in q..,^ as quotie nt 
digit and the difference in the adder is placed in A, and AQ is 



doubled to form a new partial remainder. At each doubling of AQ, 

q is shifted into q^ as veil as into a . Thus after 39 step s <t^^ 

q^ has the sign of the quotient. The dividend is used as the initial 



partial remainder; after 39 quotient digits have been generated, the 
process is complete. A similar procedure is employed if divisor 
or dividend is negative or if both are negative. 

The Illiac division has the following properties: 
(l) A rounded quotient Q is always generated. 

The roundoff is achieved by setting q = 1 

. . i^z 



in all cases. 
(2) The thirty- ninth partial remainder is left 
in A and is called the residue r .. The true 
remainder R co rresponding to the rounded quotient 
is related to r approximately as follows: ^ 

R = r 4. (2q^ - l) y, . (2.^^) 

w here y is the diviso r and q^^ is the sign 
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digit of the quotleni ;. Thus if the quotient 
is positiv e 

_^R = r - y, (2.5) 

and if the quotient is ne gat i ve 

^ = ^ -^-y- / ^ (2.6) 

flihe equations given above for the true remainder 
, are valid if the absolute value of the quo tient is 
^ less than one, and they yield results in error by 
j>ot more th an 2 . 

(3) The sign dijit of the quotient replaces the least 
s j^ificant (78th) digit of the doub le-precision 
. dividend. One effect is that the least significant 
digit of the residue C^.^q) is the same as the sign 
digit of the quotient q . 

{h) If we have a priori knowledge of the true value q 
of the quotient (such as, for example, a division 
of qy by y), the relationship between the machine 
quotient Q and the true quotient q is 

Q = q + 2'39 (1 . q^^) (;l - 2y^) (2.7) 



where q is the least significant digit of q and 
y^ is the sign digit of the divisor y.. Ec[uation 
$.7 is essentially a description of the division 
roundoff. If q^^ = l, then Q = q and the machine 
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quo t i ent is the true quotient . If g,-^ ^ = 
,.the nature of the roundoff depends upo n the 




sign of the divisor. Suppose, for example, that 
q is O.lOli then Q is either 0.101000. . .001 or 
0.100111... Ill depending upon whether the divisor^ 
was positive or negative, respectively. 
2.12 PRECISE CALCULATION OF THE DIVISION REMAIItDER . We 
define th e remainder R in relation to a quotight Q .. . a divisor y , 
and ! ^. dividend ^ hy the equation 

Qy + 2 -"^ R = d. 



or 



TR = 2^9 (d - Q.v). (2.8) 

However, the exact relationship involving the Illiac residue r and 
including the replacement of the least significant dividend digit 
do by the quotient digit q^ is 

Qy + 2"5^ [r + (2qQ - 1 )y] = d + 2-'^^(q^ - d^g). (2.9) 

Solving for R, we have 

R = 2^^ (d - Qy) = r + (2 qQ - l) y + 2"^^ (d^Q - q^). (2.10) 
Equation 2.10 gives the exact expression for the remainder R which 

n 

corresponds to the machine quotient Q. 

If the remainder R corresponding to the true quotient q 
is desired, then R is defined as 
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R = 2 (d - qy). (2.11) 

The value of R is found by substituting equation 2.? in equation 
2.9 and solving for H, which yields 

R - r 4- [2q^ - q^^ - (1 . ^^^) 2y^] y ^ 2-^9(d^Q - q,^ ) . (2.12) 

2-13 THE DIVISION HANGUP. SPECIAL CASES OF DIVISION. 
^^y,'!^^:^^ ^^^ incorporated in the Illiac for stopping the Ill iac if 
the quotien t resulting from a division exceeds one in absolut e value. 
The_sign_dig it_ of the quotient is predictable from the signs of divisor 
an^ dividend. The sign digit of the quotient is ca lculated in the 
^.JJ4 gc "by conrparing the dividend and divisor arithmet ically. Thus, 
by sensing the sign digits of quotient, divisor, and dividend, it 
is possible to detect the fact that the quotient exceeds one in 
absolute value, and stop the Illiac. 

CI The^quations derived in Sections 2'lQ_a nd_2.n are valid 
only if th ^bsolute valu e o f the di^ridend is less than the absolute 
valu e of the divisor . W^iefi_absj)l^ute^^^^^^ dividend and divisor 
_^are eq ual , the Illiac generates a quot ient -1 + 2"^^ if the dividend 
CJL^-JL^g g:^^^^ ^^^ "^^^ divisor is posit iA^e. The quotient is +1 - 2" ^^ 
i f, dividend and divis o r are both negative. I f the dividend is 
^ positive an d equal to the absolute va lue of the divis or, the Illiac 
vii;L stop . If the _divisor is ^ l^^Jhe^^llij^generate.c, a quotient 
w h^ch is the digitvise complemen t of the dividend, except for the 
quotient roundoff. 
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2.1^ MEMORY TO Q (Order Type 5) M^ STORE (Order Type U) 
INSTR UCTIONS . Instructions are provided in the Illiac for trans- 
ferring a number from a specified memory location to the Q register, 
and for transferring a number in A to a specified memory location. 

The fn-rrriRr i nstruction can he us ed to tranF^fer the tttttI tj-pl i ft ■ 

M 

+r, q before a multiplicat ion^ the latter instruction is used_jbo 

tr ansfer a result from A to the memor y . Lf a resu lt lies in Q( e . g. , 
a division quotient) it can b6 transferred to the memory by using tw o 
orders, 35 and 40 which transfer .Q to A a nd then A to the memory . 



2.15 ILLUSTRATIVE EXAMPLES}. ^^HJ-c^ 




A. The Leapfrog I Division Test. In the Leapfrog 
I division test the product xy is formed and is then divided by y. 
The sum of the quotient in Q and r esidue in A is formed and stored. 
( ^e sum of quotient and residue is then calculated independently 
and compared with the sum previously stored . 

The value of the machine quotient Q in the quotient regis - 
ter after division of xy by y is given by equation 2-7, with q = x. 
Thus 

Q = X + 2"^^ (1 - x^^) (1 - 2y^). 
The value of the residue r left in the accumulator is found from 
equation 2.12 by setting R = 0, x = q, and solving for r: 
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r = U^^ + (1 - x^^) 2yQ - 2x^] y + a'^? (x^ - x^^ y^^). 

The sum Q + r, after rearranging terms, is 

Q 4- r = X - 2"^^ x^^ y^^ " ^o (^ ' ^'^"^ ) + (l - x^) y + 
(1 - x^^) (2y^ - 1) (y - 2-^9). 



The independent calculation of Q + r thus consists of forming 

""39 -^Q 

^ " 2 X y and either adding y or (-y + 2"-^^) depending upon 

whether x is positive or negative; and finally, if x is 0, adding 

t "-^Ox 

or subtracting (y - 2 ^^) depending upon the sign of y. 

B. The "Double Precision" Division. It is 
sometimes convenient to consider a single precision divisor y as 
exact and fonn a "double precision" quotient s + 2"^t utilizing 
the double precision dividend d originally in AQ. The procedure 
used is as follows: 

1. Form d/y, yielding s + 2~^^ in Q and r + 2"^^s in A. 

2. Shift right, forming r/2 in A and s/2 in Q. 

5. Form r/2 + y leaving ^^ in Q and a residue u in A. 
4. Assemble s + 2 ^t by setting t (the sign digit of T) 
to zero, inserting s/2 in A and shifting left once. 

-39 
The precision of s + 2 -^^t can then be calculated as follows: 

Step 1 yields s and r such that (by equation 2.9) 
(s + 2-^^)y + 2"^^ [ r'+ 2"^^ s^ + (2Sq . l )y ] = d + 2"^^ (s^ - d^g). 
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sy + 2"^^ (r + 2s^y) = d - 2"^^ d^.g. (2.15) 

Step 3 yields 

Ty + 2'^^ [ u + (2tQ - l)y ] = r/2 + 2"^^^. ^ (2.U) 

The following substitutions are madej 
(a) s ^ t^ since the. sign of r/2 is the same as the sign of 

d. 
(h) From step ^, t = 2(T + s ) or T = l/2t - s . 



Substituting in eq.uation 2.1^ and solving for r, we have 
= (t - 2s^)y + 2"^^ [ u + (2s.^ . l) y ] - 2'^^ s^. 



Substitution of the value for r in equation 2.15 yields 
sy + 2"^^ [ (t - 2sQ)y + 2s^y ] + 2"^^ [ u + (2s - l)y ] 
- 2 Sq = d - 2 d^Q. 



We thus have 

(s + 2'^^t)y + 2''^^ [ 2u + 2(2s^ - l) y + d^g - 2*"^^ s^ ] = d. 

It can be shown that u + (2s - l)y < 1 so that the quantity 

within square brackets is less than 3 in absolute value. We conclude 

-59 -T6 

that the quotient 6+2 t is in error by not more than 2 

The program for forming a "double precision" quotient 
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.from a double precision diyldend and an exact single precisio n 
divisor is given in wor ds ^7 to 66 of li brary routine A k entitled, 
"1.7 Prpcision F l oating DecimaV '. A major difficulty encountered 
is the formation of •r/2 from r (Step 2 above), for r may be as large 
as 2y and may therefore exceed range. It is therefore necessary 
to set the sign digit of r/2 to that of the original dividend d. 

2.16 INTEGER OPERATIONS . * It is sometimes desirable to 
use integers for computations in the Illiac . Suppose we have an 
integer a stored in the memory or arithmetic unit. In terms of the 
formulation of previous sections, we would store 2'^^a, where a 
lies in the range -2''^ < a < 2"^^. If we wish to add or subtract 

two integers a and b, no difficulty is encountered, for 2~^^a + 2"^^b = 

-59 

2 (a + b ) indicating that the correct sum or difference lies in 

A after the ins true ion is performed. Multiplication of two integers 

a, b yields (2"^^a) (2"^^b) = 2'^^ ab . The product ab lies in AQ 

78 7ft 

and is in the range -2 ' < ab < 2 ' . If the programmer scales 

all quantities so that the product remains in the range -2^^ < ab < 2^^, 

then the kO digit signed product can be transferred to A by a left 

shift of 39 digital positions. It should be noted that the sign 

digit of Q is set to zero during the multiplication so that for 

positive products, N(Q) = ab, if < ab < 2^. 

Division of integers presents certain difficulties . An 

example is given here of a method of dividing a positive dividend 

a by a positive divisor b to yield a quotient f and remainder g. 



2-20 



The steps are as follows: 

(l) Place the dividend 2~' a in AQ . _) < a < 2 ' ' , 

(2 ) Shift left one digital position^ leaving 2a in AQ, with 

139 = °> 

(3) Divide by 2'^%, leaving (2f + 1 ) 2"^^ in Q and (2g) 

« 

2"^^ in A. < b < 2-^^, 

-39 
(h) Shift right one digital position, leaving 2 f in Q, 

-59 
and 2 g in A . 



It can be proved that bf + g = a by substitution of the 
appropriate quantities in equation 2.9^ as follows; 

(2f + 1) 2"^^ (2"^\) + 2"^^ [ 2"^^(2g) - 2"^% ] = 2"'^^(2a) 

which yields 

2"^ (2bf + b + 2g - b) = 2"^ (2a) or bf + g = a. 

38 38 

The ranges of f and g are < f < 2-^ and < g < 2 . 

2.17 SUMMARY. In the Illiac arithmetic unit are two 
registers, A and Q, which are directly accessible to the programmer. 
A single arithmetic order of a program utilizes the initial numerical 
operands in A, in Q, and in a specified memory location, and trans- 
forms these quantities to produce desired results which are left in 
the registers . The programmer must know where the operands are 
initially located and where the results are to be found. The functions 



2-21 



of the registers for the operations of arithmetic are indicated in 
Table 2.1. 

The Illiac has a fixed point arithmetic unit; the binary 
point is fixed so that any number x used in computation must lie 
in the range -1 < x < 1. The programmer must insiire that all 
quantities remain within this range during a computation. The 
sign of a numerical quantity is indicated by the leftmost of the 
forty binary digits stored in a register of a memory location. 
The sign digit is for a positive number 1 for a negative number. 
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TYPE 


INSTRUCTION 


INITIAL CONDITIONS 


FINAL CONDITIONS 


NOTES 




Specified 
Memory 
Location 


A 

, Register 


Q 
Register 


A 
Register 


Q 
Register 




L 


Add 


Addend 


Augend 


Arbitrary 


Sum 


Unchanged 


Augend can be 
set to or l/2 




Subtract 


Subtrahend 


Minuend 


Arbitrary 


Difference 


Unchanged 


Minuend can be 
set to or 1/2 


1 
ro 


7 


Clear 
Multiply 


Multipli- 
cand 


Set to at 

start of 

1 natruction 


Multiplier 


Double Precision Product 


q set to 


Rounded 
Multiply 


t> 


Set to 1/2 at 
start of in- 


II 


Rounded 
Product 


Usually 
Unwanted 


ti 




Hold 
Multiply 


ri 


Accumulant 

a 


n 

X 


xy + 2" 


•39a 


II 




6 


Divide 


Divisor 


Double Precision 
Dividend 


Residue 


Rounded 
Quotient 









Left Shift 
(n 1) 


See Note 
2 


Double Precision 

Number x 


2x 


l.q unchanged 

2. Address digits 

specify num- 
ber (n) of shifts 




' 


Right Shift 
(n 1) 


tt 


It 

X 


x/2 



Table 2,1 Use of Illiac Registers During Arithmetic Instructions 



CHAPTER 3 
THE ORDER CODE 

The Illiac is a binary computer in which the storage 
capacity of each register or memory location is ^0 binary digits. 
The orders which the machine carries out are represented by numbers 
in the machine. The relation between the order types and the corre- 
sponding numbers is called a code, and the collection of all such 
numbers is called the order code of the machine. The order code 
is interpreted by the control circuits of the machine and completely 
determines what the machine does. The machine is designed so that 
any storage location in the memory may be used either for orders or 
for numbers, the only distinction being that the control must be in- 
structed properly so that orders and numbers will be treated 
appropriately, 

3.1 THE MAKEUP OF ORDERS . An order for a digital computer 
consists of an instruction to say what to do and one or more addresses 
to say where to get the quantities to be used in carrying out the in- 
struction. In contrast to some existing computers, the Illiac uses 
what is called a one-address code. For a one-address code, the order 
is of the kind that says "add the number in memory location 12 to a 
number already in the arithmetic unit, leaving the sum in the arith- 
metic unit." 

Since a one-address code is used, it is not necessary to use 
kO binary digits to describe an order » The Illiac uses 20 digits and 
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packs two orders (an order pair ) into one location. These are the 
left-hand order and the right-hand order. Since many orders must 
refer to locations in the memory, each order contains an address . 
It is still called an address in those orders which do not refer 
to the memory. 

The electrostatic memory of the Illiac has 1024 locations 
and because 1024 =2 we require 10 binary digits for the address. 
These are the rightmost 10 digits of the 20 digits assigned to an 
order . 

Of the remaining 10 digits of each order, eight are used 
for the instruction or (function ) and the other 2 are unused. The 
digit makeup of an order pair is shown in Figure 3,1. 



8-DIGIT 
INSTRUCTION 



10- DIG IT 
ADDRESS 



8-DIGIT 
INSTRUCTION 



LEFT-HAND ORDER 
DIGITS 0-19 



-w 



10-DIGIT 
ADDRESS 



RIGHT-HAND ORDER 
DIGITS 20-39 



Figure 3,1 
Order Pair Makeup 

The instruction digits are 8 in number because of the 
convenience obtained by using a base l6 (sexadecimal) number system 
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in which k binary digits may "be represented by one sexadecimal digit. 
Thus each instruction may he coded as 2 sexadecimal digits. The 
symbols K, S, N. J, F, L are used for 10, 11, 12, 13, 1^, 1^. 

As an example of an order pair let us consider the following 
^0 binary digits o 

1111010100000001110101000000000000000110 
When divided into instruction and address digits, these digits look 

like this : w * .y 

1111 0101 00 OOOOOl'llOl 0100 0000 00 0000000110 
INSTRUCTION L.H, ADDRESS INSTRUCTION R.H. ADDRESS 
The left-hand instruction is made up of the two ^-digit numbers 
1111 0101 which are the sexadecimal digits L5. The left-hand address 
is interpreted as an integer which may go from to 1023 if we use 
decimal notation or from to 3LL if we use sexadecimal notation. 

In sexadecimal representation the left-hand address is IJ 
which corresponds to the decimal number 29. Thus the left-hand order 
is L501J where the has been supplied so that all 20 binary digits 
(including the unused 2) are accounted for. (We could have set the 
unused digits to I's and used N rather than if we had wished.) 

Similarly the right-hand order is i^■0006, and we have, in 
sexadecimal notation, t he order Dair 

L501J i^0006 
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This order pair says, "Transfer the contents of memory location IJ 
to the accumulator; store the accumulator contents at location 6 
of the memory." 

It is quite inconvenient to have to write addresses in 
sexadecimal form, and it is unnecessary if the machine contains a 
program which will take addresses written in decimal form and con- 
vert them to sexadecimal (i.e., binary) form for machine use. A 
program of this type is available for the Illiac. It is called 
the Decimal Order Input Routine (See Chapter 5)* 

3.2 EXECUTION OF 0RI3ERS . The Illiac operates by using 
orders which have been stored in the memory according to a plan deter- 
mined by the programmer. The program is begun with a particular 
order chosen by the programmer » Let us suppose that it is the left- 
hand order at location 10. (We shall refer to addresses in decimal 
notation.) The control will put into R the order pair from location 
10. Then, until something is said to the contrary the control will 
follow a fixed pattern in executing orders. It will do the left-hand 
order and then the right-hand order in R„. Then it will put the order 
pair at location 11 into R^ and again do the left and right-hand orders 
It will continue to withdraw and execute order pairs from successive 
memory locations until one of two things occurs: 

(a) One of the orders brought out says "stop," 

(b) One of the orders brought out says "change the 
sequencing," 
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Control Transfer Orders . The second Kind of order, a con - 
trol transfer order , permits the programmer to change the sequencing 
of orders and provides the flexibility required for iterative processes. 
It works in the following way. Let us suppose that after the machine 
has executed the left-hand order at location 17 the programmer wishes 
to move to a sequence of orders "beginning, say, with the right-hand 
order at location 35 , Then the right-hand order at location 17 will 
say "Transfer control to the right-hand order at location 25." The 
execution of this order will consist of arranging that the next order 
pair is "brought to R from location 35 and that the left-hand order 
is skipped. Having done the right-hand order in R , the control brings 
out the order pair from location 36 and proceeds in the usual way. 

Conditional Transfer Orders . There are two kinds of con- 
trol transfer orders, conditional transfer order s and unconditional 
transfer orders. We have just described the unconditional transfer 
order. The conditional transfer order does the same thing as the 
unconditional transfer order provided that the number in the accumu - 
lator is not negative . Otherwise it does nothing. In the example 
we have just used, if the right-hand order at location 17 has been 
conditional, then if the accumulator had held zero or a positive 
number the next order executed would have been the right-hand order 
at location 35. But if the accumulator had held a negative number 
the next order executed would have been the left-hand order at location 
18. 
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3,3 STOP ORDERS . There are four contro l transfer orders^ 
the rjght-'hand and left-hand unconditional and conditional orders . 
Ttnese have been further combined with a stop or n ot^stop order, giving 
eight orders in all. If the programmer chooses one of the four "stop 
and transfer control" orders the machine will stop before transferring 
control and must be started again with a panel switch. The switch can 
be set so that these stop orders are ignored, which of course makes 
these stop orders into ordinary control transfer orders. This is often 
useful in programming, and when it is done the machine can still be 
stopped with another kind of stop order which cannoU be ignored. 

3,1|. ORDER TYPES . The number of orders which the Illiac 
can execute is quite large (more than lOO), but not all of them are 
of general interest. What we shall do here is discuss the various 
types of orders, giving the variants of each. Following this dis- 
cussion will be a list of those Illiac orders which are of most use. 
This list is adequate for the coding of any problem, and beginners 
are urged to confine themselves to it. Experienced coders will find 
uses for other variants . 

The two sexadecimal instruction digits of an order give 
the order type and the order variant . We shall refer to them as the 
T- and V- digits, respectively. In the example L5 cited earlier the 
T- digit is L and denotes addition. The V- digit is 5 and chooses 
one of the variants of the addition order. 
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The order types are given in Table 3»1- 



T-Dlsit Order Type 

Left Shift 

1 Right Shift 

2 Unconditional Transfer 

3 Conditional Transfer 
I^ Store from A register 

5 Memory to Q register 

6 Divide 

7 Multiply 

8 Input or Output 

9 Special Input or Output 
K Increment Add from Q 

S Add from Q register 

N Not Used 

J Extract 

F Increment Add 

L Add 

Table 3,1 
Order Types 

The meaning of the address digits of the different order 
types is given in Table 3*2 
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ORDER- TYPE 


0, 1 


2, 3 


k 


5 


6 


7 


80, 82 



ADDRESS SIGNIFICMCE 

Number of shifts (cannot exceed 63) 

Memory location from which next order 
pair will come 

Memory location at which storage will 
occur 

Memory location from which word is 
brought to Q register 

Memory location of divisor 

Memory location of multiplicand 

Number of binary digits to be input 
or output 

92 Character to be punched and number of 

punchings 

^f S Address not used 

J Same as for 5 order 

F Memory location of addend 

L Memory location of addend 

Table 3.2 
Meaning of Address Digits 



It will be seen from Table 3,2 that except for the 0, 1, 
8_,_£^__K_ and S orders the address of an order always refers to the 
memory . 

3.5 ORDER VARIANTS . Let us now consider the variants 
obtained by changing the V-digit. The sexadecimal V-digit is made 
up of the four binary digits V8, V^^, V2 and VI. There are 16 possi- 
ble combinations, giving V-digits from to L, but not all are used 
The results obtained with different V-digits are as follows: 
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(a) The Dl^lt VI, If VI « 1, A will always be cleared to 
zero at the beginning of any order. If VI « 0^ A will not be cleared. 

^hus an odd V-digit means that A will be cleared. 

(b) The Digit VQ. If VI » 1 and V8 « 1, the quantity l/2 
will be put in A at the beginning of any order. This is how rounded 
multiplication is carried out, Jf V8 = 1 and VI - 0, the Illiac will 

hang u^. 

(c) The Digit V2. This digit affects all orders except 
those of types 0, 1, ^ and J. See Table 3 .3- 

(d) The Digit V^o .This digit affects all o rders except 
those of types 0, 1^ 5 and J. See Table 3.3. 

The following notation is convenient for a more detailed 
description of the orders i 

A = accumulator register 
a^ = sign digit of A 

Q = multiplier-quotient register 
q = sign digit of Q 

AQ = the 79 binary digit double register formed 
from A and Q by omitting q^, 

N(R) = contents of register R 

N(n) = contents of memory location n 

If no ambiguity Is possible, the symbols A and Q will also 
be used to denote the contents of A and of Q, 
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On the following pages there are detailed descriptions of 
the results obtained by changing the V-digit of the 15 order types 
used in the Illiac . After the detailed descriptions is an abbrev- 
iated list of orders. The orders in the abbreviated list have been 
underlined in the detailed descriptions. 
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TYPii; 


V^ 


EFFECT OF V^^ 


V2 


EFFECT OF V2 






0, 1 


— 


NONE 


— 


NONE 






2. 5 




1 


RIGHT-HAND ORDER 
LEFT-HAND ORDER 



1 


STOP IF SWITCH SET TO OBEY 
NONE 






il 




1 


NOIME 





STORE FULL WORD 




1 

H 
H 


RIGHT-HAND ADDRESS 
LE1*T-HAND ADDRESS 


1 


STORE ADDRESS 




5, J 




NONE 


— 


NONE 




6 


1 


V^ MUST BE 1 


1 


V2 MUST BE 1 






8 




1 


TAPE 
DRUM 




1 


INPUT 
OUTPUT 




9 





TAPE 




1 


INPUT 
OUTPUT 






T, K, S 

F, L 






1 


SUBTRACT 
ADD 




1 


NUMBER 
ABSOLUTE VALUE 





Table 3-5 Effect of Digits VU and V2 



^^ n Left Shift (Double) 3 Orders 

Final Stop 1 order 

If n = 0, the machine will stop. If not, repeat n times 
the operation which replaces the contents 



tQ. \> a^, . . . , a^Q, a^^; q^, q^, q^, . . „ , q^^, q^^ 



of AQ by 



a^, ag, ay ... , a^^ q^; q^, q^ q^, . . , , q^^, Q 
leaving q^ unchanged. 

The number n will be interpreted modulo 6h. 

Variants 

0, 2, k, 6 Order as described above. 

'^f 3f 5^ 7 Clear A and then execute as described. 

^f ^> N Illiac will hang up. Avoid these. 

9, S, J, L Clear A, insert I/2 in A, and then execute 
order as described above. 

F Final Stop. 



Use only 00, 01, Og for shifts 
Use OF with n = for final stop 
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IV n 



Right Shift (Halve) 3 Orders 



If n = 0, the machine will stop. If not, repeat n times 
the operation which replaces the contents 

s-Q, ^1' a^, » . . , a^Q, a^^; q^, q^, qy - - - > ^38' ^39 

of AQ hy 

aQ, a^, a^, , . , a^^, a^g; a^^, q^, qg. • • » ^ ^37' ^38' 

leaving q^ unchanged. 

The number n will be interpreted modulo 6^1-, 

Varients 

2^6 Order as described above, 

12 5 1 Clear A and then execute as described. 

8 K, N, F Illiac will hang up. Avoid these. 

9, S, J, L Clear A, insert l/2 in A, and then 

execute order as described above. 

Use only 10, 11, 19 
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2V n Unconditional Control Transfer 12 Orders 

Bring the next order pair from memory location n and 
choose the left or right hand order of this pair, stopping bisfore- 
hand or not, depending upon V. 

Variants 

Stop. The first order after starting with 

the START switch will be the right-hand order 
at memory location n. The stop can be ignored 
by setting a panel switch. 

2 Transfer control to right-hand order at memory 

location. 

h Same as except take left-hand order. 

6 Same as 2 except take left-hand order. 

1> 3} 5f 1 Same as 0, 2, k, 6 except clear A first. 

8, K, N. F Illiac will hang up. Avoid these. 

9, S, J, L Same as 1, 3> 5> 7 except also make A = l/2 

after clearing. 

Use 20, 21, 22, 23, ^, 25, 26, 2?, 29, 2S, 2J, 2L. 

Starting After A Stop . When the Illiac has been stopped 
by one of the control transfer stop orders, it is usually started 
again by moving the black switch to START, from which position the 
switch automatically returns to OBEY. 

The Illiac can also be started again by moving the white 
switch through EXECUTE to FETCH and then back to RUN. If this is 
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done, the control transfer order which stopped the tlliac will be 
ignored. The normal sequencing will then follow unless the stop 
order is a right hand order transferring control to the right hand 
side of a word. In this case, the order first obeyed after starting 
will be the right hand instead of the left hand order of the new 
order pair brought out . 

For example consider the following orders 3. 
p L5 F 

20 p+2 
p4l UO IF 

2h p 
p+2 7J 2F 
Lk 3F 

If we stop with the 20 order, the black switch will start with 
L^ 5F and the white switch with 2^ p. If we stop with the 2k- 
order, the black switch will start with L5 F and the white 
switch with 7J 2F. 
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3V n Conditional Control Transfer k Orders 



If A > Oj "bring the next order pair from memory location 
n and ohoose the left or right hand order of this pair, stopping be- 
forehand or not, depending upon V. If A < 0, go on to the next 
order. 

Variants 
0,2,^,6 IfA>0, do the same operation as for 

the corresponding 2V order. If A < 0, 
go on to the next order. 
Ij 3f 5i 7 Identical with corresponding 2V orders. 

8, K, N, P Illiac will hang up. Avoid these, 

9, S, J, L Identical with corresponding 2V orders . 

Use only 30, 32, 3^, 36. 

Starting After A Stop. The discussion given with the 
2V orders applies here to the corresponding 3V orders . 
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kY n Store 9 Orders 



Copy into memory location n all of the contents of A, 
the contents corresponding to the address of a left-hand order, or 
the contents corresponding to the address of a right-hand order, 
depending upon V. 

Variants 

0, h Replace N(n) by A. 

1, 5 Replace N(n) and A by 0. 

2 Replace address digits of right-hand order 
at memory location n by the corresponding 
digits of A. 

3 Same as 2 except clear A first. 

6 Same as 2 except take left-hand order. 

7 Same as 6 except clear A first. 

8, K, N, F Illiac will hang up. Avoid these. 

9, J Replace N(n) and A by l/2. 

S Replace A by l/2 and address digits of 

right-hand order at memory location n 
by 0. 

L Same as S except take left-hand order. 

Use only kO, _ia, kg, k^, h6, h^J , _^, i^S, ^L. 
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5V n Memory to Q 3 Orders 

Transfer N(n) to Q 

Variants 

0, 2., h, 6 Transfer N(n) to Q 

1, 3, 5, 7 Clear A and transfer N(n) to Q. 

8, K, N, F Illiac will hang up. Avoid these. 

9, S, J, L Put l/2 in A and transfer N(n) to Q. 

Use only ^, 51, 59 • 



3-18 



6V n 



Divide 3 Orders 



Divide AQ by N(n), placing the rounded quotient in Q 
(the least significant digit being 1 for the roundoff) and leaving 
a residue in A. If |a| > |N(n)| the Illiac will stop after 
dividing. If |a| = |N(n)| and if A > 0, the Illiac will stop 
after dividing; if |a| = |N(n)| and if A < 0, the Illiac will 
not stop after dividing. 

Variants 

6 As described above. 

7 Ivfeke A = 0, then proceed as above. 

L s Make A = l/2, then proceed as above. 

8 K, N, F Illiac will hang up. Avoid these. 

-) 1 p 5 These give incorrect results or results 
i^ 5, 9, J which are correct only under certain 

conditions. Avoid them. 



Use only GG , 67, 6l. 
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TV n Multiply 12 Orders 



/ V -59 
Put Q X P(n) + 2 -^^ A into AQ, the least significant 



59 digits being in Q with q, = 0. 



Variants 



P(n) = - N(n) 

1 P(n) = - N(n); A = "Unrounded Negative Multiply" 

2 P(n) = - |N(n)| 

5 P(n) = - |N(n)| ; A = 

^ P(n) = N(n) "Hold Multiply" 

5 P(n) = N(n) ; A = "Unrounded Multiply" 

6 P(n) = |N(n)| 

7 P(n) = |N(n)| ; A = 

8, K, N, F Ill iac vj JUL-JaaiQ ^ u^ . Avoid these. 

9 P(n) = -N(n) ; A = 1/2 "Rounded Negative Multiply" 

S P(n) = - |N(n)| ; A = 1/2 

J P(n) = N(n) ; A = l/2 "Rounded Multiply" 

L P(n) = |N(n)| ; A = l/2 

Use only 70, 71, 72, 73, £t, 75, 
76, TJ, 79, 7S, 7J, 7L 
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8v n Input-Output 9 Orders 

Transfer words "between A and the input tape, output punch, 
or magnetic drum. 

The address n must he a multiple of k for the tape and punch 
orders and must be 11 for drum orders. 

Variants 

Shift AQ four places left as in the 00 order 

and replace a^, a^„, a^g, a^Q by the binary- 
digits corresponding to the sexadecimal 
character being read. This is done nA times. 

1, 9 Clear A and then do as in 80 order, 

2 Punch the digits aQ, a , a^, a^ as one sexa- 
decimal character and shift AQ four places 
left as in the 00 order. This is done n/k 
times. 

3 Clear A and do as in 82 order. 

5 This is a ^0-digit order of the form 85 11 TV p, 

We distinguish two cases . 
(1) T is not 0, 1, 8, 9. In this case, A 
and Q are shifted left eleven places as 
in the 00 order and the word at drum 
location p is placed in A. Then the T 
order is executed using address p. Com- 
plete freedom is not available in drum 
addresses because p may interfere with V, 
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(2) T Is 0, 1, 8, 9, In this case after 
the word at drum location p is placed 
in A the right-hand order is skipped. 
This permits use of any drum address 
for p. 

6 This is a 40-digit order of the form 85 11 TV p, 

We again have two cases : 

(1) T is not 0, 1, 8, 9. In this case A is 
transferred to drum location p, and A 
and Q are shifted 11 places left as in 
the 00 order. Then the T order is ex- 
ecuted using address p. Complete freedom 
is not availa.ble in drum addresses be- 
cause p may interfere with V. 

(2) T is 0, 1, 8, 9. In this case the right- 
hand order is skipped after doing the 
left-hand order as in case (l). This 
permits use of any drum address p. 

7 Same as 86 except clear A first. 

S Put 1/2 in A and do as in 82 order. 

L Put 1/2 in A ana do as in 86 order. 

^' ^> ^y F Ill£ac_wil3^l2ang up,. — Awi^r-4fe^»«.^ 

^> J These are not useful. Avoid them. 

Use 80, 81, 82, 83, 8^, 86, 87, 8S, 8l, 
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9V n Special Input -Output 2 Orders 



Variants 

L Five hole input. Clear A, shift AQ four places right 
and replace a^, a , a g, a by the "binary digits 
corresponding to the four least significant holes on 
the tape. Place the contents of the fifth hole in 
position a . The address part of this order must be k 

2 Letter output. Punch on the tape a character de- 
pending upon the address digits n. Three quantities 
are defined by the 10 binary address digits s 
(l) The leftmost k digits define the usual k digit 

positions in the output tape. 
(2 ) The rightmost digit defines the 5th hole in 

the output tape. 
(5) The rightmost 6 digits determine the number 
b of times that the above-defined character 
is punched and also the number of right 
shifts executed. The number of characters 
punched will be found by dividing the number 
in the rightmost 6 digits by h and rounding 
u£ to the next integer . 

The address n may always be found from the following 
formula s 

n = Ska + U-b + c - 2 

where a is the character punched, a=0,l,2, . . ., J, F, L 

b is the number of characters punched, 1 < ^ < 1^ 
c determines the fifth hole, c = 0, 1. 
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The number of right shifts executed is ^b + c - 2 . 
Example . Punch the character 7 thirteen times. 

n ^ 6h X 1 + k X I'-j + - 2 = ^98 

It will be found that the last 6 digits contain the 
number 50 which when divided by h and rounded up gives 15 . There 
will be 50 right shifts . 

Figure 3.2 shows tlie relationship between the tape 
holes and the address digits of the 92 order. The address shown 
will print the character 7 thirteen times . 




Figure 3 '2 
Address Digits of 92 Order 
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KV n Increment Add from Q 12 Orders 



These orders are identical with the FV orders except 
that Q is used instead of N(n). 

Variants 

Add -Q - 2"^^ to A. (See note). 

1 Put -Q - 2~^^ in A. 

2 Same as if Q > 0; same as h if Q < 0, 

3 Same as 1 if Q > 0; same as 5 if Q < 0. 
h Add Q + 2 -^^ to A. 

5 Put Q + 2"^^ to A. 

6 Same as ^ if Q, > Oj same as if Q < , 

7 Same as 5 if Q > 0; same as 1 if Q < 0^ 
8, K, N, F jn_i-igo yjli ^^^P, ^T — -A^^iid-Jili^se . 

9 Put -Q - 2~^^ + 1/2 in A. 

S Same as 9 if Q > 0; same as J if Q < 

J Put Q + 2"^^ + 1/2 in A. 

L Same as J if Q > 0; same as 9 if Q < 



Use KO, Kl, K2, K3, K4, K5, K6, K7, K9, KS, KJ, KL 



-59 
NOTE: -Q - 2 is the digitwise complement of Q. 
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SV n Add from Q 

Q is added, subtracted, etc., to A. These orders are 
identical with the LV orders vith Q used instead of N(n). 

Variants 



Subtract Q from A. 

1 Put -q in A. 

2 Subtract |q| from A 
5 Put - |Q| in a. 

^ Add Q to A. 

5 Put Q in A. 

6 Add |q| to A. 

7 Put |q| in A. 

^' ^' ^' ^ £l-liac_will hang; up. Avoid thf^ ^. 

9 Put 1/2 - Q in A. 

S Put 1/2 - |q| in A. 

J Put 1/2 + Q in A. 

L Put 1/2 + |Q| in A. 



Use SO, SI, S2, S5, SU, S5, S6, S7, S9, SS, S.T, SL 



12 Orders 
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jV n Extract 5 Orders 



If two corresponding digits of N(n) and Q are both 
I's, put 1 in that place of Q. Otherwise put 0. This order gives 
the logical product of N(n) and Q. 

Variants 

0, 2) h, 6 As described above. 

1, 5, 5, 7 Clear A and do JO order. 

8, K, N, F Illiac will hang up. Avoid these . 

9, S, J, L Put 1/2 in A and do JO order. 

Use only JO, Jl, J9 
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FV n Increment Add 12 Orders 

Programmed Stop 1 Order 



Variants 

Add -N(n) - 2"^^ to A. (See note) 

1 Put -N(n) - 2'^^ in A, 

2 Same as if N(n) > 0; same as k if 

N(n) < 0. 

3 Same as 1 if N(n) > 0; same as 5 if 

N(n) < 0. 

h Add N(n) + 2"^^ in A. 

5 Put N(n) + 2*^^ in A. 

6 Same as h if N(n) > 0; same as if 

N(n) < 0. 

7 Same as 5 if N(n) > 0; same as 1 if 

N(n) < 0. 

8, K, N Ill iac vill hang up. Avoid these. 

9 Put -N(n) - 2"^^ + 1/2 in A. 

S Same as 9 if N(n) > 0; same as J if 

N(n) < 0. 
J Put N(n) + 2"^^ + 1/2 in A. 

L Same as J if N(n) > 0; same as 9 if 

N(n) < 0. 

F Stop. This stop is used to indicate 

failure of a programmed check. 



Use FO, Fl, F2, F3, F^^, F^, F6, F^J, F9, FS, FJ, FT, FL. 
NOTE: -N(n) - 2"^^ is the digitwise complement of N(,n). 
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LV n Add 12 Orders 



N(n) is added, subtracted, etc., to A. 

Variants 

Subtract N(n) from A. 

1 Put -N(n) in A. 

2 Subtract |N(n)| from A. 
5 Put - |N(n)| in A. 

k Add N(n) to A. 

5 Put N(n) in A. 

6 Add |N(n)| to A. 

7 Put |N(n)| in A. 

8, K, N, F Illiac will hang up. Avoid these. 

9 Put l/2 - N(n) in A. 

S Put 1/2 - |N(n)| in A. 

J Put 1/2 + N(n) in A. 

L Put 1/2 a |N(n)| in A. 

Use LO, LI, 12, L3, L^, L5, L6, L?, L9> LS, LJ, LL^ 
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ABBREVIATED ORDER LIST 



ORDER 



00 


n 


09 


n 


OF 





10 


n 


19 


n 



20 n 



22 n 

2k n 
26 n 



DESCRIPTION 

Shift AQ left n places, 1 < n < 63. 

Make A = I/2 and then shift AQ left n 
places, 1 < n < 63. 

Final stop. 

Shift AQ right n places, 1 < n < 63 . 

Make A = I/2 and shift AQ right n places 

„ -1 
so that AQ contains 2*" " , l<n<65. 

Stop. The first order after START will 
be the right hand order at location n. 
The stop can be ignored by setting the 
panel switch to IGNORE. 

Transfer control to the right hand order 
at location n. 

Same as 20 except take left hand order. 

Same as 22 except take left hand order. 



30 


n 


32 


n 


3k 


n 


56 


n 


40 


n 


41 


n 



If A > do as in the corresponding 2V 
order . 

If A < go on to the next order. 



Replace N(n) by A. Do not change A 
Replace N(n) and A by zero. 
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ORDER DESCRIPTION 

k2 n Replace address digits of the right hand 

order at location n by the corresponding 
digits of A. Do not change A. 

k6 n Replace address digits of the left hand 

order at location n "by the corresponding 
digits of A. Do not change A. 

ij-9 n Replace N(n) and A by l/2 . 

50 n Replace Q by N(n). 

66 n Divide A + 2"^^ (Q + q^ ) by N(n). 

70 n Put -N(n) Q + 2"^^ A into AQ. 

71 n Multiply -N(n) by Q, putting result in AQ. 

72 n %Put - |N(n)| Q + 2 -^^ A into AQ. 

73 n Multiply - |N(n)| by Q, putting result in AQ. 
^h n Put N(n) Q + 2~^^ A into AQ. 

75 n Multiply N(n) by Q, putting result in AQ. 

76 n Put |N(n)| Q + 2 -^^ A into AQ. 

77 n Multiply |N(n)| by Q, putting result in AQ. 

79 n Put the rounded product -N(n) Q into A. 

7S n Put the rounded product - |N(n)| Q into A. 

7J n Put the rounded product N(n) Q into A. 

7L n Put the rounded product |N(n)| Q into A. 

80 n Input n/k sexadecimal characters from the tape 

81 n Clear A and input n/^ sexadecimal characters 

from the tape . 
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252M DESCRIPTION 

^^ ^ P^ch n/h sexadecimal characters on the tape 

^5 " Repla.ce A by the word at drum location n. 

See page 3-21. 

^^ " Replace the word at drum location n by A. 

See page 3-22 . 

'^^ ^ ^^ad one five-hole character from the tape. 

92 131 Punch one carriage return and line feed 
character . 

92 ^^2 Punch one + sign (or K). 

92 706 Punch one - sign (or S). 

92 769 Advance film. 

92 963 Punch one space character. 

^ "^ Put - Q - 2"^^ in A. (This is the digitwise 

complement of Q). 

Add Q + 2"^^ to A. 

Put Q + 2"^^ in A. 

Subtract Q from A. 
Put - Q in A. 
Subtract |q| from A. 
Put -|q| in a. 
Add Q to A. 
Put Q in A. 
Add |q| to a. 
Put |q| in A. 
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Ki^ 


n 


K5 


n 


SO 


n 


SI 


n 


S2 


n 


S3 


n 


Sk 


n 


S5 


n 


S6 


n 


S7 


n 



ORDER DESCRIPTION 

JO n Replace Q by the logical product of Q and 

N(n). 

Fl n Put -N(n) - 2"^^ in A. (This is the digit- 

wise complement of N(n). 

Add N(n) + 2"^^ to A. 

Put N(n) + 2"^^ in A. 

Subtract N(n) from A. 
Put -N(n) in A. 
Subtract |N(n)| from A. 
Put - |N(n) I in A. 
Add N(n) to A. 
Put N(n) in A. 
Add |N(n)| to A. 
Put |N(n)| in A. 



Fk 


n 


F5 


n 


LO 


n 


LI 


n 


L2 


n 


L5 


n 


L4 


n 


L5 


n 


l6 


n 


L7 


n 
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CHAPTER h 



ROUTINES 



One of the first tasks in programming a computation is 
to break it down into a number of small or medium size operations, 
each of which forms a fairly distinct logical step in the whole 
calculation. Such a step may be the evaluation of a function such 
as a square root, logarithm or cosine. Each step must be carefully 
defined so that it fits correctly with all the other steps. Having 
decided exactly what is required of each step, it is possible to 
proceed with the coding of the individual steps. 

This procedure has three important advantages. Firstly, 
it enables the programmer to concentrate on one part of the job at 
a time. Secondly, the coding of each step can be tested separately 
before incorporating it in the program. Thirdly, certain steps 
are common to a very large number of different calculations; these 
have been coded and tested once for all, and may be used by anyone 
in the Laboratory, 

The coding for one step constitutes a routine . The 
collection of routines for performing common operations in the Illiac 
is called the library of routines. It includes routines for evaluating 
many simple functions like those mentioned above, also for printing 
numbers in various layouts, for integrating differential equations, 
for solving simultaneous equations, and for many other operations. 
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A routine, when used to perform part of the work of another 
routine, is called a subroutine . The simplest type of subroutine 
is merely a string of order pairs which can be inserted in the appro- 
priate place among the other orders of a program. This is called an 
OEen subroutine. However, for various reasons another type, called 
a closed subroutine, is more commonly used. 

^.1 CLOSED SUBROUTINES' A closed subroutine is also 
a string of consecutive words but these do not have to be placed 
among the other orders of the program. Instead, they may be put 
in any convenient place in the store. Each time the subroutine 
is to be executed, control must be transferred to it in a certain 
special way (the subroutine is said to be entered ). It is so arranged 
that when execution is complete, control is automatically returned to 
the point from which the subroutine was entered, so that execution of 
the rest of the program may continue. 

In this way a program is seen to consist of several dis- 
tinct, self-contained blocks, namely the various subroutines and 
the part of the program (usually called the main program or master 
routine) which makes use of its subroutines by sending control to 
them. Internal rearrangement of a routine is usually difficult, 
but the routines comprising a program can be shuffled about very 
easily, and this fact makes the coding of large problems much simpler. 

It is not necessary to enter each subroutine directly from 
the master routine; there is nothing to prevent one subroutine being 
entered from another. A subroutine may itself have subroutines. 



4-2 



k.2 ENTERING A CLOSED SUBROUTIN E. The following two 
orders must "be inserted in the master routine at the point from 
which a subroutine is required to be entered. 

k Apiy 

k+1 26 y Subroutine starts at m. 
Any 

These orders must be, as shown, in the right-hand half of one word 
and the left-hand half of the next (it may be necessary to waste an 
order to do this). The 50 order contains the address of the memory 
location in which it is itself contained. The 26 order contains the 
address of the memory location containing the first word of the sub- 
routine (it actually transfers control to the left-hand half of this 

word ) . 

The subroutine, after execution, automatically transfers 
control to the order following these, i.e., to the right-hand order 
in memory location (k+l), so that there is no break in the logical 
continuity of the master routine. There is no need for the user to 
know how the subroutine does this, but it is not difficult to under- 
stand and it helps to complete the picture. 

if. 3 RETURNING CONTROL TO THE MASTER ROUTINE . It will 
be seen that the effect of the 50 order above is to place the order 
pair, of which it is the right-hand half, in the Q register. Thus 
whenever a closed subroutine is entered, the Q register contains a 
pair of orders of which the right-hand order contains the address 
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of the memory location from which the pair came. Control must be 
returned to the right-hand order in the memory location following 
this. 

The subroutine uses the information in the Q register to 
set the address in a transfer of control order (a 22 order) called 
the link , which ultimately causes the transfer back to the master 
routine. It does this by means of the orders shown in Table k.l. 



m 



m+1 



K5 

k2 




m+n 



These orders "plant" 
the link 



m+n 



22 (Ojj 



link 



Table h.l 
Forming a Closed Subroutine Link 

The K5 order transfers the order pair from the Q register 
to the A register, increasing the right-hand address by 1 as it 
does so. This address is now the address to which control is 
eventually to be returned. The k2 order places this address in 
the link. The link is situated so that it will be encountered 
by control when the execution of the subroutine is complete. Note 
that the link should be a right-hand order so that its address may 
be inserted from the right-hand side of the A register. 

It may sometimes be convenient to use a 32 order as a link; 
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this order may then perform a test within the subroutine and return 
control to the master routine only if some condition is satisfied. 
Also the above procedure for "planting" the link may be varied if 
desired, so long as the result is the same. In particular note the 
variants in sections U.5 and k.6 below. 

k.h PLACING THE ARGUMENT . All subroutines operate on 
at least one number somewhere in the machine, and there must be 
agreement between the subroutine and the master routine on the 
placing of these numbers, and also on the placing of the results 
of the subroutine operation. 

If only one number is operated upon, it is convenient 
to use the A register to hold this number when the subroutine is 
entered (since the A register is not being used for anything else 
at this moment). Similarly, if only one number is produced by the 
subroutine it is convenient for the subroutine to leave this number 

in the A register. 

For example. Library Routine R 1 is a closed subroutine that 
finds the square root of the number given in A, and leaves this 
square root in A. Suppose we wish the Illiac to find the largest 

root of 

i|. 3 2 
X +2x +x -c=0 

where c = N(10) and satisifies < c < l/2 . It can be shown that 

the root is 

X = -1/2 + /l/^ + v^ • 
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Suppose that the constants l/2 and l/i^- are given in 20 and 21 

respectively, and that x is required in 11. Routine R 1 will be 

used to find the square roots; let its first word be in 100. Then 

the master routine might contain the orders of Table k.2 

Put c in A 

Enter Code R 1 

to form /c~ 

Add l/k to /-C" 

Waste order 

Enter Code R 1 

to form ^ 1/h + \r~c 

Subtract l/2 to form x 

Put X in 11 

Table k .2 
Master Routine Using A Closed Subroutine 

The subroutine itself must use a slightly different method 
or link planting from that given in section h.3, which would des- 
troy the argument in A before using it. It is here necessary to 
rescue the argument and hold it in a storage location until the 
link has been planted, as in Table k.^ 

m ^0 q Put argument in q 

K5 
m+1 k2 m+n Plant link 



50 


L5 


10 




50 


50 


51 


26 


100 




L4 


21 


52 


22 


52 




50 


52 


53 


26 


100 




LO 


20 


5^ 


40 


11 



m+n 



22 (0) Link 

Table 4.3 
Modified Link Planting 
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^•5 PROGRAM PARAMETERS . Sometimes a subroutine is 
made to carry out slightly different operations on different occasions 

For example, one subroutine might be made to print numbers to any 

th 
number of figures. Another might form the n root of a number, 

where n is any positive integer. The value of n, or the number of 

figures to be printed, is called a parameter of the subroutine. 

The value of a parameter must always be specified if the 
subroutine is to operate correctly. There are two standard ways of 
specifying parameters; parameters specified in the way now to be 
described are called program parameters . (The others are called 
"preset" parameters - see below, section ^.T)* 

A program parameter is specified in the first half of 
the word containing the 50 order which is required on entering any 
closed subroutine. Let us take for example Library Routine R 2, 
which is a closed subroutine for replacing N(A) by its p root. 
Here p is a program parameter, to be specified as follows: 

k 50 p Program parameter 

50 k 
k+1 26 m Subroutine starts at m 

Any 

To illustrate this, suppose we wish the machine to com- 
pute the positive real root of 

x^ - 5x^ + 3x^ + 2c = 0, (1/2 < c < 1) 
which is given by the formula 

v: 



X = + ^ 1 + ~\'/ 1 - 2c 
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Suppose that c = N(10) and that x is reqxiired in 11; suppose also 
that N(20) =-■ -1, and that Routine R 2 starts at 100. Then the master 
routine might contain the program given in Table k.k 

N(A) = (1/2) - c 

N(A) = 1 - 2c 

Program parameter 

Link 

To Routine R 2 



50 


L9 


10 




00 


1 


51 


50 


3 




50 


51 


52 


26 


100 




LO 


20 


55 


50 


2 




50 


55 


^^ 


26 


100 




ho 
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N(A) = 1 + "V 1 - 2c 

Program parameter 

Link 

To Routine R 2 

X to 11 

Table k,k- 
Use of Program Parameter 

It will be noticed in Table k-,k- that the order contain- 
ing the parameter (50 3 or 50 2 ) is actually obeyed by the machine 
before the subroutine is entered . This is Just a waste of time 
which is not worth avoiding; it does no good and, provided the order 
containing the parameter is suitably chosen, no harm. It is usual 
to use a 50 order here; this puts rubbish in the Q register, but 
the latter is then immediately reset by the following 50 order and 
no damage has been done . 

A more elaborate example of the use of a program parameter 
is to be found in Library Routine P 1 which prints A as an integer 
or fraction. The entry is by means of the orders of Table ^.5. 
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k XY d 

50 k 

k4l 26 m 

Table k.^ 
Multiple Program Parameters 

In Table U.5 x = 5 or J; if x = 5 and A is negative the printed 
number is preceded by a minus sign, otherwise it is preceded by 
a space. Also Y = or 2; if Y = 0, A x 2-"^ is printed as an 
integer; if Y = 2^ A is printed as a fraction (correctly rounded 
off). The layout of digits is given by d = lOq + s; q digits are 
printed with a space after the first s. 

It will be noted that all combinations of X and Y form 
harmless orders. , y .A ', 

h,6 IITTERPRETIVE ROUTINES .' There is a type of sub- 
routine which, instead of executing a single distinctive operation, 
carries out a whole series of operations. Each operation requires 
a parameter for its specification, so that the master routine contains 
a string of parameters, one for each operation. The string may be 
of indefinite length. 

Such subroutines are called interpretive routines. Their 
use lies in programs that involve actual operations on elements which 
are not numbers stored in the usual form but may be numbers stored 
in some special form, or different mathematical entities altogether 
such as expansions in Boolean algebra. The commonest application 
is to numbers stored in the so-called "floating point" form (see 
Chapter 6). There are certain advantages in storing numbers in this 
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form; however, such numbers cannot be added or multiplied in a single 
Illiac operation. A routine is required to handle simple arithmetic, 
and for this purpose an interpretive routine is used (e.g. Routine Al). 
Each parameter of the interpretive routine corresponds to one arith- 
metical operation, just as each Illiac order corresponds to one opera- 
tion in the Illiac. Thus, for example, to place in register 6 the 
sum of the numbers in registers k and 5 we write the parameters 

8i^ 5 
8S 6, 

which act in a similar way to the Illiac orders 

L5 ^ 
Lk 5 
ho 6. 

Owing to the close similarity between the parameters and 
the ordinary Illiac orders, the parameters are themselves often re- 
ferred to as interpretive orders, or merely as orders, "orders," or 
orders . The "order code" of the interpretive routine may be des- 
cribed in a similar way to the order code of the machine; it involves 
reference to an "accumulator" which behaves, for floating point num- 
bers, like the A register of the Illiac. 

An interpretive routine is entered in the same way as a 
closed subroutine, but the parameters (i.e., the interpretive "orders") 
follow the orders causing entry. Finally a special parameter may 
be used to cause control to the Illiac to be transferred out of the 
interpretive routine and back to the master routine. Thus if we 
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imagine that the above example is a complete set of operations to 
be carried out by the interpretive routine, the master routine 
would contain the program of Table 4.6, 

k Any 

50 k Enter interpretive routine 

k+1 26 m 

85 h Parameters (interpretive orders) 

k+2 Qk 5 

8S 6 

k+3 8 J kf4 Send control of Illiac to L.H. 

00 '^^^ °^ ''+'' 

etc . 

Table k.6 
Program for Interpretive Routine 

k.J PRESET PARAMETERS , By making use of a li'brary routine 
a programmer avoids not only the necessity of writing out the orders, 
but also the labor of punching them, since a master copy of the tape 
is kept on file in the Teletype Room to be copied when required. 

If the routine has one or more parameters associated with 
it, it can meet a wide variety of requirements. However, the use 
of program parameters consumes both computing time and storage space. 

A program parameter has the property that it can be varied 
from one application to another within the same program. Frequently 
it happens that, although the ability to choose a value of a para- 
meter to suit any particular program is desired, the ability to change 
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the value during the execution of a program is not needed. The 
value can be set before execution begins; there is no need for the 
orders which, in the case of a program parameter, set the value 
afresh each time the routine is entered. A parameter whose value 
is set before execution begins is called a preset parameter. 

The master copy of a library routine must be valid for 
all values of any preset parameters involved. The fixing of a pre- 
set parameter for a particular program must therefore be done after 
(or during) the copying of the master tape; in fact it is done as 
the program is read into the Illiac. This operation is carried 
out by the Decimal Order Input Routine which is described in 
Chapter 5» 
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CHAPTER 5 
THE DECIMAL ORDER INPUT 

The Illiac is a binary machine intended for University 
research. Since it will be used by a large group of people with 
various scientific backgrounds, it is important that its use be 
made as simple as possible. To this end a library of subroutines 
has been organized and aids in using it have been devised. One of 
these aids is the Decimal Order Input, Library Code XI. Tliis code 
has two principal purposes: 

(1) To make it possible for programmers 
to use decimal notation in coding, 

(2) To make use of the library easy. 

The use of subroutines is discussed more fully in Chapter 
h, but one way to characterize a subroutine is to say that it is 
essentially an extension of the order code of the machine. It is 
a group of orders used to carry out one or more operations which 
may be as simple as a square root or as complicated as a complete 
program for executing the details of floating point arithmetic. 
In any case the interior structure of a subroutine will depend upon 
its location in the memory. For example, the left-hand address of 
the second word in the square root routine (Library Code Rl) refers 
to the eighth word of the routine. If the routine begins at loca- 
tion 10, this address must be 18. But if the routine begins at 
location 97, this address must be 105^ The use of subroutines is 
very awkward unless problems of address changing can be easily 
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handled . 

Another principal nuisance in coding is the number system 
used by the Illiac. In the sexadecimal system using Illiac notation 
(see Chapter 2), the 1024 addresses of the memory are represented 
by numbers between and 3LL. It is much more convenient if addresses 
can be expressed in decimal form with the Illiac doing the necessary 
converting to the binary system. 

With the Decimal Order Input every number corresponding 
to the address part of an order must be written in the decimal system. 
These decimal numbers are always converted to sexadecimal (i.e., binary) 
form by the Decimal Order Input* 

The Decimal Order Input overcomes the obstacles of changing 
addresses in subroutines and of converting from decimal to binary 
form. It also provides certain other useful services. 

5.1 REIATIVE Am) FIXED ADDRESSES . If we can devise some 
way of using addresses in decimal form, we never need the 6 characters 
K, S, N, J, F, L in an address and can use them for other purposes. 
Let us consider a program beginning in location 11 of the memory, as 
s-hown in Table 5<.l. 

11 40 1 
S5 1 

12 lA 13 
k6 20 

13 51 1 
10 1 

Table 5-1 
Program Beginning at Location 11 
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The same program, if begun at location 20, would read aa in Table 
5.2. 

20 kO 1 
S5 1 

21 Lk 22 
kS 29 

22 51 1 
10 1 

Table 5.2 
Program Beginning at Location 20 

The second vord .has changed. But if we mark each address which de- 
pends upon the location of the routine with the symbol L and each 
address which is independent of the location of the routine with 
the symbol F, we have the program in Table 5*3 • 

40 IF 
S5 IT 

1 lA 2L 
k6 9L 

2 51 ^ 
10 IF 

Table 5-3 
Program With Relative and Fixed Addresses 

The addresses ending in L are relative to the location of the first 
vord of the routine. 



5-3 



We can now write each subroutine as if it begins at location 
zero provided we add to each L- terminated address the location of the 
first word of the routine when we store the routine. We shall do this 
by introducing a new syinbol K, 

5.2 DIRECTIVES. Orders which are written for the Illiac 
use two sexadecimal function digits followed by an address. We 
follow this convention even in the case of pseudo orders. A directive 
is such a pseudo order. When we want a group of orders to be placed 
in memory locations n, n-fl, n+2, . , , , we punch on the tape the 
directive 

00 nK ^^ , 
The Decimal Order Input will recognize this as a directive and will 
place the orders following the directive in pairs in the locations 
nj, n+1, n+2, . „ . ^ starting with the left-hand side of location n. 
It will add n to the address of any order terminated with L before 
placing it in posit ion o The directive 00 nK is not placed in the 
memory. Thus if we have 00 IIK followed by the orders of Table 5,3 
we will place in the memory the code given in Table 5-1? while the 
code of Table 5.2 would be obtained by using the directive 00 20K. 

5.3 ASSEMB LING OF ORDERS . Since each address ends with 
an alphabet character and since there are always two function digits, 
the Decimal Order Input can distinguish between address digits and 
function digits, Ths Decimal Order Input uses the fixed storage 



We shall use small letters to represent decimal quantities and 
capital letters to represent sexadecimal quantities » 



5-^ 



locations 0, 1, and 2 as temporary storage. After two function digits 

have been read and shifted to a right-hand order position in location 

-39 
1 the routine converts the decimal address n to binary form n x 2 

and adds it to the function digits. The order pair in location 1 is 

then stored. Next the contents of 1 are shifted left 20 places and 

the next order is formed as before . Again the order pair in location 1 

is stored in the same address as the previous time , but now the correctly 

assembled order pair has been stored. The address of the store order 

which is assembling the program is increased every other time . Thus 

the left-hand 20 digits of location 1 always contain the previous 

order while the next 8 digits of location 1 contain the function digits 

of the order being put in. 

When the address ends in K (a directive) it Is stored in 
location 2 and is added to each L terminated address in location 1 
before the appropriate order is stored. 

^.h DECIMAL ADDRESSES . The address is converted one 
digit at a time as it is read from the tape. This kind of con- 
version makes it unnecessary to write non-significant zeros, so 
orders may have the form L5 7L, ^0 1021F, 26 L, etc. Indeed, 
we are not restricted to addresses smaller than 102^ and may use 
anj^-thing we please. If the address exceeds 4095 It will add to 
the function digits (remember that there are two unused digits 
between the function digits and the address), but even this may 
be made use of as we shall see below in Section 5.6. 
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5.5 STARTING THE PROGRAM . After the Decimal Order 
Input has placed a program in the memory we must start the pro- 
gram. This is done by using the terminating symbol N. The symbol 
N causes the order which it follows to be obeyed. It is used with 
a control transfer order which must appear on the tape as if it were 
a left-hand order . It will never be stored in the memory, other 
than as the right-hand order in location 1. 

For example, the order 26 pN will be obeyed and will trans- 
fer control to the left-hand order at address p. Any unconditional 
transfer order may be used. The order 20 qN will stop the computer, 
and when the computer is started again control will go to the right- 
hand order at address q. The addresses p and q are fixed addresses. 

Control can also be transferred using a relative address, 
for the previous order with its adjusted address is always in the 
left-hand side of location 1. Therefore the two orders 22 rL 
26 IN will cause control to be transferred to address r relative 
to the last directive if 22 rL is the right-hand order of the last 
order pair, i.e., if the phase is correct. 

5-6 INPUT OF DECIMAL FRACTIONS . In section 5-^ it was 

pointed out that an address is formed and added to the function 

39 
digits. This address could be as large as 2 - 1 without getting 

into the sign digit, and any positive integer n could be input as 

-39 
n X 2 by letting the left-hand order and the right-hand function 

digits be zero. For example, the "order pair" OOF OO2896P would 

-39 
appear as 2896 x 2 . Hence a 12 digit right-hand address smaller 
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than 2 (about 5.5 x 10 ) could then he converted to a decimal 

^9 / 12 
fraction by multiplying it by 2-^^/10 . 

This is done when the terminating symbol J is used. For 
example, the characters OOF OO2969 0000 OOOOJ would cause the 
quantity O.2969 to be placed in the memory. Zeros at the end can- 
not be omitted, but preceding zeros can be omitted. Do not omit 
function digits . Remember that this is essentially an integer 
input. The range can be extended by using theleft-hand function 
digits, for they are simply a number to which the right-hand 
address is added after being multiplied by 2-^^/10 . The left- 
hand function digits i+0, 80, NO reprecent I/2, -1, -l/2. Thus 
numbers could be input as follows : 

-0.8888 8888 8888 as 80F 00 1111 1111 1112 J, 
0.785^1 3216 0000 as UOF 00 285^ 32l6 OOOOJ. 
This is not an efficient way to read decimal fractions into the 
machine, but it is convenient for occasional numbers scattered 
through the program. 

5.7 PRE- SET PARAMETERS . MODIFICATION OF ORDERS . The 
remaining unused alphabetic sexadecimal character is S. This is 
used to modify orders by using pre-set parameters. The symbol S 
differs from the other terminating symbols in that it is always 
followed by another single character which may be any of the 13 
sexadecimal characters 3 to L. The termination SD on an address 
causes the content of location D = 3, ^, . . . , F, L, to be 
added to the order while it is in a right-hand position as 
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described in Section 5*3 s^nd "before it is placed in the memory. 
This is for either right or left-hand orders of completely assembled 
order pairs. For example, if location J contains ^ x 2 , the 
order pair L5 20SJ ^0 3OSJ will be modified to read L5 27F hO 
37F before being stored. This facility is very convenient for using 
parameters with a program because the program can be written with 
orders of the form L5 S3, L^ S^, and if locations 3 and h have 
been previously set the parameters will be added appropriately as 
the program is read into the machine. 

Many examples of pre-set parameters can be found in the 
library. In Illinois Code P6, Single Column Print, the parameter 
OOF OOmF is used to specify the printing of m decimal digits and 
must be placed in location 3 before Code P6 is read. Thus if we 
wish to place Code P6 in locations beginning with 50 and to print 
7 decimal places, the pertinent part of the program tape would 
read 

00 3K Directive 

00 F Place 7 x 2"^^ in 

00 7F location 3 

00 5OK Directive 

Code P6 Place Code P6 beginning 
with location 50. 

5.8 EXAMPLE OF USE OF DECIMAL ORDER INPUT . The decimal 

Order Input has 25 words and occupies locations 999 to 1023 in the 

memiory. It is placed in these locations with its own bootstrap 

input which then transfers control to the Decimal Order Input 
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so that it can take over the control of program input. Let us con- 
sider the following simple example; 

Compute to 10 decimal places the square roots of 
jc/lO and e/lO using the Square Root Routine (Code 
Rl)and the Single Column Print routine (Code P6). 
We have the following data: 

(a) Code Rl - 9 words, closed, finds 
square root of argument placed in 
A and places answer in A. 

(b) Code P6 - 1^ words, closed, prints 
words in A to ra decimal places followed 

by carriage return and line-feed. Con- 

-39 
tents of address 3 must contain m x 2 

as Code P6 is input. 

In this program, the code for which is given in Ta.ble 

5.)+ on the following page, we have scattered the orders through 

the memory to indicate how directives are used. Notice that the 

arrangement on the tape is arbitrary after the Decimal Order Input 

except that the parameter in address 3 must be in place when Code 

P6 is being input. The parameter is used in word 6 of Code P6, 

this word being 

19 63S3 
50 F 

Thus the address of the 19 order was set to 9 (73 = 9 mod 6h) for 

counting the number of digits to be printed. 

Notice also that the directive 00 56OK can be changed 
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MEMORY 


LOCATION 


PROGRAM TAPE 


REMARKS 


999 - 


1023 


Decimal Order Input 


Routine XI 






00 lOK 


Directive 


10 - 


18 


Square Root Routine 


Routine Rl 






00 3K 


Directive 


3 




00 F 
00 lOF 


Parameter 






00 30K 


Directive 


30 - 


k3 


Print Routine 


Routine P6 






00 50K 


Directive 


50 




00 F 

00 314159265359J 


7f/lO 


51 




00 F 

00 271828182846J 


e/lO 






00 56OK 


Directive 


560 




L5 5OF 


jt/lO to A 






50 L 


Link 


561 




26 lOF 


To Routine Rl 






50 IL 


Link 


562 




26 3OF 


To Routine P6 






L5 5IF 


e/lO to A 


563 




22 3L 


Waste Order 






50 3L 


Link 


564 




26 lOF 


To Routine Rl 






50 kL 


Link 


565 




26 3OF 


To Routine P6 






OF F 


Stop 


566 




22 6L 


Waste 






26 L 


Start Program at 
relative locatio: 






2k IN 


after stop 



Table 3.k 
Use of Decimal Order Input 
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so that the wnrds following it are placed differently hvt that no 
other change need be made to move these worac, because the program 
was started using a relative address. The waste order at location 
'^66 is required so that 2k IN will have a left-hand location 

5-9 USE WITH INTERLUDES . RETAINED DIRECTIVE . An 
interlude is a computation performed during the input of a program, 
the input being interrupted and then resumed. 

A tape bearing a library routine may begin with an inter- 
lude which is. placed in locations destined eventually to hold the 
routine itself. When the words of the interlude have been read 
control is directed to it, the interlude is executed* and then in- 
put is resumed. The next part of the tape carries the routine it- 
self which is written over the interlude. The purpose of the inter- 
lude is usually to prepare some orders or constants required for 
the routine. For example, a printing routine, where the number of 
digits printed is determined by a preset parameter, may use an 
interlude to compute the roundoff constant. 

Input is resumed after the interlude by transferring con- 
trol to the left side of location 999 (3F7 sexadecimal). Either 
the first word on the tape must contain the needed directive 

-"50 

00 mK or Q must contain m x 2 . 

If, upon resuming inpvit, it is desired to retain the last 
used directive , control should be transferred to the right side of 
lOlii- with m X 2''-^^ in A. The next words on the tape will be placed 
in m, m+1, . . . , retaining the previous relative address. The 
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Interlude must not use location 2 . 

5.10 STOPPING THE TAPE . The order 20 1019N on the 
tape will stop the computer and will have no other effect. Upon 
being started the tape will continue being read from where it stopped. 

5.11 PLACIKG THE DECIMAL ORDER IMFUT . BOOTSTRAPS . 

Up to this point we have not said how the Decimal Order Input is 
itself put into the Illiac. It occupies locations 999 to 1023, 
the last 25 positions of the electrostatic memory, and it is placed 
in the memory with a bootstrap input routine . 

With panel switches we can clear the control counter to zero and 
place the sexadecimal equivalent of the order pair 80 40F 40 F in the 
order register. The code in Table 5.4, which must be written with sexadeci- 
mal addresses, will then place the Decimal Order Input in the memory. 

80 028 
i^O 001 

Bo 028 
40 002 

19 026 
26 000 

80 028 V 

4o(ooo) "^^ 

lA 001 
40 001 

80 028 
h0(3F6) 

Table 5.4 
Tape for Bootstrap Input Routine 



2 



Parentheses are often placed around addresses which change during 
the course of a program. 
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This bootstrap actually places the 3-word routine (shown in both 
sexadecimal and decimal forms) of Table 5 • 5 in the memory and starts 
it 

Lk 001 LU if 



kO 001 ^0 IF 

80 028 Bo ^OF 

1 1 

ho (3F7) ^0 (999)F 
19 026 19 38F 

2 2 

26 000 26 F 

Table 5.5 
Memory Contents for Bootstrap Input Routine 

at location 1. Clearly it will take the next words on the tape 
and start putting them at location 999F. To stop the input we 
place the order pair 22 3LS 00 001 (This is 22 IOI9F 00 IF) 
on the tape so that it gets put into location 0. The control 
will be transferred to the right-hand side of location IOI9F and 
the Decimal Order Input will be started. 

This bootstrap input may be used with any code and it 
or something like it must be used to input programs whenever the 
Decimal Order Input has been overwritten. 

The terra bootstrap start is often used for tapes which 
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are started by setting the order register to 80 ^OF ^0 F and 
the control counter to zero. 
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CHAPTER 6 



SCALING 



With the convention adopted for the Illiac (See Chapter 2), 
only numbers which lie in the range -1 < x < 1 can he held in the 
registers. Since most problems require numbers outside this range, 
some scaling process is usually needed to fit a problem to the 
machine . It is necessary that each number at every stage of a 
calculation lie within the capacity of the machine. The organization 
required to assure this is sometimes trivial, but in many instances 
it is the very essence of the problem. 

6.1 SCALING BY SHIFTING . Although the number 2 lies out- 
side the range of Illiac numbers, we can multiply and divide numbers 
by powers of 2 by shifting. Thus the left shift order 00 lOF will 
c^use AQ to be multiplied by 2^^ = 1021+ . Similarly, the right shift 
order 10 9F will divide AQ by 2^ = 512 • A knowledge of the use of 
the shift orders is essential to an understanding of scaling. 

6.2 NU^^^ Wira TIffi BIMRY FOOT Sm^ Let us con- 
sider the problem of computing with numbers in which the binary point 

has been moved 10 places to the right of its Illiac position. We 

-29 
are then dealing with numbers in the range - 102^+ < y < 102^+ - 2 

Let these numbers be designated ^V .\q' Then we have H^^im) = 2 N(m) 

We can formulate rules for doing arithmetic with the numbers 

N . Addition and subtraction are simple. If N(q) = N(m) + N(n), 
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then 



2^^ N(q) = 2^^ N(m) + 2^^ N(n) 



and 

N^o(q) = ^qM + \qM . 

Thus the Illiac addition rules hold. 

Multiplication requires a shift to the left of 10 places 



We want 



Thus we require 



and 



\oM = N^Q(in) X N^Q(n) . 



2^^N(q) = 2^^N(m) x 2-^°N(n) 



N(q) = 2-^^N(ra) x N(n) . 



We may consider several simple routines to carry out the 
multiplication. The method given by (b ) is probably the best. 

(a ) The shortest method merely multiplies and shifts . 
Notice that a 75 order must be used rather than a 7^ order because 

the least significant digit of N^q^^^ ^^ ^°^ ^"^^ ^^^ ^^'^ 2"^^. 

-50 
A bias of -2 is introduced by the absence of a roundoff. The 

program is given in Table 6.1. 

(b) This method rounds of f by adding 2'^'"^ to N, ,(p), 

10 '■-^ ' •' 

giving an unbiased result. The program is given in Table 6.2. 
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50 mF N(m) to Q 

75 nF N(m) x N(n) to AQ 

00 lOF 2^^ N(m) x N(n) 
kO pF to p. 

Table 6.1 
Multiplication with Binary Point Shifted 

19 lOF 2""'""^ to A 
50 mF N(m) to Q 

-•SO 
7^ nF N(m) x N(n) + 2 ^ to AQ 

00 lOF 2"^*^ N(m) x N(n) + 2* ^ to A 

ho pF to p. 

Table 6.2 
Unbiased Multiplication with 
Binary Point Shifted 

(c) A roundoff similar to that of division is obtained 
with the program given in Table 6.3. 

50 mF N(m) to Q 

75 nF N(m) x N(n) to AQ 

00 9F 2^ N(m) x N(n) to AQ 
50 tF 2"-^ Q = 1 
00 IF This makes 2*^^ A = 1; 
N(A) = 2^^ N(m) X N(n) 
ho qF to q. 

Table 6.3 
Division- Type Roundoff in Multiplication 
with Binary Point Shifted 
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In Table 6,3 the order 50 tF is used simply because the word at 
location t has a 1 in the proper place. 

In division we also need an extra shift to restore the 
quotient to the proper range. But here the shift precedes the divide 
order and no special arrangements for roundoff are necessary. We 
want 

N^o(^) = N^o(ra)/N^Q(n). 
Hence 

2^S(q) = N(m)/N(n), 
and 

N(q) = 2"^^ N(m)/N(n). 

The following program in Table 6,h will carry out the 
required operations: 

L5 mF N(m) to A 

10 lOF 2"-^° N(m) 

66 nF 2"-^^ N(m)/N(n) 

S5 F to A 

ho qF to q. 

Table 6.k 
Division with Binary Point Shifted 

It can be noted that by using such relations as 

2"^ 
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where 10^ < 2^ we can use decimal scaling although it will be slower 
and clumsier to handle because of the factors 10^ / 2 . 

6.3 SCALING A FULL PROBLEM . There are conceptually two 
ways in which we can approach the scaling of a problem; both give 
the same program. 

(a) We can alter the problem using such 
substitutions as :^ = lOOx or p = 32p, 
so that' the modified problem has all its 
variables less than one but retaining full 
significance. 

(b) We can use scaled numbers inside the machine 
to represent the variables. That is, we can 
use x/lOO instead of x and p/52 instead of p. 

The final result must be such that the variables lie 
within machine range and retain sufficient accuracy. Constants 
greater than one can be represented by numbers less than one in 
conjunction with scaling factors. For example, multiplication by 
5.63 can be done by multiplying by 5-63/8 and shifting left 5 places 

Let us consider the following simple problem: 

Example 1. Program the Illiac to compute the quantities 

X 

e 

y = 



l + x^ 



at intervals of 0.01 from x = to x = 6. 

We note that we must scale x. Since the largest value of 
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X is 6, let us use x/8 inside the computer because this will minimize 
the loss of significant figures. 

What scaling factor is needed for y? A rough estimate shows 
that y does not exceed e / (l + 36) ~ 11, so that we shall compute 
and store y/l6. 

Library Routine S 2, the exponential routine, will give 
valid results only if x is negative. We therefore write 

.V8 .... -1 + x/8 



so that 



e X e 



X f x/8v8 
e = (e M . 



Now y/16 = 



16 X Q^ 



l/6i^ + lAJGh 



X 



and 



16 X 64 



= {[(e-^^^/^xeAf X2]2>,2p . 



We shall store the constants e/i»- and l/6i^. Each part of the computa- 
tion is within machine capacity. We proceed as follows: 



(a 
(t 
(c 
(d 
(e 



Sqimre x/8 and add 1/6^1-, 

Call the result R and store it. 

Evaluate e 

Multiply by e/U and square. 



-1 + x/8 

' with Library Routine S 2, 



Double and square. 



e^e 



(f ) Double and square again. Call the result P. 

(g) Form y/l6 by dividing P by P.. 

Accuracy in the Result . Let us now consider how much 
accuracy we obtain. If x/8 is near imity, numbers remain large 
during the calculation and we do not lose significance by sub- 
tracting nearly equal numbers or by other ill-conditioning. 

However, when x/8 is small we notice that we form e"^ / (l6 x 6h) 
quite accurately and then divide by R ~ l/64 which is about equiva- 
lent to a left shift of 6 places and losed 6 binary digits on the 
right. We can prevent this if we perform the division while the 
double-length P is in AQ. 

Accuracy in the Argument . We should also consider the 
accuracy of the argument x/8. It can be formed either by successively 
adding the increment O.OI/8 = 0.00125 or by counting and multiplying. 
The adding method is not so good because the quantity 0.00125 will 
not be stored exactly and the accumulated roundoff error obtained 
by the time x is 6 (which requires 6OO additions) may be troublesome. 

The counting method avoids this trouble. And x/8 is 
n X 0.01/8. Hence we have 

x/8 = n X 0.01/8 

= (n X 2"^^)x 0.64 X 2^^ 

-59 
If we store 0.64 and count to get, n x 2 , we can use the following 
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set of orders to get x/8 (where we have used arguments in place of 
addresses ): 



50 


-39 
n X 2 ^^ 


75 


0.61I- 


00 


30 


40 


x/8 . 



The progiram for the entire calculation of y/l6 is given 
in Table 6,^ on page 6-9 where we have again used arguments instead 
of addresses . 

Example 2; Solve the pair of simultaneous equations 

ax + by + c =0 
dx + ey + f = 

where the coefficients are in absolute values less than l/2 and 
where the answers are known to lie within machine range. Retain 
as much accuracy as is reasonably possible. 
We can distinguish two cases: 

(a) If |a| > |d| , then 



X = - 



dc/a - f 


bd/a - e 


c + by 
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ORDER 


ARGUMENT 


F5 


n X 2'^^ 


1^0 


(n + 1) X 2'^^ 


50 


(n + 1) X 2"^9 


75 


.6if 


00 


30 


ko 


x/8 


50 


x/8 


7J 


x/8 


Lk 


1/61^ 


ho 


l/6k + x^/6i^ 


L5 


x/8 


L^ 


-1 


50 


(link address ) 


26 


to S 2 


ko 


-1 + x/8 



ORDER 


ARGUMENT 


50 
TJ 


-1 + x/8 
e X 2"^ 


ko 


eV8 / ^ 


50 


e^/8 / k 


75 


e-/8 / 1. 


00 


1 


i^O 


e^/^ 1 8 


50 


e''/'* / 8 


75 


e^/^ / 8 


00 


1 


ko 


e^/^ / 52 


50 


e^/2 / 32 


75 


e^/2 / 32 


66 


l/61^ + x^/64 


S5 


(any address ) 


ko 


y/i6 



Table 6.5 



Calculation of e / (l + x ) 
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(b) If |d| > |a| , then 
b - ae/d 



X £-t-S 



In this program we shall follow the conventions of 
the Decimal Order Input (See Chapter 5)- Notice that divisions are 
always made using a full 78 digit dividend to retain as much assuracy 
as possible. The program treats the two cases (a) and (b ) separately, 
distinguishing with the 36 order at IL. Locations and 1 in the 
memory are used as temporaty storage. The program is given in Table 
6.6 on pages 6-12 and 6-I5. 

6,k ADJUSTABLE SCALING FACTORS . It is not always 
possible to arrange a program so that a single scaling factor can 
be used throughout the calculation. Then it is necessary to make 
tests at appropriate places to discover when variables are becoming 
too large or too small and to make proper adjustments in the scaling 
factors . For many problems it is advantageous to have the variable 
less than I/2 . Then two numbers can be added or multiplied without 
exceeding capacity. By using the LL n order we place* I/2 + |N(n)| 
in A. Hence A is positive if |N(n)| < 1/2 and A is negative if 
|N(n)| > 1/2. 

6.5 CONTINUOUS SCALING . FLOATING POINT ROUTINES . For 
calculations in which continual tests are required to maintain accuracy 
floating point routines (See also Section ^.7) may be used. These 
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routines represent numbers as y = a x 10 and store a and b. Thus 

they can represent accurately the numbers in some large range such 

-65 63 
as, for example, 10 < y < 10 where y has 30 significant binary 

digits. There are two such routines in the Illinois Program Library. 

The first. Library Routine A 1 is as described above. The second. 

Library Routine A ^ is a multiple precision floating point program 

-153 155 
in which numbers lies in the range 10 < y < 10 with y 

having 20 significant decimal places. 

Floating point routines are slow because numbers are 

scaled at each step of the calculation. Certain conveniences have 

been programmed in, however, and these to some extent compensate for 

the extra time required and also simplify the programming. 
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L7 22L 
L2 25L 


1 


36 15L 
50 27L 


2 


75 22L 
66 25L 


5 


SI F 
Lif 24l 


h 


i+0 F 
50 22L 


5 


75 26l 
66 25L 


6 


SI F 
Lk 25L 


7 


i+0 IF 
L5 F 


8 


66 IF 
SI F 


9 


ho 29L 
50 a9L 


10 


75 26l 
Lk 27L 


11 


66 25L 
SI F 


12 


kO 28l 
OF F 


15 


50 23L 
75 25L 


11+ 


66 22L 
■S5 F 



|a| - |d| 

|a| < |d| ; f to Q 
fa/d 

-f a/d + c to 
ae/d 

-ae/d + b to 1 

(-fa/d + c) / (-ae/d + b ) = -y 

y to 29L 

(ey + f ) 

-(ey + f ) / d = X 

X to 28L 
STOP 

|a| > |d I ; b to Q 
bd/a 

Table 6,6 
Solution of Two Simultaneous Equations 
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15 LO 26L 

i^-O F bd/a - e to 

16 50 25L d to Q 
75 2i^L 

17 66 22L dc/a 
SI F 

18 Lk 27L 

66 F (-dc/a + f)/(bd/a - e) = y 

19 S5 F 

i^O 29L y to 29L 

20 75 23L by 
Lk 2i^L 

21 66 22L (by +c)/a = -x 
22 IIL Control to IIL 

22 a 

23 ^ 
2ii- c 

25 d 

26 e 

27 f 

28 X 

29 y 

Table 6,6 (Continued) 
Solution of Teo Simultaneous Equations 
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CHAPTER 7 
MACHINE METHODS AND CODING TRICKS 

There are usually a number of special techniques which 
can be used on any particular digital computer and which will simplify 
programming. Some of these techniques are applicable on many 
different computers but usually, as is the case in those which follow, 
they result from particular orders or combinations of orders which 
are peculiar to an individual machine. This chapter is concerned 
with a number of unrelated sections having to do with operations 
which frequently arise in programming. 

7.1 THE SUMMATION OF PRODUCTS . We often need to form 
sums of products, and on the Illlac this cannot be directly done in 
the accumulator. The accuracy can often be enhanced by performing 
a summation either exactly or with only^^l3tie round-off error. This 
is comparatively easy to do using the 7^ order. All that we need 
to do is to place the least significant half of the partially summed 
products into the accumulator before performing the ^h order. In 
fact, this can usually be done by an 85 order because the quotient 

register will usually hold the last single half of a summed product. 

-39 
Then, since 7k n giyes N(n) Q + 2 A, we obtain the double-length 

product in AQ. Of course, the most significant part of the pre- 
viously summed products needs to be added using an L^^ order. 

Using similar schemes we can also arrange to add or 
subtract products with double-length accuracy in a program. As a 
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first example we shall sum 50 double-length products, assuming we 
do not exceed capacity. 

Example 1: Place the rounded sum 

Z N(100 + i) N(l50 + i) 

i=0 

in location 0. The program is given in Table 7.I 



m 


kl F 
26 IL 


Clear location for sum 
Wasted order 


m+l 


50 8L 
L5 IIL 


Put 1/2 in Q 
Set i « 


m+2 


^0 3L 






S5 F 


Bound off (See Section 7, 


w+^ 


50 ( )F 
7^ ( )F 


N(100 + i) N(l50 + i) 


m+4 


lA F 
i^O F 




m+5 


L5 9L 
Lk 3L 


Increase i by 1 


m+6 


i^O 3L 






LO lOL 


Test for i > 150 


m+7 


32 2L 
OF F 


Re-enter loop 
Stop 


m+8 


il-O F 
00 F 


Roundoff constant =» 1/2 


m+9 


00 IF 
00 IF 


Increment 


m+lO 


JO 150F 

74 200F 


End constant 


m+ll 


50 lOOF 
7^ I5OF 


Starting constant 



Table 7.I 
Program for Example 
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7.2 REVERSING THE CONTROL TRANSFER o FURTHER DISCUSSION OF 
EXAMPLE 1 . There is a second coding trick in Example 1. The end 
constant, instead of being 50 150F 7^ 200F has had -1 added to it, 
making the first order JO I5OF. The effect is to reverse the sense 
of the following 32 order (in location m + 7) so that we transfer 
control to re-enter the repetitive loop. If this had not been done 
a 22 order following the 32 order would have been required and the 
technique is equivalent to having a conditional transfer order act 
when the accumulator is negative. 

Another coding trick might have been used to save a full 
word. The left-hand order at m + 1 puts I/2 in Q so that It may be 
used to round off on the first step of the sum, this being the sole 
roundoff. Instead of storing l/2 in m + 8 we could have used the 
order pair in m + 2 as the roundoff constant. This order pair is 
1/2, plus at most 2", and would serve quite well. 

As a second example we consider a summation of two products 
with a single roundoff. 

Example 2: Given 

cos © in 10, 

sin © in 11, 

x in 12, 

y in 13> 

place the rounded quantity (x cos © + y sin ©) in 
location 20 and the rounded quantity (-x sin © + y cos ©) 
in location 21. The program is given in Table 7.2. 
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-ko 

m 50 lOF X cos © + 2 
7J 12F 

m+1 ko F Most significant half to location 
S5 F Least significant half to A 

m+2 50 IIF y sin + 2*^^ (l,s. half of x cos + 2' ) 

jk 13F 

. -Uo 

m+3 L4 F Add most significant half of x cos + 2 

^0 20F Store x cos + y sin + 2"^^ 

m+^ 50 IIF -X sin + 2 
79 12F 

m+5 ^0 F 
S5 F 

m+6 50 lOF y cos + 2"^" (l.s. half of -x cos 0+2" ) 
7^ I3F 

m+7 Lk F _^ 

40 21F Store -x sin + y cos 0+2 

Table 7.2 
Program for Example 2. 

7.3 BINARY SV7ITCHES . It is sometimes necessary to do two 
different operations alternately. This can be done by changing the 
sign of a number each time we pass it so that it will be alternately 
positive and negative. Usually it is not necessary to use a special 
number for this because some number or order pair in the rest of the 
program may be used. To accomplish the switchman LI order followed 
by a ko order is used to change the sign of the number and put it 
back in its location with sign changed. A conditional transfer order 
may then be used to decide which of two sequences will be performed. 

A variation is the requirement that an order (or order pair) 
take on two different values alternately. This can be accomplished 
by using the identities 
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b =? (b + a) - a, 
and a = ("b + a) - b. 

Thus, if the current value of an order (or order pair) is subtracted 
from the sum, the other value is obtained. 

Example 3: Arrange a program to alter the address 
of the left-hand order at (ra+2) so that it takes on 
the values and 5 alternately. A program for this 
is given in Table 7.3» 



m any order 

L5 pF Put sum of orders in A 

m+1 LO 2L Form alternate address 

i|.0 2L Store alternate address at (m+2) 

m+2 L5 (0)F 

Li|- 12F Normal program, order 

p FK 6F Sum of L5 OF Lk ,12F 
F8 2^F and L5 5F lA 12F 

Table 7-3 
Binary Switch 

In Table 7-3 only a single address is taking on alternate 
values and it is also possible to carry out the switch using the 
program of Table 7-^. In Table 7-^ the 50 5F order at location 
ra is provided for the switch. If some order which needed address 
5F could be used here, we should be very well off indeed. 
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m 50 5F 

L5 mF 

m-f-l LO (m+2)F 

he (m+2)F 

m+2 L5 (o)F 

U (12 )F 

Table 7.4 
Binary Switch 



7.4 TESTS FOR AND -1 . In order to test for a particular 
number value held inside the machine it is generally necessary to 
use two tests. However, the numbers and -1 can be tested for using 
absolute value orders and a single test. In machine l anguage is 
the onl y numb er_v£hnse -negative abfinlute value is positive ,. and -1 
is the only number whose positive absolute value is negative . 
Thus we can test for using an L3 order followed by a conditional 
transfer order, and we can test for -1 using an L7 order. Similarly, 
we can test for -2 -" and 1 - 2 '^'^ by using F3 and F7 orders, res- 
pectively. 



Example ^l-: Transfer control to location 200 if A 

is zero but transfer control to location 3OO if A 

is non-zero. Two ways to do this are given in Tables 

7.5 and 7.6. The program of Table 7.5 has only two 

words. The program of Table 7.6 has four words but 

will be faster than the other if A is usually negative. 
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Moreover, the program of Table 7.6 can be used 
to transfer control to any of three locatbns de- 
pending upon whether A is positive, negative or zero. 

m ifO F A to 

L3 F - |n(0)| to A 

m+1 36 200F To 200 if A > 0, i.e., if A = 
26 3OOF To 300 if A ;^ 

Table 7.5 
Testing for Zero 

m 36 IL 

26 5OOF To 300 if A < 

m+1 LO pF A - 2"^^ 

36 3OOF To 300 if A - 2"^^ > 0, i.e., if A > 

m+2 26 200F To 200 otherwise, i.e., if A = 

P 00 F Constant 2"^^ 
00 IF 

Table 7.6 
Testing for Sign 

7-5 USE OF ORDERS AND ADDRESSES AS CONSTAM'S . SV and 
KV orders do not use their addresses, so these addresses can often 
be used for other purposes. For instance, they may be used to store 
a starting address taken by a cycling order or an increment which 
is used to change an address. In such cases we naturally use k2 or 
h-6 orders to make certain that the function digits do not become 
altered . 
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In the following example the address of a K5 order is 
used as a counter. 

Example 5s Given the positive number a in A, 
write a closed subroutine which will furnish 
the positive integer m such that l/2 < 2 a < 1. 
The program is given in Table T-T* 



Store a at location 3 
Form link 

Plant link 
Clear counter 



Shift a 



Test to see if 2a > 1 



m 


ho 


F 




K5 


F 


m+1 


k2 


h-L 




h-5 


L 


m+2 


L5 


F 




00 


IF 


m+5 


ho 


F 




% 


5L 


m+^ 


ih 


L 




22 


( )F 


m+5 


F5 


L 




ho 


L 


m+6 


26 


2L 




00 


F 



Counter to A 
Link 

Count 



Re-enter loop 
Waste order 



Table 7-7 
Address Use in K5 Order 

7.6 RESETTING AM) STARTING OF CYCLES OF ORDERS . In 
many cases we have cycles of orders, some of which are being modified 
by the same increment. In such cases the variable addresses can all 



be modified by modifying one order and then deriving the other orders 
(or addresses) by adding the constant difference between the variable 
orders. When this is done it is economical to use the same orders 
to set these addresses when the cycle is begun. The following example 
with the program given in Table 7.8 illustrates this: 

Example 6: For i = 0, 1^ . . . , 99 place in 

location 200 + i the sum N(lO) + N(lOO + i) + N(200 + i). 

m L5 7L 

26 i^L Set i = 

m+1 L5 lOF Form N(l0) + N(lOO + i) + N(200 + i) 
L4 F and put in 200 + i 

m+2 Lk F 
i^O F 

m+3 L5 8L 

Lk 2L Increase i in m + 2 

m-f4 kO 2L Increase i in ra + 1 and test for end 
LO 9L 

m+5 h2 IL 
36 IL 

m+6 OF F Stop 

00 F Waste order 

m+7 Lk 200F Starting constant 
kO 200F 

m+8 00 IF Increment 
00 IF 

m+9 Jk 30OF End constant 

00 lOOF Constant to change i in 2L 

Table 7.8 
Resetting of Addresses 
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In Example 6 the addresses in location m + 2 have "been 
changed and then that in m + 1 has been obtained from one of them. 
Notice that the end test has been combined with the second address 
change and that the end test constant is not Lk but 7^1- so that the 
36 order at location m + 5 will cause re-entry to the repetitive loop. 

When two orders have to be varied in a single cycle, it is 
advantageous to let these form a single order pair as in Example 6 
so that they can both be modified simultaneously by the same orders. 
This arrangement is not always possible and the second best arrange- 
ment is to place the variable orders on the same side of their res- 
pective order pairs so that the orders required to modify them will 
be as simple as possible. 

For very simple operations it is sometimes advantageous to 
do three or four operations in a single cycle. This saves time 
although the advantage is bought at the expense of moi;e orders. This 
is illustrated in Example 7. 

Example 7: Add the numbers in memory locations 
10 to ik, putting the sum in location 15 . It is 
simpler and faster, both in coding and in machine 
operation, to use the program given at the top of 
page 7-11 than to write a repetitive code which counts. 
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L5 lOF 

Lk IIF 

Lil- 12F 
Li^ 13F 

L4 lifF 
i^O 15F 

7-7 USE OF THE QUOTIENT REGISTER FOR INTERCHANGES . In 
many programs we wish to replace the number in a certain storage 
location and yet use the value that is there to continue with the cal- 
culation. In such cases the old value can be placed in the quotient 
register before the new value replaces it in the memory. Thus, the 
old value is available in the quotient register for further computation, 

Example 8: Store A in location 10, but use the old 
N(10) as a dividend to form N(10) / N(11). The program 

is then 



50 lOF 

ifO lOF 

S5 F 

66 IIF 



7-S TESTING IF NUMBERS JlEE GREATER THAIf ONE-HALF . When / 
scaling numbers it is very often necessary to test when numbers are 
larger in magnitude than one-half. This can easily be ddne with the 
appropriate L or S order. For example, the order LL n will cause the 
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accumulator to be negative if the magnitude of N(n) is greater than 
or equal to l/2 . 

Example 9: If |n(10)| > l/2 replace it hy half 
its value. The program is given in Table 7.9' 



m 



LL lOF 1/2 + |n(10)| to A 
52 (m+2 )F 



m+1 L5 lOF 
10 IF 

m+2 ho lOF 

Table 7-9 
Scaling by Testing for One-half 

Y.9 COWERGENCE CRITERIA . When iterations or repetitive 
calculations are carried out ve frequently want to stop when we have 
achieved the maximum accuracy. In some cases it is difficult to 
specify in advance the tolerances which can be used as end criteria 
because we have to compromise between achieving the greatest accuracy 
and yet assuring that we terminate the processes (i.e., don't loop). 

In such cases it is worthwhile to use more complicated 
criteria which will give us maximum accuracy but which will not loop. 
One such criterion is to terminate the process if e^ > e^^^-j^ where e^ 
is some positive number which tends to zero as the process converges 
with increasing n. This criterion will terminate the process only 
when either e. is the same for two successive iterations or when the 
roundoff error has actually caused it to increase. 
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7.10 MARKING . It is often possible to use marking 
techniques instead of the more usual counting processes. The 
simplest illustration of such a technique is Library Routine N3 
where a sequence of numbers read from the tape is automatically- 
terminated by the mark N. In this code instead of counting up to 
some predetermined number we test each character as it is read from 
the tape until it is N. 

A binary digit is sometimes shifted as a marker; this is 
illustrated in the next example. 

Example 10: Using a print routine stored at locations 

beginning with 50 print the 7 numbers in locations 10-16. 

The program is given in Table 7. 10. 

o 

19 7F 2" to IF. This is the marker 
kO IF 

1 L5 IF 

00 IF Advance marker by shifting 

2 ho IF 

32 3L Test for end 

3 OF F Stop 
L5 lOF 

h 22 hh Waste order 

50 i^L 

5 26 5OF Enter print roptine 
F5 3L 

6 ^0 3L Increase address of number to be printed 
26 IL re-enter loop 

Table 7-10 
Use of a Marker 
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In the program of Table 7-10 the marker is shifted into the sign 

digit to Indicate the end of the repetitive process. 

This technique can be used in a similar way for completely 

internal programs. For example, when we are dealing with a group 

of numbers in the memory we may arrange that the storage location 

following the group contains some uniqiie number such as or -1. 

Then the code merely has to test for the presence of one of these 

numbers rather than for a predetermined count. 

7.11 REMAIKDER IN INTEGER DIVISION . In the general case 
it is difficult to compute the remainder from the residue that is 
left in the accumulator after a division. However, if we are dealing 

with positive Integers less than 2 "^ in magnitude we can do this 
quite readily. We place twice the dividend integer in the quotient 
register, clear the accumulator and divide by the divisor integer. 
The accumulator then contains twice the integer remainder and the 
quotient contains twice the integer quotient. We store an integer 

-39 
m as m X 2 

Example 11: Divide the positive integer in location 
20 by 10. Place the quotient in location 11 and the 
remainder in location 12. The program is given in 
Table 7.11. 
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m 51 20F in x 2"^^ to Q 

m X 2 -^ to AQ 

m+1 66 k-L (m x 2*^^) / (lO x 2"^^) x 2"-'- = m/lO 



m+2 40 12F Store remainder 



51 


20F 


00 


IF 


66 


hi. 


10 


IF 


40 


12F 


S5 


F 


ii-O 


IIF 


26 


pF 


00 


F 


00 


lOF 



nn-3 ho IIF Store quotient 



m+4 



Table 7-11 
Remainder in Integer Division 

7.12 BINABY CHOPPING . This is the method of repeated 
subdivision of an interval. It is easy to code although it may be 
wasteful of memory space and it is slow because it will usually take 

the full 59 steps to go from an interval of length unity down to 

-59 
one of length 2 . 

If we use binary chopping to find the zero of a function, 

we proceed as follows . We choose bounds for the zero, perhaps -1 as 

a lower bound and +1 as an upper bound. Then we bisect the interval 

and compute the function at the midpoint. Depending upon whether the 

sign of the function at the midpoint agrees with the sign of the lower 

or upper bound, we substitute the midpoint for the appropriate bound. 

-59 
After 59 steps the difference between the two bounds will be 2 

and the zero will be determined . 
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The code will be simpler if the signs of the upper and 
lower bounds are known so that a comparison with the sign of the 
midpoint is not needed at each step. This is the case in Example 
12. 

Example 12s Find the square root of a = 0. 269^1-5 

by using a binary chopping technique on the function 

2 
a - a . 
n 

Here we choose initial upper and lower bounds of 

1 - 2 ^ and 0. Instead of counting 59 steps, 

2 -37 
a test has been included so that if a - a < 2 

n 

the code stops. The program is given in Table 7 •12. 

7.13 EVALUATION OF POLYNOMIALS . Polynomials are best 
evaluated by use of a recurrence solution. Given the polynomial 

, , 5 n-i n n-1 
f(x) = L a. X = a„x + ax + . . • + ^v, 1 n ' 
i=0 ^ 

we can express it in the form 



^i+l = ^ s^ + ^1 . 



So = 0, 



Sn+1 = ^(^) • 
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L5 

10 


12L 
IF 






1 


J^o 

L5 


F 
13L 




a/2 to location 


2 


10 


IF 
F 




(a + a )/2 = a to l^L 


3 


50 


i4l 

l^l-L 






k 


7J 
LO 


li^L 
IIL 






5 


J+0 
52 


15L 




2 
a - a to I5L 

n 


6 


L5 
ho 


lii-L 
15L 




Change a 


7 


22 
L5 


8l 
li^L 






8 


40 

L7 


12L 
15L 




Change a 


9 


LO 
56 


i6l 
L 




2 
Test for a - a < 2 
n 


10 


OF 
00 


F 
F 




Stop 


11 


00 
00 


F 
26945OOOOOOOJ 


a = 0.269i^5 


12 


7L 
LL 


i|-095F 
]+095F 




i = 1 - 2-59 


15 


00 
00 


F 
F 




a = 


li^- 


00 
00 


F 
F 




a 
n 


15 
16 


00 
00 

00 
00 


F 
F 

F 


Table 


2 
a - a 
n 

2 ^' 
7.12 



Square Root by Binary Chopping 
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If, as is often the case, the coefficients a^ are the 

i 

quantities N(m+i) the recurrence relation becomes 

^i+l = ^ S^ + N(m+i ) , 



Vl = ^(^) 



Example 13 shows how a polynomial may be evaluated. In practice 
the summation of products should be done with a Jk order as in 
Example 1, but we do not wish to obscure the general idea here. 

Example 13: Given x in location k9 and coefficients 
in locations 50 through 67 evaluate the polynomial 

17 T7 , 

E x"-'"^ N(50 + i) . 
i=0 

The program is given in Table 7.13 . 

7-1^ SHIFTING A LEFT WITHOUT SHIFTING Q . The contents 
of A may be shifted to the left (doubled) by adding A to itself. 
This is slower than using a shift order but it has the advantage 
that the Q register is not altered (as it would be by a shift order) 
The order pair i^O nF Li^- nF will shift the contents of A one 
place to the left for each use. 
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p 


1+1 


F 


Sq = 




50 


^9F 


X 


p+1 


7J 


F 


X S^ + N(50 + 




lU 


50F 




P+2 


ko 


F 






F5 


IL 


Increase 1 


P+3 


i^O 


IL 






LO 


5L 


Test for end 


p+i^ 


32 


L 






OF 


F 




P+5 


LJ 


F 


End test con£ 




LU 


68F 








TalDle 7-13 


Evaluation 


of a Polynomial 
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CHAPTER 8 
CHECKING METHODS 

To obtain a correctly working program from a vritten one 
the coder must find and remove all of the coding errors . To the 
person unfamiliar with coding it might seem that a careful inspection 
of a program before it is run on the Illiac would remove all of the 
coding blunders but unfortunately this is not true. It is necessary 
to check the program on the machine to remove the remaining blunders 
from a code, and it is unusual to find all blunders in one run on 
the machine. As a general rule it will be necessary to have several 
checking runs on the machine before a code is correct. 

The number of coding errors in a program depends upon a 
number of factors, including complexity, length of untested program, 
programmer, and the care with which the code was inspected. It is 
possible to remove most coding blunders from a program by carefully 
inspecting it and it is best to obtain a person other than the coder 
to inspect the code. In many cases this cannot be done, and then the 
code should be set aside for as long as possible before being inspected. 
This to some extent prevents a grooved mind from missing the flaws in 
the program, for the coder remembers the general nature of the program 
but not the details where most errors occur. 

8.1 COMMON BLUNDERS . A great many of the coding mistakes 
made by programmers are familiar ones which are made over and over again, 
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The list given in Table 8„1 is one of common errors. It should be 
pointed out, however, that this list is not complete, and that 
programmers develop their own idiosyncrasies and should be on guard 
for their own pet blunders. 

^•2 THE DISCOVERY OF EREOTRfi. Errors are found by running 
the program and comparing the actual performance with the designed 
performance o That is, to find an error it is necessary to obtain 
Information about the way a program runs. It is almost useless to 
have a program run through a long calculation and then print out the 
result. If the result is wrong, no information is usually available 
to explain why. Therefore, when checking a program we have to print 
out more information about the intermediate results than is required 
in the actual running of the program. There are many ways in which 
this can be done and the remainder of this chapter describes some of 



them. 



Because we have to compare the intermediate stages of the 
calculation with those estimated by other means, it behooves us to 
choose parameters and variables so that the initial calculation is 
as simple as possible. Further trials will probably have to be run 
with more complicated calculations and with values designed to test 
the special cases and boundary limitations. Simple blunders should 
be detected and removed before an attempt is made to look for more 
subtle errors. Localization of mysterious errors can be effected 
by continually printing more and more information from smaller and 
smaller parts of the program until the error is found. 
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1. L5 ordere used instead of L4 orders 

2. S5 orders omitted after divisions. 

3. Orders terminated by L Instead of F and vice-versa. 

h. The renumbering of a code not complete after a 
modification has been made. 

5. Rounded multiplication used when dealing with 
integers . 

6. Control transfers to the wrong address or wrong 
order of an order pair. 

7. Accumulating storage registers not cleared before 
a cycle of orders is entered. 

8. The end condition for a cycle of orders not correct. 

9. Allowing the temporary storage of a subroutine to 
erase useful data. 

10. Using a kS order instead of a U2 order and vice-versa. 

11. . Omitting directives and starting orders on the 
program tape. 

12.. Incorrectly remembering the specifications of a 
subroutine. 

13. Forgetting to reset addresses when coming back to 
a cycle of orders . 

1^. Making corrections incorrectly. 

15. Using the same relative addresses on correction 
words although the preceding directive is different 
from that of the program. 

16. Overlooking the digits shifted from the quotient 
register to the accumulator on a left shift. 

17. Attempting to convert fractions greater than one- 
half by using the J terminating symbol. 

Table 8.1 
Typical Blunders 
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8.3 SOURCES OF INK)RMATION , 

The Punch . The punch is the most effective way by which 
data can be extracted from the Illiac. However, there are other ways 
to obtain information, particularly when a program does not run far 
enough to punch any data at all. 

The Reader o If the input tape stops before the entire 
program has been read into the memory, then an examination of the 
characters punched on the tape just ahead of the place where it stopped 
will often provide an explanation. 

The Order Register and Control Counter . If the program 
is read in correctly but comes to a sudden unexpected stop, then the 
order register R„ will exhibit the order on which the program stopped 
while the control counter will contain a number one greater than the 
storage location from which the order pair came. With this information 
the programmer can often discover the cause of failure of his program. 

If the Illiac loops, observSition of the slave tube will give 
an indication of how extensive the loop is. Then, if the machine is 
stopped and the contents of the order register and the control counter 
noted, we usually have enough data to identify the loop in our program. 

8.^ MODIFICATIONS TO PROGRAMS . Temporary modifications 
often need to be made to a program while it is being checked. These 
may be for the purpose of correcting the program, for temporarily 
arranging the punching of extra data, or for some other purpose. In 
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such cases it is undesirable to repunch the entire program tape. The 
modification can be carried out by "overwriting" the original program, 
that is by replacing some of the orders already placed in the memory 
by other orders. To do this a short extra tape is punched with suit- 
able directives and words so that, when it is read into the memory 
by the Decimal Order Input, the appropriate words of the original 
program are replaced by those on the tape. For example, 00 523K 
L5 ^F i+0 2F will cause the word at 523 to be replaced by L5 ^F kO 2F, 

In order to be able to use this technique of modifying a 
program we must: 

(a) Prevent the original program from being started 
before the modification is made. 

(b) Transfer control to the Decimal Order Input so 
that the modification tape can be read into the 
memory . 

(c) Start the program after the modification has been 
made. 

This can be done efficiently if the original program ends 
with a stop transfer of control to the Decimal Order Input followed 
by a transfer of control to the program. An example might be 2^ 999N 
26 93N if it were necessary to transfer control to location 93 to 
start the program. 

If we wish to modify the program before starting it, all 
that is necessary is to place the modification tape in the reader 
before moving the switch to START. This will cause the modification 
tape to be read instead of the 26 93N. The modification tape will 
naturally have to end with 26 93N in order to start the program. 
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If the original program ends with a stop transfer of control 
to itself instead of to the Decimal Order Input a slightly more compli- 
cated modification tape is necessary. An example might be 2i<- 93N. 
Then the sexadecimal order pair 26 3F7 00 000 is punched at the 
head of the modification tape. When this tape is placed in the reader 
the machine is completely restarted by setting the order register to 
82 ^OF i^O F and the control counter to zero without, however, clearing 
the memory. This results in the first order pair on the tape (namely 
26 3F7 00 000) being transferred to location zero and being obeyed, 
so that control is transferred to 3F7 sexadecimal or 999 decimal (the 
Decimal Order Input), and the rest of the tape is read in the usual 
way. This is a bootstrap start. See Section 5. 11. 

Corrections. As each coding blunder is found, a modification 
tape should be repunched to include all of the corrections. It is not 
worthwhile repunchlng the entire program tape until all or a large 
number of coding errors have been found. 

8.5 BLOCKTHG ORDERS . This is the name given to control 
transfer orders which are used to replace normal program orders, so 
that some printing or checking can be done at the point of replace- 
ment. Before control is restored to the program the replaced order 
of the program is executed and the contents of the arithmetic registers 
are restored. Thus the original program is unaffected, but the 
extra orders that are obeyed can be utilized to do printing of desired 
data . 
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An example may illustrate this: Let us suppose that 

(a) we have a program In which we wish to 
print N(19) after the left-hand order of 
the order pair ^i-O 9F L5 29F at storage 
location 100^ has been obeyed, 

(b) we wish to preserve the quotient register 
but not the temporary storage of the print 
subroutine, 

(c) locations 800-80^ are unused by the program, 

(d) the program print routine starts at 200, 

(e) the original program ends with 2k 999N 26 nN. 
We then prepare the following modification tape: 



MODIFICATION TAPE 



00 


lOOK 


ko 

26 


9F 
BOOF 


00 


800K 


S5 

ho 


F 

kL 


L5 

50 


19F 
IL 


26 
L5 


200F 
29F 


50 
26 


kL 
lOlF 



26 nW 



Directive 

Plants blocking order in progr 

Directive 

Stores Q at 8o4f 
Enters print subroutine 

Does omitted order 

Restores Q 

Control back to program 

Starts program 



8.6 TYPES OF CHECKING ROUTmES . There are several types 
of checking routines. One prints out the contents of certain 
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memory locations after a program has stopped. We call this a post 
mortem routine . 

A second type takes a given program and allows it to be 
obeyed order-by-order while printing out information about the course 
of the program. It is often called a sequence checking code. 

A third type arranges for information to be printed out at 
specified points of a program. It is usually known as a check point 
or blocking order routine. 

8.7 POST MORTEM ROUTINES . Library Routines C3, C^, C^ . 
These are printing routines and have been arranged so that they can be 
used with little or no preparation on the part of the programmer. They 
are read into the memory with bootstrap input routines and are located 
in storage locations at the end of the memory, the longest occupying 
locations from 986 to 1023 and using locations 0, 1, and 2, as temporary 
storage. The end of each of the post mortem tapes contains 100 two- 
decimal-digit numbers. These numbers are used to specify the locations 
from which printing will occur. Suppose, for example, that we wish to 
know the order pairs in memory locations 5^0 to 5^9. Then we read in 
Routine 05. When it stops we find the number 5^ on the end of the tape 
and place it in the reader. When the START switch is moved, the order 
pairs in locations 5^0 to 5^9 will be printed and then the Illiac will 
stop. If we start again we will get the order pairs in locations 550 
to 559, Codes 03 and Ch perform similar functions for decimal fractions 
and decimal integers. 
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The programmer should keep in mind the storage locations 
used by the post mortem routines so that he will be able to take 
full advantage of these important checking routines. 

8.8 POST MORTEM VERSION OF THE DECIMAL ORDER INPUT . Library 
Routine CI . This is a very important checking routine and is usually 
the first one used after a program failure. It compares the contents 
of the memory with the contents of the input tape. Only discrepancies 
are printed out, enabling programmers to discover which orders of a 
program have become altered while the program was in the memory. This 
is an aid in making sure that orders which should have been modified 
have been modified correctly and that no order has become modified 

accidentally. 

This routine is used in the following manner: The Post 
Mortem Version of the Decimal Order Input is read into the memory 
in the usual way with a bootstrap input routine. It occupies storage 
locations 962 to 1023. The program tape is then placed in the reader 
after its copy of the Decimal Order Input, so that the Decimal Order 
Input is not read into the Illiac. Then when the machine is started 
the program tape will be read and the words created from it compared 
with the corresponding ones in the memory. When a discrepancy is 
found it will be printed on a line giving, first, the location at 
which the descrepancy was found, second, the word read from the tape 
printed as an order pair and, third, the word found in the memory, again 
printed as an order pair. Thus a typical line of printing might be 
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3^5 L5 000 i^O 35^ L5 5^6 ko 35!^ 

indicating that the left-hand address of the order pair in location 3it-5 
had become modified, taking the value 5^6. 

It should be pointed out that closed subroutines which have 
been used will usually have their links printed. This gives an indi- 
cation of the part of the program from which they were last called in. 
Interludes cause a large amount of printing because the contents of 
the interlude locations are changed twice in the course of input and 
are printed out both times . 

The Post Mortem Version of the Decimal Order Input causes 
the memory to be changed to the original state as it is being executed, 
so that if the program is started the original program will be performed 
again. 

If desired, only selected parts of the tape need be compared 
with the contents of thememoryo However, when doing this it is 
necessary that all the preset parameters pertaining to that part 
of the program be input and that the selected part of the program 
begin with a directive » 

All post mortem codes occupy locations at the end of the 
memory and all of them use storage locations 0^ 1, 2 as temporary 
storage, so that it is desirable when coding not to use these storage 
locations for constants or numbers which may need to be printed out. 
It is worthwhile to note that if addresses which are to be changed 
are initially read in with their final values then no printing will 
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take place on the post mortem if they have been modified correctly. 

8.9 THE ADDRESS SEARCH ROUTIME . Library Routine C2 . 
Programs sometimes fail because of a transfer of control to an order 
which causes the machine to stop. The usual order causing the stop 
is a zero left shift order because the memory is normally cleared to 
zeros before a new program is read in. 

Under such circumstances the offending transfer of control 
order may be searched for with the aid of Library Routine C2. The 
search routine is read into the memory, occupying storage locations 
normally occupied by the Decimal Order Input. Nex.t the address to 
be searched for is I'ead into the machine as a three character sexa- 
decimal address. The routine then searches the memory (exclusive 
of itself) for order pairs containing this address. When found they 
are printed out (in sexadecimal form) together with their location 
(in sexadecimal form). 

The routine naturally has other similar uses. For example, 
if it is known that some number becomes modified but it is not known 
why, then the store order which does the damage can be sought for in 
the above manner. 

8.10 SEQUENCE CHECKING CODES . Library Routines D2 and D3 . 
These routines control a program order by order and print out suitable 
information about the execution of each order. This enables the action 
of a program to be traced order by order. They use a blocking order 
technique which enables a selected part of the program to be checked 
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In this manner o These codes are very slow on account of the printing 
involved and should not be used blindly . 

Routine D3 prints out the function digits of the orders which 
are actually obeyed, starting a new line of printing whenever a control 
transfer order is obeyed » This enables the sequence of orders obeyed 
by a program to be traced. 

Routine DQt prints in full each order that is obeyed, having 
one order pair per line of printing. After each store order or address 
order, the number transferred to the memory is also printed. This 
enables every step of a program to be completely checked. However, 
the amount of printing is such that this routine is very slow. It 
should be used only in the final stages of tracking down an elusive 
fault . 

8.11 CONTROL TRAKSFER CHECK . Library Routine D4 . This 
routine takes charge of a program and allows it to be obeyed order 
by order. Each transfer of control that is obeyed is placed in a 
list kept in a specified piace. The list is cyclic, that is to say, 
the later entires overwrite th« earlier ones in a cyclic fashion. At 
the end of a program, the list can be printed, so that it can be dis- 
covered how the program reached its final end. There is no printing 
during the execution of the program, so that this routine allows the 
program to be obeyed at speeds much greater than those of the routines 
which print. 
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6.X2 THE CHECK POIHT ROUTIME . Library Routine Dl. This 
routine Is designed to print out Intermediate Information about some 
other program. It uses the hlocicing order principle, and the programmer 
prepares a specification tape to descrlhe the kinds of Information he 
wishes to Obtain. It Is possible to go through Iterative loops and 
print results on various passages through the loops. Data can be ob- 
tained as an order pair, a right-hand address, a left-hand address, 
a 10 character sexadecimal word, a signed integer, a signed 12 decimal 
place fraction or a signed 5 decimal place fraction. 

Library Routine Dl Is a very powerful checking routine be- 
cause of the great latitude glven-the programmer in choosing where 
and how he will obtain information and because it utilizes the pro- 
grammer 's own knowledge of his code. 

8.13 COHCLUSION . The fact that it will be necessary to 
check a program should be kept in mind when writing it. Its storage 
locations, Including those of the temporary storage, should not inter- 
fere with routines likely to be used in diagnosing its faults. The 
program tape should end with a stop control transfer to the decimal 
order input as explained in Section 8.1*. 

Occasionally a blunder turns up in the blind spot of the 
programmer and it appears to be impossible to find it. It is not much 
comfort to him to point out that the blunder can be found by continually 
narrowing down the section of the code in which it is known to be. This 
is tedious but eventually all blunders yield to this technique. A final 
„ord to this chapter might be that all programmers make mistakes but good 
programmers find theirs first. 
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CHAPTER 9 



TAPE PREPARATION 



9.1 THE ILLIAC INPUT , The input unit of the Illiac 

is a photoelectric tape reader that transfers binary digits from a 
punched paper tape to the A register. The tape prepai-ation equip- 
ment is used to translate instructions and data from the programmer's 
manuscript into a hinary-coded punched tape acceptable to the Illiac, 

9.2 THE ILLIAC OUTPUT . Output from the Illiac is 
usually in the form of punched tape. This tape may be printed on 
any of the Teletype page printers that are a part of the tape jpre- 
paration equipment. The page printer performs a conversion from 
the binary- coded representation on the tape to the sexadecimal or 
decimal characters on the printed page. 

9.3 THE PERFORATED TAPE . SEXADECIMAL TAPE CODE CHARACTERS , 
The tape preparation .equipment presently in use consists of Teletype 
equipment which has been modified to correspond to the Illiac tape 
code. This tape code is shown in Figure 9.1- The paper tape can be 
punched in any one of 6 positions across its width. One of these 
positions is always punched. This is the feel hole or sprocket hole 
position. Peed holes are smaller than the other holes. 

Four of the remaining positions are used to represent the 
2^ sexadecimal characters. There are 16 keys on the Teletype tape 
punch labeled through 9 and K, S, N, J, P, L. When one of these 
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0123^56789KSNJFL 




Figure 9.1 
Sexadecimal Tape Codes 



Delay- 



Space 




Tabulate 



Carriage Return 
and Line Feed 



Figure 9.2 
Format Tape Codes 
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keys is pressed a corresponding pattern of holes is punched across 
the width of the tape. Since a hole or absence of a hole is a binary 
affair, we speak of the character punched in the k positions across 
the tape as a binary tape code for the corresponding sexadecimal digit 
If we speak of only the 10 sexadecimal characters 0-9 and their 
corresponding binary codes, we call the characters punched in a tape 
binary-coded decimals . 

The Illiac tape code for sexadecimal characters is very 

easy to memorize since the hole positions across the tapfe simply 

3 
coi'respond to powers of 2 from 2 to 2 . 

Page Printer Format Characters . Punched tapes are con- 
verted to a printed page by means of a page printer. Four addi- 
tional tape codes shown in figure 9.2 are used for format control, 
causing the printer to space, return the carriage and advance to 
the next line, tabulate, and delay, i.e., perform one printing 
cycle without printing anything. The functions are similar to 
those of a typewriter, but the tabulate feature occurs on only- 
one printer. Successive tabulate characters set the carriage to 
10, 25, 40 and 55 spaces from the left margin. The carriage 
return-line feed and tabulate characters must always be followed 
by; a delay character to allow sufficient time for the carriage to 
move. 
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Notice that these codes all have a hole In a position that 
is never punched for a sexadecimal. This is called the ^th hole 
position . When a punched tape is placed on the Illiac tape reader 
the circuits cause any character with a 5th hole punched to be skipped 
over, i.e., not read into the computer, when reading in the normal 
way. (However, a special input order will read such a character; 
see the order code.) This means that in the preparation of an 
instruction tape, these format codes can be interspersed with sexa- 
decimal order digits in any desired way. Then when an instruction 
tape is printed these tape codes will control the printer. The usual 
method in preparing instruction tapes is to follow each order pair with 
a carriage-return and line feed code. This produces a single column 

print of words . 

Tape Codes for Letters . The printer will also print all 

the letters of the alphabet for identifying headings and the like. 

Tape codes are provided for these and are shown in the complete tape 

code list. Figure 9.3- To accomodate all of these symbols, there 

are also numbers shift and letters shift tape codes, analogous to the 

shift- lock on an ordinary typewriter. 

A list of instructions to be used with the Illiac and the 

Decimal Order Input (See Chapter 5) to produce any of the tape codes 

shown in Figure 9-3 is given in Table 9.1. 
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P 







Q 
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R 




■h 


S 




- 


T 




5 


U 




7 


V 




y 


W 




2 


X 




1 


Y 




6 


Z 




X 



MACHINE ORDER 
(DECIMAL ADDRhiSS) 

92 963F 
92 13 IF 
92 515F 
92 259F 
92 707F 



92 


387F 


92 


195F 


92 


83 5F 


92 


67F 


92 


194F 


92 


898F 


92 


579F 


92 771F 


92 


514F 


92 


834F 


92 


642F 


92 962F 


92 6/^3F 


92 


770F 


92 


578F 


92 


2F 


92 


66F 


92 


258F 


92 706F 


92 


322F 


92 


450F 


92 


323F 


92 


130F 


92 


/v$lF 


92 386F 


92 


899F 



Tahle 9-1 
Coded Tape Operation of the Teletype Printer 

■K)nly one printer has a tabulation mechanism. 
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9.1^ CLASSIFICATION OF OPERATIONS . The operations in 
preparing tapes may be listed as follows : 

1. Punching a tape from a manuscript by 
means of a keyboard tape punch. 

2. Punching a tape from one or more previously 
punched tapes. This operation is called 
reperforating and is used for duplicating 
tapes, joining short tapes, and making 
corrections . 

3. Printing a tape on a page printer. 

k. Comparing two tapes to see if they are 
identical. 

5. Comparing the manuscript- with a printed 
copy (proof-reading). 

It is necessary to punch the main part of the code from 
a manuscript by means of a keyboard. It is not necessary to punch 
routines that are in the routine library. 

When all the parts of the code that have been punched from 
manuscript are correct, obtain any necessary routine tapes from the 
routine library. Using the reperforator, these tapes can be assembled 
into one long problem tape » 

Every time a tape is reperforated it is wise to check the 
copy by comparing it with the original tape. 

9.5 DESCRIPTION OF EQUIPMENT o Several Teletype units 
are provided for making instruction tapes . These include a keyboard 
perforator unit, a page printer unit, a reperforator unit, a combination 
unit, and a tape comparer unit. 
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Keyboard Perforator Unit . This unit is used to punch tape 
from a manuscript; see Figure 9.^. 'Turn on the power switch under the 
table by pulling it forward. The power supply takes about a minute to 

warm up. 

The sixteen sexadecimal keys that will be used most are 
in the center of the keyboard and arranged to be operated with one 
hand. It is recommended that a touch system be learned both for maxi- 
mum accuracy and maximum speed. 

If a wrong key is pressed, operate the backspace lever and 
then the space bar. This will punch the space symbol (all five holes) 
over the incorrect symbol and so be skipped on the Illiac reader if 
the usual 80 instructions are used. 

Page Printer Unit . This unit is used to make a printed 
page corresponding to a punched tape; see Figure 9-5. Turn on the 
power switch. The power supply takes about a minute to warm up. 

In order to operate the printer mechanism the page printer 
motor switch should be turned on. 

The function of the transmitter-distributor is to read paper 
tape and translate the binary tape codes into electrical signals. 
These signals are sent to the page printer. The binary tape code 
is sensed by 5 pins which are periodically pressed against the tape. 
A hinged lid holds the tape down against the pins. The tape is moved 
by a sprocket wheel which engages the feed holes punched in the tape. 
Place a tape in the reader and close the lid. Then turn the start- 
stop switch to start. The tape will now be transmitted to the page 

printer. 
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Figure 9.5 
PAGE PRINTER UNIT 
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There is an end of tape stop device that stops the transmitter- 
distributor when the end of the tape passes under the lid. A tight tape 
stop bar will also stop the unit and keep the tape from being torn should 
the tape become snarled or tight and raise the bar. 

Reperforator Unit . This unit is used to duplicate tapes; 

see Figure 9.6. 

Turn on the power switch by moving it towards the back of the 
table. After about a minute, the power supply will warm up and the 
reperforator motor switch should be turned on. 

The transmitter-distributor has been described in the pre- 
ceding paragraph. Place a tape in the reader and turn the start-stop 
switch to start. A copy of this tape will be made by the reperforator. 

When a few corrections need to be inserted in a tape, punch 
a second tape containing only the corrections. Then use the reperforator 
unit to duplicate the correct portion; stop; insert the second tape in 
the reader and reperf orate the correction; then replace the first tape 
and continue with the correct portion. 

Combination Teletype Unit . This unit, shown in Figure 9-1, 
combines the functions of those units already described. 

Underneath the top of the table are two switches. Switch 
10 turns on the dc supply which must be on before any of the other 
units will operate. On some of the tables it takes about a minute 
for the dc supply to warm up. 

Switch 9 turns on the motor of the transmitter-distributor. 
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Figure 9.6 
REPERFORATOR UNIT 
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Figure 9.7 
COMBINATION TELETYPE UNIT 



The transmitter-distributor has already been described. It sends 
signals to both the page printer and the reperforator and will operate 
either or both of these units if the corresponding motor switch is on. 

Turn on the reperforator motor by means of switch I3. The 
reperforator will now duplicate any tape that is run through the 
transmitter-distributor . 

Pressing lever 11 will cause the reperforator to punch 
space tape codes. Pressing button 8 will cause the reperforator to 
punch zeros. However, one unpredictable character will be punched 
at the end of the zeros. Zeros are useful at the beginning of tapes 
to provide space for writing identifying information. They should be 
followed by spaces to make it easy to set the tape in the Illiac 
reader. 

To operate the page printer turn on motor switch 7B. The 
page printer will now print any tape that is run through the trans- 
mitter-distributor. The carriage may be returned manually by pressing 
lever 2. Crank 1 will turn the platen. 

It is possible to operate the page printer and the reperforator 
at the same time. However, it is possible for the reperforator mechanism 
to make an error that would not appear on the printed copy. If the 
printed copy that is to be proof-read against the manuscript is pre- 
pared from the final tape, it is quite certain to show any errors in 
the tape as well as any possible errors in the printing mechanism. 
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Set switch 7A to "tape". Tape can now be perforated from the 

keyboard . 

It is also possible to type on the keyboard, and perforate 

type in the reperforator at the same time. To do this turn TA to 
keyboard, 7B on, and 15 on. T^ere is no way to back-space tape that 
comes from the reperforator. It is also more difficult to read the 
last character that has been punched on the tape than when the keyboard 
perforator is used for punching tape. The transmitter-distributor 
may be used to reperf orate correct portions of a tape. The tape should 
be stopped before the error by means of switch 5- Type the corrections 
from the keyboard. TTaen move the tape to the next correct section 
and continue reperf orating . 

It is also possible to print with the page printer at the 
same time a tape is being punched from the keyboard. To do this turn 
on the printer motor by means of switch TB and set switch TA to 
"keyboard and tape". A disadvantage of this method is that the keys 
cannot now be operated as fast as when punching tape alone. There is 
no way to change errors on the printed copy that are corrected on the 
tape with the backspace lever. Since the printed copy is made by a 
separate mechanism at the same time that the tape is punched, there 
is also a small possibility that it may not correspond to the holes 
in the tape. This copy should not be used for proof-reading. 

Switch TB can be turned on and off without affecting anything 
being reperforated if switch TA is placed in the "tape" position. 
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FIGURE 9.8 



PHOTOELECTRIC PAPER TAPE COMPARER 

CONTROL PANEL 



before TB is operated, and kept there while 7B is off. 

Comparing Tapes . A tape comparer consists of two tape readers 
connected to checking circuitry. Two kinds of tape comparers are avail- 
able in the Laboratory, one kind using mechanical readers and the other 
using photoelectric readers. When tapes are compared in the mechanical 
readers, it is necessary only to place the two tapes in the readers and 
to turn the power on. The tapes will advance in synchronism as long as 
the binary codes in the two tapes agree and will halt if a disagreement 
is detected. On of the machanical units will halt one row of holes 
beyond the disagreement and the other will halt on the holes which dis- 
agree. When the tapes stop, the power should be turned off and the tapes 
inspected visually. They may then be marked if in error, replaced in 
the readers and the comparison continued. One of the two mechanical 
units has push-buttons for advancing either tape one step at a time 
while they are stopped on a disagreement . 

The tape comparer using photoelectric readers is more versatile 
than those with the mechanical readers . By means of a switch it is 
possible to choose three modes of operation in which the reader will 
"skip no characters", "skip all 5th hole characters" or "skip all spaces." 
A set of neons is associated with each tape reader and when the tapes 
stop on an error the neons indicate the character read from each tape. 
As in one of the mechanical readers, it is possible to advance either 
tape one step at a time by means of switches . A full description of 
the operation of the photoelectric tape comparer is provided in Internal 
Report 62, and a view of the control panel is shown in Figure 9-8. 
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Short tapes may be checked more rapidly visually by super- 
imposing them and holding them up to a strong light. 

9.6 GENERAL REMARKS. The misreading of even a single 
hole in a long tape can cause enormous changes in the behavior of 
Illiac. It is important to keep tapes clean. Dirt from tapes 
accumulates in the transmitter-distributor tape reading pins and 
ca^lses errors. Worn or torn tapes will also cause errors. 

Do not mark tapes with a waxed pencil: use a lead pencil 
or ink. 

Library tapes and other tapes which need to be handled 
frequently are on a grey parchment stock which is heavier than normal 
tape stock. One of the reperforators is mechanically adjusted to 
punch this heavy stock. 

It is generally more convenient to make corrections to be 
made on the tape itself, as well as on the printed copy. 

If corrections are needed at the end of a long tape, it is 
sometimes more convenient to reperforate the tape tail-end first, 
so that the corrections can be made at the beginning, and then the 
reperforator left running unattended . 

A hand punch is available and can often be used to change 
one or two characters and so avoid reperf orating a long tape. Holes 
punched with it should be inspected carefully to see that they line 
up with the other holes in the tape. If in doubt, it is better to 
have the extra hole slightly oversized so that it is certain to be 
read by the Illiac reader. 
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9.7 HIGH SPEED PRINTER . This printer performs the 
same function as the slower printers except that it operates at 
10 characters per second. The turn-on switch is on the lower 
front side. Carriage return and line feed buttons are provided 
on the front. After printing is complete the paper should be 
advanced before shutting the printer off. Otherwise unwanted 
characters will be printed on the page as the unit stops. 

9.8 HIGH SPEED REPERFORATOR . The high speed reper- 
forator uses two photoelectric readers to provide a flexible 
means of preparing tapes. Four switches and a push button con- 
trol the operation as follows. 

Switch S4, Stop-Run . This switch has only two positions 
which determine whether the unit is in a standby position or 
operating in accordance with settings of the other switches. 
It should be at STOP to change any switch except S2. 

Switch SI, Reader Selection . This switch has three posi- 
tions, allowing a tape to be read from the front reader or the 
back reader or allowing carriage return-line feed characters to 
be punched on the new tape while neither reader operates. 

Switch S3, Rode Control . This switch, in conjunction with 
S2, determinesthe mode of reading a tape and punching another. 
It may be set as follows: 

a. "Single Character Punching" to punch one character and stop. 

b. "Stop Punching on CR & LF" to punch up to and including the 
next GR & LF character and then stop. 

c. "Single Character Bypassing" to skip one character and stop. 

d. "Stop Bypassing on CR & LF" to skip all characters up to and 
including the next CR & LF and then stop. 
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Switch S2, Start, Obey , Ignore , This switch obeys the 
instruction of S3 when set to OBKY. If moved to START after a 
stop it will cause the instruction of S3 to be executed once. 
If set to IGNORE it will cause all of the stops specified by 
S3 to be ignored. 

Set Button . This is a reset control which should be 
pushed when the unit is first started and v-fhenever- a tape is 
manually shifted in either reader. 

Neon Indicators . A row of neons on the panel indicates 
the last character read and either punched or bypassed. 
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CHAPTER 10 



CALCUIATION OF RUNNING TIME 



Each programmer should estimate as well as he can the 
length of time his program will run. This information is needed 
for efficient scheduling of machine time and to enable the computer 
operator to decide whether a program may have failed to stop when 
it should have. 

The estimation of running time consists of summing up the 
times for executing all of the orders of the program, taking into 
account that many orders are executed more than once. Thus, the 
time is calculated for one passage through an iterated loop and then 
multiplied by the number of trips through the loop. The number of 
passages through a loop is not always known (as in the square root 
code where the number of passages depends upon the quantity whose 
square root is being calculated) and the programmer may have to 
make rough estimates for such cases. He can make use of the code 
checking periods to help him here . 

10.1 ORDER TIMES. The progreunmer must know how much 
time the Illiac requires for each order. The order times are given 
in Table 10.1, the time required for getting the order from the 
memory being included . 

The values given in Table 10.1 are in some cases not 
exact. The 70« microsecond value for multiplication is a maximum, 
the time depending upon the multiplier. It could be as low as 625 
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microseconds . 



ORDER TYPE 

On, In 

2, 5, h, 5, J 

3 not executed 

6 

7 

80, 81 

82, 92 

K, S, F, L 



TIME 

16 n microseconds 
55 microseconds 
l8 microseconds 

800 microseconds 
700 microseconds 

h- milliseconds/character 

17 milliseconds/character on punch 
1 millisecond/character on display 

90 microseconds 



Table 10.1 
Order Times 

The time required to execute an 80 or 8l order so as to 
read one charactei^ from the tape is k milliseconds. Since most 
programs use conversion routines for input it will be found that 
the time is closer to h l/2 milliseconds for decimal input. 

1^-2 EXAMPLE OF RUMING TIME CALCULATION. Let us con- 
sider a simple example. The program given in Table 10.2 will trans- 
fer the contents of memory locations 100 to 199 into memory locations 
200 to 299- How long will it take? 

The program in Table 10.2 consists of 7 words. There is 
a loop consisting of the 9 orders beginning with the left-hand order 
at location and ending with the left-hand order at location k. 
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This loop is executed 100 times. There are no other orders executed 
except the stop order. Thus we have the following orders: 

00 50K 

L5 (lOO)F 
kO (200)F 

1 L5 L 
LO 5L 

2 32 hL 
L5 L 

3 Lk 6l 
UO L 

U 26 L 
OF F 

5 L5 199F 
UO 299F 

6 00 IF 
00 IF 
2^^ 50N 

Table 10.2 
Repetitive Program 

Five L orders at 75 microseconds. 
Two h orders at 55 microseconds. 
One 2 order at 55 microseconds. 
One 3 order at 18 microseconds. 

The time T is then given as 
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T - (5 X 75 + 2 X 55 f 55 + l8) 100 microseconds, 
= 53>800 microseconds, 
= .05^ seconds. 

10.3 A SIMPLE RUMING TIME FORMULA . A rule which is 
accurate enough for most calculations is the following one: 

Let Nq = number of orders obeyed, 

\ = number of multiplication orders obeyed, 
N^ ^ number of division orders obeyed, 
N = number of shifts . 

Then the running time T, exclusive of input and output orders, is 

rr - ^ ^ s m d 

-^ - 1^ "^ "^O" 2 milliseconds. 



In the example given above N^ is 90O while N = N ^ W = 

^ s m d 

We thus have T = 9OO/16 milliseconds or about ,056 seconds. 
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CHAPTER 11 
PREPARATION OF A COMPLETE PROGRAM 

The problems solved on the Illiac vary in complexity from 
programs requiring only the punching of parameters and data for use 
with a library routine to programs involving hundreds of words re- 
quiring the entire memory capacity of the machine. There is, of 
course, no "typical" program; however, the problem described in the 
following sections illustrates features common to many problems. 

The problem described is one that arose in the electron 
tube research laboratory of the University of Illinois.-^ The problem 
is first described qualitatively, and a general m3,thematical formu- 
lation is given. Approximations are then made to simplify the 
equations . The resultant equations are reformulated so that a 
library routine can be used for the solution, and reformulated again 
in terms of scaled variables so that machine range will not be ex- 
ceeded. The general organization of the program is then described, 
and the details of the coding are given. Those interested primarily 
in the programming may omit the introductory description. 

11.1 GENERAL STATEMENT OF THE PROBLEM . The purpose 
of the problem is to solve for the trajectories of electrons in an 



■^ The staff of the computer laboratory is indebted to Mr . Irving 
Kaufman of the electron tube research laboratory for permission 
to describe his problem and for his assistance in the preparation 
of this chapter . 
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electron cyclotron. An electron cyclotron is a device that accelerates 
electrons to energies of several million electron volts. The accelera- 
ting mechanism is a radio-frequency electric field in a microwave 
cavity. According to the well-known Lorentz force equation, the force 
on a charged particle is given by: 



Here 



f = qE + qv X B , (ll.l) 

f = vector force (newtons)^ 

q = charge (coulombs), 

E = vector electric field (volts per meter), 

V = vector velocity (meters per second), 

B = vector magnetic flux density (webers per -square meter). 



It is seen that the component of force produced by the 
electric field E is parallel to E ; that produced by the magnetic 
field B is perpendicular to both v and B. The former causes accelera- 
tion or increase in energy; the latter merely changes the direction 
of motion. This combination gives rise to the mechanism of the 
electron cyclotron shown in figure 11.1. 

Here a microwave cavity, externally supplied with electri- 
cal energy, is immersed in a uniform steady magnetic field (shown 
by the dots) whose direction is perpendicular to the paper. The 
cavity has been constructed so that near the axis it has a time- 
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varying electric field E that is almost entirely parallel to the 

axis X-X' and whose magnitude is given by E = E^ sineiOt. The inagni' 

tude of E is on the order of 5 x 10' volts per meter, so that it 
m 

is possible for an electron passing through the 1 cm. gap of the 
microwave cavity to acquire an energy of 5 x 10 electron volts . 
(This corresponds roughly to the relativistic rest energy of the 
electron). 




Microwave 
Cavity 



Figure 11.1 
Electron Cyclotron 
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The^action of this cyclotron is as follows: 

1. The very high value of E forces some electrons out of 
the cavity wall during some parts of the cycle. (Let us consider 
here the trajectory of electrons emitted from surface R-R). 

2. The E-field accelerates the electron toward sux-face 
S-S. 

3. A fraction of the electrons emitted pass .through th* 
orifice in S-S. The magnetic field as well as space charge pro- 
duces a slight bending of the trajectory. 

h. The magnetic field external to the cavity causes 
the electron to move in circular orbits. The orbit radius is pro- 
portional to the product of electron speed and total (relativistic ) 
energy. Consequently, a fraction of the electrons ejected are 
returned to the orifice in cavity surface R-R5 the rest collide with 
the walls . 

5. Those particles that re-enter the cavity when the field 

is in its accelerating phase (recall, E - E sincOt) acquire 

m ^ 

additional kinetic energy during transity through the gap. They conse. 
quently execute larger orbits externally. 

6. For those electrons that maintain the favorable phase 
conditions of (k) and (5), above, an indefinite number of orbits 

is possible. For a 1 cm. cavity gap, each orbit corresponds to a 
step in energy of approximately I/2 million electron volts . 
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11-2 jg_A.:JECTORY EQUATIONS . The combination of the 
reiativistic equation of motion for a particle and the Lorentz forc^ 
equation yields, for electron motion in the cavity, the equation 



d_ 
dt 



m (iv + jv ) 

Ox y 



1 - 



2 2 

V + V 



1/2 



= e [ i (E^ sincot - Bv ) + j (v B) ] (ll 



Here: 

m = rest mass of electron, 

v^ = dx/dt = component of velocity in X-direction, 

V = dy/dt = component of velocity in Y-direction, 

c - velocity of light, 

i == unit vector in X-direction, 

j = unit vector in Y-direction, 

t = time, 

e = magnitude of electronic charge, 

E^ = peak electric field intensity, 

<~o = angular frequency of cavity field, 

B = magnetic flux density. 

For motion exterior to the cavity, the term E sin^t 

m 

should be deleted from equation 11.2. The two resulting vector 
equations are then to be solved for electron trajectories. 

To reduce the complexity for a first order solution, 
three approximations are made: 
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1. The effect of the magnetic field inside the cavity 
is neglected. 

2. The effect of v in the gap is neglected. 

5. Particles exterior to the cavity are taken to move 
in circular arcs. The arc length is taken as the difference be- 
tween a full circle and the cavity gap length, (in other words, 
the chord length of the gap is set equal to its arc length.) As 
before, the radius of the circle is proportional to the product 
of velocity and total energy. These approximations leave only one 
differential equation to be solved; for now, after some manipulations, 
in the gap 

d [ m^ V (1 - p2) -1/2 ] / dt = e E sina;t , 



m 
while outside of the 



(110) 



0. _.. = K^ [ 1 - K^ (1 - p2) ^/^ /p ] (1 - p^) -1/2 



i,n+l 



^-©e,n - (2«) (n-^1) . (H.i,) 



The symbols used in 11. 3 and 11. i^ and not previously 
mentioned are: 

V = electron speed (meters per second), 

^ = v/c, 

n = number of the orbit in which moved; i.e., for 

first cavity crossing and exterior orbit n = 1, etc 
®e,n " "^^^ phase angle (radians) in the microwave RF cycle 
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^. 



i_,n+l 






at which the electron leaves the gap in orbit n. 
(See Figure 11.2 ). 

the phase angle (radians ) in the RF cycle at which 
the electron re-enters the gap after completing 
orbit n. (See Figure 11.2). 

6.51599 

0.092910 



LJ 




i,n+l 



Figure 11.2 
Phase Angles of an Electron 
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^^- • ^ ^O CESSING OF EQUAT IONS FOR SOLUTION BY THE ILLI AC . 

&:£:^.?iL.IrM£2_:tory_m^^_t^ The equation of motion of 
the elctron in the gap is given by equation 11. 3, which is reduced 
to a set of first order ordinary differential equations so that the 
solution can he obtained using library routine F 1. 

Using as variables Q =ujt, u = v/u), we can reduce equation 
11.3 to the pair of first order equations: 

dx/dO = u, 

du/dQ =-. (eE^ / m^co^) [1 - (uTu^ ) / c^] ^/^ sin ^ . 

The sine function my be e^^almted by using two additional 
equations with appropriate initial conditions: 
dCsin Q) / dQ = cos Q or df/dQ = g 
d(cos Q) / dO - -sin 9 or dg/d9 = -^f 

The equation of motion of the electron in the gap has thus been reduced 
to the following set of four first order differential equations: 

dx/dO = u 

du/dO = (eE^ / m^u?) [l - (Ju^ / c^ ) ] 5/2 ^ 

df/d9 = g 
dg/d9 = -f 

lEi"^^22£^ill2Bg_i^^ggiP .Traject ory Equations . An 
electron is introduced into the gap at x = with negligible velocity 
V =uAi at a certain phase 9^ =uJt . The initial conditions are then. 
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for arbitrary Q^, 



X = 0, u = 0, f - sin ^Q, g = cos ^,-^. 

The effect of variation of the parameter Q^ will he studied. 

Motion of the Electron E xterl^r_to_theJ^vit^. The elec- 
tron is accelerated while in the gap and leaves the gap at x = 1 cm. 
with a velocitycOxi^ ^ and phase . It then follows a circular 
orbit at constant velocity and re-enters the gap at x = 0. The 
phase at the time of re-entry can he calculated by equation 11. i^, 
in which ^ is equal tocUxi^^^ / c. 

Re-en-ory Conditions for Gap T'rajectory Equations . The 
re-entry conditions for the gap trajectory equations are: 

^ = 0, u = u^^^, f = sin Q^^^^3L' S = ^°' ^i,n-fl' ® = ^i,n+l' 

where Q . is calculated from 9^ ^ as indicated above. 

Information to be Printed . For a given value of the 
initial phase Q^, the following information is desired for each of 

ten orbits: 

The orbit index n. 



(a 
(b 
(c 
(d 
(e 
(f 



The phase Q of the electron leaving the gap, 

e J n 

The phase 9, . of the electron re-entering the gap, 
i,n+l 

The relative velocity ^ of the electron in the orbit, 

. 2 1 -1/2 
The quantity [ 1 - p J y 

The energy of the electron in the orbit. 
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The calculation is to be repeated vith different values of the para- 
meter 0^. The parameter Q^ is printed for identification of results. 

11-^ SCALING OF THE QUANTITIES FOR ILLIAC SOLUTION . 

Magnitude of Quantities Involved . For the rnage of in- 
tegration (electron trajectory in the gap), the phase angle ^ ;^s 
known to be in tho range 1 < ^ < 2.5 from a previous desk calculator 

solution. The design constants and physical constants are E , uJ 

m-* ' 

e, m^, and c, which are, respectively, the peak electric field 
intensity in the gap, the angular frequency of the field, the charge 
of the electron, the rest mass of the electron, and the velocity of 
light. The gap length of 1 cm. is the maximi;in value of the variable 
X. The maximum value of u is easily caluclated, since oPu/c is 
the velocity of the electron relative to the velocity of light and 
is less than one. The variables f, g, and their derivatives are 
the well-known and well-behaved sine and cosine functions. The 

maximum value of the derivative du/dQ is simply the constant 

/ 2 
eE / m,u; . 

m ' J 

The ranges of the unsealed variables and values of 
the constants are indicated in Tables 11.1 and 11.2. The scaling 
of the variables is indicated in Table 11. 3. 
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SYMBOL 



u 



dx/dO 
du/dQ 

df/d9 
dg/d© 



DESCRIPTION 



Phase 



Distance 



f = cos 9 



g = sin 



dx/d9 = u 



NUMERICAL VALUES 
Min . Max . 



1 radian 





-1 
-1 



df/dO = - sin 9 -1 
dg/d9 = cos 9 -1 



2 .5 radians 



0.01 meter 

3.017 
meter/ radian 

1 

1 



O.JIT m/rad 

-o.02970ii- 0.02979^ 



SOURCE 



Previous desk 

calculator 

solution 

Design constant 



u)u < c 



sin 9 



(1- P^) 



2n3/2 



< 



Table 11.1 



Table of Variables and Derivatives 
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e/m^ 



DESCRIPTION 

Ratio of charge to 
rest mass of an 
electron 



NUMERICAL VALUE 

1.7592 X 10-^-^ 
coolomtos/kg. 



50LTE{CE 



Physical 
constant 



Velocity of light 



2.99776 X 10 
m/sec . 



8 



Pnysical 
constant 



ra 



lO 



Peak electric 
field intensity 



Angular frequency of 

the field 



5.2266 X 10'^ 
volts/meter 

(2«)(2.8)(lo9) 
rad/sec . 



Design 
constant 



Design 
constant 



X 



nmx 



Gap length 



0.01 meter 



Design 
constant 



e E 



Id 



mo) 



0.02970^ 



Derived 
constant 



K. 



6 . 51599 



K. 



0.092910 



Table 11.2 
Table of Constants (raks units ) 
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SYMBOL NUMERICAL VALUE RELATION TO UNSCALED VARIABLE 

Min . Max . 

fj 0.1 0.25 = 0.1 

y 0.001 y = 0.1 X 

p 0.017 p = u 

r 0.005 r = 0.005 f 

s 0.005 s = 0.005 g 

dy/d0 0.017 dy/d^ = dx/dQ 

dp/d0 -0.2970^ 0.2970^ dp/d0 = 10 du/dQ 

ds/d0 -0.05 0.05 ds/d0 = 0.05 df/d^ 

dr/d0 -0.05 0.05 dr/d0 = O.05 dg/d9 



Table 11. 5 
Table of Scaled Variables and Derivatives 



The Scaled Equations . The scaled equations become 

dy/d0 = p , 

dp/d0 = E [ 1/2 - u^/G ] ^^''^ (2s/0.0l) , 

dr/d0 = -10s , 

ds/d0 = lOr . 

The quantity in brackets in the equation for dp/d0 has 
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been scaled by l/2 . The constants G and E are; 

G = 2^loj = 0.00058078 ; 

E = 20 \r2~ eE /m.u)^ = 0.8^0l6 . 
m u 

The scaling for dp/d0 and for the bracketed expression 
have been absorbed by the constant E. 

Scaling for External Orbit Recycling Equation . The time 
phase of an electron re-entering the gap is calculated as follows: 

®i,n+l = K^ [ 1 - K2 (1 -a;^u2/c^)^/2 / (uAi/c ) ] [ 1 - (u>u/c f ] "^/^ 

where Yi^ = 6.51599 and K^ = O.O929IO. Since = q/io and n can be 
as great as 10, 

^i,n+l = ^\A0) (1/0.1) [ 1 - K^ (1 - u^p^/c^)^/^ / (I0u>p/c) ] 

[ 1 - (arp/cf ] "-^/^ - 2it(n + l) / 100 + 

' '^e^n 

where 0^ ^ and p are values at the time of exit after the nth tra- 
jectory of the gap. After performing the indicated subtraction, the 
quantity within the brackets <j / is on the order of 0.01. 

Additional Calculation for Printed Results . The additional 
calculations below are to be performed . 



6 =u>u/c (relative velocity) 
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done: 



U^/Uq = l/l - R (ratio of total energy to rest energy) 

T = m^c^ [ (1 - p^) -^/2 . 1 ] ^ (511. 2i^) (10^) [U^/U^ - 1 ] 
(electron kinetic energy) 

After scaling we have 
j3 =top/c , 

0.01 u/Uq = Q-Q^ ^ , 

[ 1/2 - p2/G ] ^/^ 

Q 

10"°T = 0.5112^^ [ 0.01 U /Uq - 0.01 ] . 

11.5 CODING OF THE PROBLEM . 

Organization of the Program . The following steps are to be 

(1) Read parameter from tape and print. 

(2) Set initial conditiofts for integration in gap. 

(3) Perform one stop of integration. 

(h) Test to see if last step has been performed^ if not, 
repeat step 5' 

(5) Perform external orbit recycling calculation. 

(6) Calculate and print data desired for each orbit. 

(7) Test for last orbit; if not, return to step 2. 

(8) Stop, then return to stop 1. 

Use of Library Routines . The following library routines 



were used: 
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(1) Decimal Order Input X 1 

(2 ) Constant Listing Auxiliary X 3 

(3) Differential Equations F 1 
(h-) Sine-Cosine T 1 

(5) Decimal Number Input N 5 

(6 ) Print p 1 

(7) Square Root r 1 

Use of Parameters . In order that the coding could be 
independent of the allocation of memory space, parameters in 
addition to those required by the differential equations F 1 were 
assigned for memory locations of the subroutines. The assignment 
of parameters is as follows: 

S5 location of first variable, 

S^ location of scaled derivative of the first variables, 

55 location of first word of temporary storage for 
differential equations routine, 

56 number of differential equations to be solved, 
ST location of first word of auxiliary subroutine, 

58 temporary storage, 

59 square root subroutine R 1, 

SK initial conditions setting program, 

SS main program, 

SN print routine P 1, 

SJ decimal number input routine N 5, 
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SF sine-cosine routine T 1, 

SL differential equations routine F 1. 

The parameters S3-S7 are those required for the differential 
equations routine. 

Use of Differential Equations Routine. In using the 
differential equations routine, a number of choices must be made, 
as indicated in the specification sheets . These include: 

(a) selection of increment length h and the parameter m, 

(b ) method for handling the independent variable . 

It has been noted that two first order equations are required for 
solution of the second order differential equation of motion of the 
electron in the gap. Two additional equations were included to 
evaluate the sine function required. The range of integration is 
from y = to y = 0.001 where y is one of the dependent variables. 
Since it is not possible to predict the relationship between y and 0, 
it is not possible to determine an increment h of in such a way 
that integrating over an integral number of equal increments h will 
result in y assuming its final vpilue. The technique used is one 
of integrating beyond the range desired, then integrating in the 
opposite direction, with a decreased value of h. The integration 
then oscillates around the end of the range, with decreasing values 
of h, until y lies within some predetermined interval containing its 
final value 0.001. The initial length of increment h is determined 
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from considerations of accuracy and time. The value of h chosen 
was O.JOl. For this value of h, the estimated n™,ber of steps of in- 
tegration is 150 and the accuracy Is on the order of lo"". The value 
Of m selected is the largest possible such that none of the scaled 
derivatives exceeds range, namely m = 9, or 2"h = 0.512. 

For homing on the final value y^ of the dependent variable 

y, the value of the scaling factor 2^ is changed in the following 

manner . 








i^i+2" i^i+l+h 



K^F^^ 



Figure 11. 3 
Adjustment of Interval Length 



When y^^^ - y^ differs in sign from y^ - y 2^ is changed to a new 



value 



2\- '"^^^1^1 - ^e ) 



^i ■ ^i+1 
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The integration then proceeds with the smaller interval h' (of opposite 

sign from h) until two successive values of y again lie on opposite 

sides of y . If, at any time |y - y J < 10 , the integration is 
e 1 " 

complete. 

In the differential equations of motion of the electroii 
in the gap, the independent variable j) does not appear. The value 
of is calculated by using a counter in the main routine, rather than 
integrating the equation 0' = 1. The number of equations used is thus 

Details of the Coding . The Initial Conditions Setting 
Program. Words to % of this part of the code are devoted to 
reading a value of the parameter 0^ from tape (using the decimal number 
input) and arranging the format and printing 0^ in both degrees and 
radians. (Print routine P 1 is used for printing). Since the scaling 
of is by a power of 10, printing is arranged so that the decimal 
point appears in the printed results to yield the correct values for 
the unsealed quantity 0^. 

The remaining words are used for the following: 

(l) Set initial value of orbit index n, 

(2 ) Set initial values of p and y to zero, 

(5) Set temporary storage locations for differential 
equations routine (S5 f f • ) to zero, 

{h) Calculate and set initial values of r and s. 

(The sine-cosine routine is used for the calculation), 

(5) Set initial value of 2^1, 
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(6) Set initial value of y. ., - v , 

(7) Transfer to main program. 

Details of the Coding . The Auxiliary Subroutine . The 
auxiliary subroutine is a closed routine whose function is the calcu- 
lation of the scaled derivatives from the values of the variables in 
accordance with the equations of motion of the electron. The variables 
y, p, s, and r are in locations S3, 1S5, 2S5, and 5S3, respectively. 
The scaled derivatives 2% y', etc., are to be placed in locations 
Sk to 5S^. Constanta necessary, such as G, E, and negative powers 
of ten, are obtained by using the constant listing auxiliary (Library 
Routine X 3)- Results of intermediate calcualtions are placed in 
temporary storage locations 1S8 through iksQ, for reference in checking 
the program for errors . The calculations are arranged so that none of 
the intermediate results exceeds capacity; otherwise the program itself 
is self-explanatory. The auxiliary subroutine uses the square root 
routine (entry from 6l), the first word of which is in location S9. 

Details of the Coding . The Main Program . The main program 
can be subdivided into three parts | the supervisory program for inte- 
gration in the gap, (words OL to 18l), the external orbit recycling 
calculations (words 18l to 30L), and the orbit counting and printing 
routine (words 30L to 6lL ) . 

The first part of the main program performs the following 
functions: 
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(1) The independent variable / is increased by h. 

(2) Values for the ith step of the variables y, p, 

s, r and of (y. - y ) are stored for future reference. 

(3) The i+1 step of the integration is performed by 
entering the differential equations routine; i.e., 
new values of y, p, s, and r are calculated. 

(^) y-.i - y 13 calculated and its sign is compared 

with the sign of y. - y . If the signs differ, 2"^ 

is adjusted; if the signs agree, 2 n is left 

unchanged . 

/ \ I 1-9 

(5) If |y.,T - y \ > ^0 , another step of the inte- 

1 . -9 
gration is performed. If y. - y | < 10 , the 

3tJL e 

external orbit recycling calculation is begun. 



The second part of the main program performs the calcula- 
tions indicated by the external orbit recycling equation in such a 

way that none of the intermediate results exceeds capacity. It should 

2 1/2 
be noted that the quantity [1 - p/G] ' is necessary for this calcu- 
lation. To obtain this quantity, the auxiliary subroutine is entered 
directly from the main program. The auxiliary subroutine leaves the 
desired quantity in location 3S8, where it is then available for 
further calculation by the main program. It might be thought that the 
quantity could be obtained directly without re-entry from the main 
program; it should be noted, however, that the auxiliary routine was 
last used with intermediate values of the variables . For accurate 
results, it is therefore necessary to re-enter the auxiliary subroutine 
with the final values of y, p, s, and r in locations S5 through 5S5. 
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The last part of the main program, beginning with the right- 
hand order of word 50L, is devoted to orbit counting, and to calculations 
for the format of the printed results . The print routine P 1 is entered 
in such a way that decimal points are placed to correspond to the un- 
sealed values of the results. The format used is indicated in Table 
11. 14- 
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H 
H 
I 

ro 



degrees radians 












050 


.0000 


.872665 














9 


n 


Q 








n 


radians 


degrees 


radians 


degrees 


/^ 


V^O 


T 


01. 


2.07513 


118.896 


1.36207 


078.041 


0.863718 


01.98 428 


00.5^^3205 


02, 


2.00608 


114.940 


1.37242 


078.634 


0.942114 


02.98271 


01.013640 


03. 


1.98559 


113.766 


1.40726 


08 0.630 


0.967971 


03.98364 


01.525355 


04. 


2.00941 


115.131 


1.45802 


083.538 


0.979635 


04.98157 


02.035536 


05. 


2.05493 


117.739 


1.47955 


084.772 


0.985871 


05.97197 


02.541871 


06. 


2.07354 


118.805 


1.44881 


083.011 


0.989611 


06.95861 


03.046279 


07. 


2.04100 


116.941 


1.40483 


080.491 


0.992050 


07.95137 


03.553819 


08. 


1.99582 


114.352 


1.39111 


079.705 


0.99373 


08.95100 


04.064869 


09. 


1.98125 


113.517 


1.41901 


081.303 


0.99493 


09.95242 


04.576833 


10. 


2.00854 


115.081 


1.46613 


084.003 


0.995812 


10.95028 


05.08698 



Table 11 A 
Prlnted Results 



LOCATION 


ORDER 




30 3K 


3 


00 F 




00 320F 


k 


00 F 




00 52i<-F 


5 


00 F 




00 528F 


6 


00 F 




00 kF 


7 


00 F 




00 200F 


8 


00 F 




00 5i^0F 


9 


00 F 




00 6kF 


K 


00 F 




00 22 OF 


S 


00 F 




00 250F 


N 


00 F 




00 12 5F 


J 


00 F 




00 T^F 


F 


00 F 




00 95F 


L 


00 F 




00 15 3F 



NOTES 



number of equations 



Auxiliary Subroutine 

00 200K 
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LOCATION 


ORDER 




) 


K5 F 
^2 i8l 




1 


50 1S3 


P 




TJ ISB 


2% 


2 


^0 s^ 


2%: 




50 IS 3 


p 


5 


75 IS 3 
66n oof 005 
8078 0000 J 


p 


1+ 


S9 F 

^0 2S8 




5 


32 5L 
50 5L 




6 


26 S9 
^0 3S8 




7 


50 2S8 
7J 3S8 




8 


i+0 i^s8 
50 ^S8 




9 


75 2S3 
66n OOF 00 






100 0000 0000 


J 


10 


7JN ^OF 00 5^01 




6000 0000 J 






^0 F 




11 


50 F 
75 1S8 




12 


00 IF 




13 


50 S8 
L5 3S3 





NOTES 



Set link 
2%' = 2% 



a = 1/2 - p /g 

Waste 



[ 1/2 - p^/G ] ^/^ at 5S8 



[ 1/2 - p^/G ] ^/^ at ii-S8 



0.01 



E' 



2%' = 2"^ (20Cs) E- [a]^/^ 
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LOCATION 


ORDER 




NOTES 


Ik 


66n oof 00 1000 

0000 0000 J 


0.1 






75 1S8 


2\ 




15 


ho 2Si^ 
50 S8 • 


2^s' 


= 2"^ (lOr) 


16 


LI 2S3 

66n oof 00 1000 

0000 0000 J 






17 


75 1S8 








ko 3S4 


2V' 


= 2\ (-10S) 


18 


32 18l 

22 ( )L By 0' 


Waste 





5 
if 

5 
6 



00 220K 

ifl S8 

92 149F 

92 513F 

92 96IF 

ho 5s8 

50 2L 

26 SJ 

L5 5S8 

52 75F 

50 kL 

26 SN 

L5 5S8 

50 S8 

66n i4-0F 00 729 

5779 5151 J 



Set initial conditions 



D.N.I. Input 



Print routine Print (degrees ) 



1.8/rt 
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LOCATION 


ORDER 


Vi^XCiU 




7 


S5 F 
1^0 5S8 






8 


52 TIF 
50 8l 






9 


26 SN 
92 12 9F 


Print in radians 




10 


92 513F 
19 5F 


Set n + 1 = 2 X 2"'^ 




11 


kO 6s8 








i^-1 1S5 


Set p to 




12 


i^l S5 
i^-1 S5 


Set y to 




13 


kl 1S5 
ill 2S5 






li^ 


ill 5S5 








L5 5S8 







15 


50 S8 
00 2F 






16 


Li| 5S8 5 in A 
50 i6l 






IT 


26 SF -* 


Sine cosine routine 


Set 




kO TS8 




starting 


18 


TJN OOF 00 100 




values 




0000 0000 J 


0.01 






kO 2S5 


s = 1/200 sin 10 




19 


50 TS8 
TJN OOF 00 








100 0000 0000 J 


0.01 




20 


^0 5S3 

L5N OOF 00 5120 


r = 1/200 cos 10 






0000 0000 J 


2% 





11-2T 



LOCATION 


ORDER 


21 


kO 1S8 




LIN OOF 00 10 




0000 0000 J 


22 


40 8S8 




26 SS 



NOTES 



Set 2% 



0.001 



Set (y^ - y^) 



MAIN ROUTINE 




00 250K 





L5 1S8 




10 9F 


1 


Lk 5S8 




ho 5S8 


2 


L5 S3 




4o i^-S5 


3 


L5 1S5 




40 5S5 


i^ 


L5 2S3 




40 6S5 


5 


L5 3S5 




40 TS5 


6 


L5 8s8 




40 9S8 


7 


00 9F 




50 7L 


8 


26 SL 




L5 S5 


9 


LON OOF 00 . 




OOOD 0000 J 




40 8S8 



2^ 
h in A 

Replace by /S + h 



Store old values of y, p, s, r 



Store (y^ - y^ ) 



Differential equations Routine F 1 



y = 0.001 
e 
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LOCATION 


ORDER 






10 


36 12L 
L5 9S8 






11 


36 13L 
26 17L 






12 


L5 9S8 

36 17L 






13 


L5 ^S5 




^i 




LO S5 




^i+l 


li+ 
15 


kO 1DS8 
50 8S8 
75 1S8 
66 10S8 




2% 


16 


S5 F 
kO 1S8 






17 


L7 8s8 








LON OOF 


00 






1000 J 






18 


36 L 
50 18L 






19 


26 S7 




—> 




L5N OOF 


00 kk6 






6080 0000 J 




20 


50 S8 
66 3S8 






21 


S5 F 
kO 11S8 






22 


LIN OOF 


00 






10000 0000 J 






50 S8 






23 


66 1S3 

S5 F 







NOTES 



If y^ - Yg and y^^^ - y^ agree in 



sign, integrate again 



If y^ - yg and y^^^ - y^ disagree 
in sign, form 2 n' 



2^. = 2^ (y.^^ - ye) / (yi - yi^i^ 



Is 



V - y 

l-^i -^e 



-9 
< 10 ^? 



Auxiliary subroutine 



M 



10 



-3 
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LOCATION 


ORDER 




2i^ 


^0 12S8 

50 6S8 (n+l)2"'^ 


-10"^ / P 


25 


7JN OOF 00 1000 






0000 0000 J 


0.1 




ifO F 




26 


50 F 

75N 40F 00 1285 






1800 0000 J 


2Tt/lJ 


27 


00 7F 
10 11S8 




28 


LO 12 S8 

66n oof 00 1000 






0000 0000 J 


0.1 


29 


SI F 




30 


Lk 5S8 
40 15S8 
92 12 9F 



^i,n+l 


51 


92 515F 
L5 6s8 




52 


LON OIF OOF 
10 52F 


n X 2~^ 


55 


JO 22F 
50 55L 




3h 


26 SN 
L5 5S8 


Print n 


55 


52 6IF 
50 55L 




56 


26 SN 

5ON i<-OF^OO 729 

5779 5151 J 


Print 

^e,n 



NOTES 



in raaians 
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LOCATION 

37 

59 

i+0 
Ul 

^2 

hh 
hi 



1+8 
1^9 



ORDER 

7J 5S8 

22 38L 

52 65F 

50 38L 

26 SN 

L5 13S8 

i^0 5S8 

22 ^IL 

52 6IF 

50 i^■lL 

26 SN 

50N ^OF 00729 

5779 5151 J 

7J 13S8 

22 ^^L 



Waste 



Waste 



52 63F 

50 i+^L 

26 SN 

56 ^6l Waste 

L5 1S5 

50 S8 

66n oof 00 170 

il-100 0000 J 
7JN OOF 00 1000 
0000 0000 J 

52 71F 
50 ^8l 

26 SN 

L5N OOF 00 70 

7106 7812 J 



NOTES 



Waste 

Print in degrees 
'^e,n 



Set for next integration 



Print 0. ^ in radians 
1.8/n 



Print 0^ ^^^ in degrees 



Calculate and print 
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LOCATION 


ORDER 


50 


50 S8 




66 3S8 


51 


S5 F 




^0 I4s8 


52 


52 72F 




50 52L 


53 


26 SN 




L5 1^S8 


5^ 


LON OOF 00 100 




0000 0000 J 




ho F 


55 


50N ^OF 30 112 




4000 0000 J 




7 J F 


56 


52 82F 




50 56l 


57 


26 SN 




L5 6s8 


58 


L^N OIF OOF 




ho 6SQ 


59 


LON ONF OOF 




54 SK 


60 


26 12SK 




00 F 



NOTES 



Calculate and print 



I 1 



- p ] 



-1/2 



CalcLialte and print energy (mev . ) 



Increase ri 



If n -^ IJ. s^.op, read in ne-, 
If n < 11, integrate again 



Temporary Storage Starting at Memoi-y Location ^ho 




2\ 
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J'.OCA.TION ORDER NOTES 

2 1/2 - p^/G 

[ 1/2 - p2/G ]^/2 

]| [ 1/2 - p^/G ]^/^ 

6 (n+1) 2"'^ 

7 1/2 cos 10 

6 (y^^i - y^ ) 

9 (y^ - ye) 

10 . y^ - y^^i 

11 M [ 1/2 - p^/G 3 -^/^ 

l;2 -10'^ / p 



1, 



2.k [ 1 - P^ ] "^/^ 



Mlocatlon of Memory Space . Memory space was allocated 

as f ulJ,c''"3 : 

j-l'^ para.meters 

1.6 --bj available for list of constants for 

constant-listing auxiliary 

6h~'J'5 sqimre root routine R 1 

■jk'<^k decimal number input routine N 5 

95-] 2*4- sine-cosine routine T 1 

125-152 print routine P 1 
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153-193 differential equations routine F 1 

194-199 unused 

200-218 auxiliary routine 

219 unused 

22O-2I+2 initial conditions setting program 

243-2i)-9 unused 

25O-311 main program 

320-323 variables for integration routine 

324-527 derivatives for integration routine 

328-331 temporary storage for integration routine 

332-335 previous values of variables for integration 
routine 

336-339 unused 

3^0-35^ temporary storage 

355-9^0 unused 

9^1-961 constant-listing auxiliary X 3 

962-998 available for post mortem D.Q.I. (Librarv 
Routine C 1) ^ J 

999-1023 decimal order input X 1 

^"2 temporary storage 

Tape Preparation. The tape was prepared in two steps . The 
non-standard parts of the program - namely the list of parameters, 
the auxiliary subroutine, the initial conditions setting program, and 
the main prograon — - were punched, printed, and visually checked for 
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errors . Us ing library routines , the tape was prepared as follows 

Decimal order input X 1 

00 9^1K 

Constant listing auxiliary X 3 

00 5K 

List of parameters 

00 6kK 

Square root subroutine R 1 

00 T^K 

Decimal number input routine N 3 

00 95K 

Sine -cosine routine T 1 

00 125K 

Print routine P 1 

00 155K 

Differential equations routine F 1 

00 200K 

Auxiliary routine 

00 220K 

Initial conditions setting program 

00 250K 

Main program 
2k 999N 
26 220N 
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ko6n 

K055N 

K05N 

The tape was then checked for reperf orating errors, using a tape 
comparer. The directive 2k 999N vas placed on the tape so that a 
correction tape could be placed in the reader before the program was 
started. Initially, however, the stop of the 2h 999N directive was 
by-passed, and control was transferred to the initial conditions 
setting program by the directive 26 220N. The quantities K06n, 
K055N, and KO5N at the end of the tape are values of the parameter 
0Q to be read from tape by the decimal number input. 

11-6 CHECKING THE PROGRAM. The electron cyclotron pro- 
gram was typical in that a number of blunders in coding were made. 
A description of the sequence of events in checking the program follows 

When the program was first read into the Illiac, the 
quantity 

07680. 0000 00528 59832 06685 86695 k9Q3 
was punched continually. The temporary storage was printed (P.M. 
routine C 5) indicating that 0^ had been read from tape and stored 
(at 5S8). 0Q was to be initially expressed in degrees scaled by 
100 and was later to be converted to radians scaled by 10. Since the 
conversion had not occurred, the difficulty was isolated to the 
region 5L-7L of the initial conditions setting program. A subsequent 
printing of the orders of the print routine (with P.M. routine C 5) 
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Indicated an incorrect link and the difficulty was isolated to the 
print routine entry (word ^L ) which was 

52 75F 50 2L rather than 52 7,:^ 50 ^L. 
When the program was read into Illiac again, noises in- 
dicative of integration were heard, and after an appropriate interval 
of time punching occurred. The results were, however, not entirely 

correct. Values of ^ and ^. were incorrect, although some small 

e,x 1^^ 

comfort was gained from the fact that Q^ and n were printed correctly. 



The behavior of the variables /), y, p, s and r for the first ten 
steps of integration were then observed by using check point routine 
D 1 with a blocking order placed at location 256 (6l in the main routine ) 
ahead of the entry to the Runge-Kutta routine. The initial values of 
0f Yi V) s, and r were correct, exonerating for the moment the initial 
conditions setting program, but behaved peculiarly on successive 
steps. Rather than increasing uniformly by increments of O.OJl, 
increased by 0,301, then by 0.02^4-6, by 0.0?j06, and eventually de- 
creased slightly. It was finally found that the location 5S8 used 
for storage of was also being used for temporary storage by the 
auxiliary routine. This blunder was cured by replacing the orders 
^0 5S8 50 5S8 in words 10 and IIL of the auxiliary routine by 
hOY 50F. 

With a somewhat longer correction tape, the program was 
again checked on the Illiac . The results of the first integration and 
orbit cycle were correct, but the results for the second orbit were 
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all incorrect with the exception of the orbit index n. The re-entry 
into the initial conditions resetting program was checked and it was 
found that p was incorrectly being reset to zero; i.e., the program 
was robbing the electron of all the velocity it had gained during its 
first orbit. The correction tape increased in length; word 6lL of the 
main program became 26 12SK 00 F rather than 22 IISK DO F. The 
orders for clearing q^ and p were interchanged in words IIL and li^L 
of the initial conditions setting program. 

After the following code check, the results for the first 
orbit were again incorrect. After mutterings of "something is wrong 
with the computer", it was discovered that a terminating symbol had 
been omitted on th^^ correction tape. The correction process then 
converged and correct results were obtained. 

11.7 CONCLUSION. The preparation of a problem for 
solution on a digital computer is by no means a completely objective 
process. No two programmers would prepare the same problem for solu- 
tion in an identical way. Furthermore, a second coding of the problem 
by the same programmer would differ from the first preparation. 

The personal preferences of the programmer affected the 
electron cyclotron program in a number of ways. Memory space other 
than location 0, 1, and 2 was assigned for temporary storage; it 
was felt that such an arrangement might aid in code checking. The 
S terminating symbol was used for designation of the locations of 
first words of subroutines, so that the detailed coding could be 
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completed before memory locations were assigned to the subroutines. 
Thus, the location of a subroutine affects only one of the S para- 
meters rather than the addresses of a number of orders in the program 

Among the changes which might be made if the program were 
rewritten is a change in the handling of the variable 9. If were 
expressed in revolutions, the overflow properties of the computer 
would simplify the recycling computation. 
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CHAPTER 12 



THE CATHODE RAY TUBE DISPLAY 



12.1 GENERAL DESCRIPTION , Results of certain types of 
calculations may "be most conveniently obtained from the computer by 
use of the catlode ray tube display. Whenever results of a cal- 
culation can be represented pictorially as a graph or diagram, the 
cathode ray tube output from the computer provides a rapid and 
elegant method for obtaining these results. If results would other- 
wise be graphed manually the cathode ray tube saves the human time 
required for this additional processing and also reduces the machine 
time required to present the results to the user. 

Another use of the cathode ray tube is to provide inter- 
mediate results in such a form that they may be analyzed subjectively 
during the course of a program. In this way the programmer may 
discover errors in his program at an early stage in the calculation. 
He may also obtain immediate results which may be used to determine 
the way the remainder of the program is to be run. Even if the 
final results must be presented with greater accuracy than is 
possible using a graphical display, some advantage may be obtained 
by supplementing the digital results with a graph or diagram. The 
cathode ray tube output is called into use by turning the output 
switch on the input-output rack to the cathode ray tube position. 
When the switch is in this position the cathode ray tube is the only 
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form of output available to the user. The display on the cathode 
ray tu^e is limited to a 3" x 3" square area centered on the face 
Of the tube. Within this square area, chosen points may be brightened 
on a 256 x 256 regular square raster. A finer raster is unnecessary 
because of the limited resolution of the cathode ray tube. This 
limited resolution makes the display of digital information such 
as numbers or figures somewhat unhandy. It has been found that no 
more than about 200 fully legible numbers or letters may be dis- 
played in a single frame. In general the requirements of the format 
will restrict it to somewhat less than this. This disadvantage, 
however, is somewhat offset by the possibility of using characters 
having any shape which the programmer desires and by the greater 
speed of output. A subroutine has been written which will display 
characters at the rate of about h^ per second. 

In order to retain a permanent record of displayed re- 
sults it is necessary to photograph the face of the cathode ray 
tube. Two cathode ray tubes are installed in the rack. They 
are driven so they operate in parallel and hence display the same 
results. One is used for visual observation and the other is 
equipped with a semi-automatic camera to permit photographing 
the display. The film advance mechanism is controlled by the 
program while insertion, removal, and development of film must 
be done by hand. 
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12.2 ORDERS COIJTROLLING THE CATHODE RAY TUBE DISPLAY . 
Orders which would normally cause characters to he punched or printed 
will affect the cathode ray tube when the output switch is in the 
cathode ray tube position. The conventional output order to use for 
cathode ray tube display is 82 16F. When this order is executed 
it will cause one spot to be brightened on the face of the cathode 
ray tube. The position of the brightened spot depends upon the 
contents of the l6 leftmost binary digits 2° - 2'-^^ of the accumulator. 
Let these digits be designated by the symbols a^, a^, . . . , a^^. 
Assume an origin or coordinates in the lower left-hand corner of 
the square raster and let d represent the length of one side of 
the square raster (d =3"). Then the coordinates of the brightened 
spot will be: 

y = d (a^ 2"-^ + a^ 2" + . . , + a^ 2" ) 

o 
X = d (ag 2"-^ + a^ 2" + . . . + a^^ 2.' ) 

These formulae say, in effect, that the ordinate is obtained from 
the first eight binary digits of A, regarded as a number, and that 
the abscissa is obtained in the same way from the second eight digits. 

During the execution of the 82 16F order the AQ register 
will suffer a left shift of l6 places. This output order takes 800 
microseconds. A variation of this order permits a spot to be 
brightened in only ^^00 microseconds . When an 82 8F order is executed 
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the ordinate Is determined in the same manner as it is for an 82 16f 

"LP o 

order, y = dCa^ 2 + a^ 2 + . . . + a 2* ), but the abcissa vill 

be the same as that of the spot produced by the last 82 16F order 
which was executed. In this way time may be saved when several points 
having the same abscissa are displayed. 

A letter output order is used for operating the semi-automatic 
camera. The following order is used for this purpose: 

92 769 Advance the film one frame. 

The camera has no shutter, so that programs having long 
periods of computation (say I5 minutes) between displays of results 
should be arranged to advance the film before and after each display 
in order that the results will not be fogged from filament light in 
the cathode ray tube. 

Exposure of a single frame of the film takes place during 
the period of time that the spots that make up the picture are 
brightened. Amount of exposure is therefore controlled by the 
Intensity setting of the cathode ray tube. Exposure of a spot may 
also be increased by multiple brightenings . If calculations are 
made while the frame is being exposed the resulting photograph will 
in no way be affected although the image on the visual cathode ray 
tube may fade partially before the completion of the frame. Ad- 
vancing the film one frame takes one second and the Illiac will 
wait until this operation is completed before executing other out- 
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put orders . If the orders f olloving the film advance do not in- 
volve output, their execution will be begun after only 200 microseconds 

No output orders other than those mentioned are of practical 
use in the operation of the cathode ray tube display. A 92 order 
having an odd address will affect the camera mechanism, while any 
other output order will brighten a spot. 

12.5 PROGRAMfl:NG FOR TH E CATHODE RAY TUBE DISHAY. As 
with typewritten output, the cathode ray tube display is usually 
programmed by means of subroutines. If digital resiats are desired 
the number of orders concerned with output will be greater when the 
cathode ray tube display is used than when the pimch is used. In- 
formation needed to form a character (such as a letter or decimal 
digit) may be stored in a single word in the memory, but the 
mechanism required to decipher this information rapidly may take 
as many as 55 words. Such a device is used in library program 
CO 1. This program uses 65 words to display information similar 
to that printed by library program C 5 which has only 51 words. 

Graphical display is also most conveniently accomplished 
by means of subroutines . Program 1 is an example of a subroutine 
which enables one to plot axes and points on the cathode ray tube. 
Any method of graphical display will almost inevitably require 
scaling of the coordinates of the points to be displayed. Since 
it is desirable to utilize the full screen area of the cathode 
ray tube, up scaling may be necessary in some cases instead of down 
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scaling. The full screen is regarded as having a range of 1 for 
both coordinates when Library Routine 1 is used. The full range 
of all variables in this case should be scaled so as to correspond 
to 1. If more than one graph is to be displayed on a single frame, 
however, it is sometimes better to avoid entering the point plotting 
subroutines more than once. Consider the following examples The 
points (x^, y^) and (x^, y^ ) are both to be displayed relative to 
the origin. Let us assume that axes have already been displayed 
by program 1. The quantities x^ and y^ have been supplied to 
program 1 when the axes were plotted so as to locate the origin. 
They represent the x and y coordinates respectively of the center 
of the screen in the coordinate system used. Let the necessary 
parameters be those given in Table 12.1: 
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LOCATION 


PAEAMETKH 


10 


^0 


11 


^0 


12 


^1 


13 


^1 


1^ 


72 


15 


scalin 


16 


scalin 


IT 


1/2 



Assume these coordinates have been 
properly scaled before being stored. 



Table 12.1 
Parameters for Cathode Ray Tube Display 

A program to plot the tvo points might be the one given in Table 
12.2: 

p 50 13 scale y^ 

TJ 15 
P+-L JO 12 plot (x^, y^) by use of subroutine 1 

50 p+1 
p+2 26 (to program l) 

50 l4 scale y^ 

p+3 TJ 16 

Lit- IT prepare 2(y2 + l/2 - Yq) 

^+k LO 11 to be displayed 

00 1 
p^5 82 8 plot (x^, y^) without use of a 

subroutine 

Table 12.2 
Program for Cathode Ray Tube Display 
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In the program of Table 12.2 the second ordimte y^ is first scaled 
and then translated so as to be present in the accumulator in the 
proper form for display purposes. It is then displayed by use of the 
82 8 order, thus saving machine time and also saving the programmer 
the trouble of re -introducing the abscissa x . 
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CHAPTER 13 
PROGRAMMED CHECKS 

13.1 GENERAL PRINCIPLES OF CHECKING . Any person using 
the Illiac as an aid to his research is likely to question the 
reliability of results of calculations carried out by the machine. 
Although the Illiac is quite reliable in its operation, one is 
almost never Justified in using Illiac results directly for publi- 
cation or for further research without assuring oneself in some 
way of their correctness. If the calculations carried out by the 
computer are to be treated with the same care and thoroughness 
that is usually required of results obtained in other ways some 
method of checking must be found. Naturally the Computer Laboratory 
cannot guarantee the correct operation of the computer and it is 
the responsibility of the person preparing the problem to provide 
adequate checking of his results. 

The subject of checking is by its nature, a difficult 
one to treat in generality since the most suitable checking method 
usually depends upon the nature of the calculations. An obvious 
checking procedure is available in some cases and not in others . 
For example, one may check the solution of a set of algebraic 
equations by substituting the results into the original equations. 
Whenever a checking process of this sort exists it is preferable 
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to checking methodo which only check some of the operations Involved 
m the calculation since It provides an overall check on the operation 
of the computer from the beginning to the end of the problem. If the 
results Of a calculation may be plotted as a graph one may often verify 
the correctness of the results by noting whether the points follow a 
smooth curve. Such a simple check as this may be used only when the 
various points are calculated independently and should then be used In 
conjunction with other checking processes. Some types of calculation 
may be closely approximated by much simpler formulae which may be 
evaluated independently. -When such approximations exist they are 
Often adequate for checking purposes. A checking process must thus 
be tailored to fit the problem and its effectiveness will often depend 
upon the ingenuity of the programmer. 

Let us assume that a program is divided into distinct runs 
Of length t^ and that a checking period of length t^ follows each 
run. If a malfunction or machine error occurs either during the run 
or during the checking period It will be detected at the end of the 
checking period and the entire run will have to be repeated. The 
average time wasted due to errors and checking and the optimum length 
Of run may be calculated using these assumptions. Take the time T as 
the mean time between machine errors, i.e. the reciprocal of the error 
frequency. It is not strictly correct to speak of a fixed time T 
alnce the error frequency depends upon the nature of the calculation 
and other time dependent factors, nevertheless, we shall assume a 



13-2 



fixed T and show later that our conclusions do not depend upon an 
accurate knowledge of T. A reasonable choice might be T = 8 hours. 
The average useful time per run may be taken as the prob- 
ability that no error will occur during the run multiplied by the 
length of run. Assuming a Poisson distribution of errors one obtains 

Average useful time per run = t exp[-(t + t )/T] (l) 

Since the total time taken by the run is t^ + t we find that the 
fraction of the total time which is useful is 



The first factor differs from 1 because of the time taken by the 
checking process while the second factor differs from 1 because 
of the likelihood of errors occurring during the run. If we 
maximize- (2 ) with respect to t we obtain 

^1 ■*■ V2 - V = ^' (3) 

giving 

t^/t^ =Vl/^ + T/t^ -1/2, (k) 

If we as Slime that T » t^ we obtain the approximation 

\=VV"- (5) 
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If the same assumption is made in formula (2 ) we obtain for the 
fraction of the total time which is useful 

F = 1 - t^/\ - t^/T . (6) 

Substituting (5 ) in (6 ) we obtain 

F = 1 - 2VyY (7) 

as the largest fraction obtainable. It is important to notice 
that if T is much larger than t^ formula (6) will remain nearly 
1 over a wide range of t^. In effect, this means that our selec- 
tion of t^ is not critical provided it is much more than t and 
much less than T. For example, if the time for checking, t , is 
5 minutes we obtain for the optimum length of run 

t-j_ =V5 min. x 480 min. = h9 minutes . 

Yet as seen from formula (6), if the length of run were doubled 
or halved the wasted time 1 - F would not be much greater. 

In formula (l) if t^ is allowed to increase we see that 
the function reaches a maximum when t = T. Therefore if a run 
is permitted to exceed T the useful time achieved by the run (on 
the average) begins to decrease. Under no circumstances will any 
advantage be gained, therefore, by making t greater than T. From 
(h) we see that this limit is approached as t increases. 



13-^ 



An alternative system of checking may be used in which 
the checking process is incorporated in the calculation and the 
program is arranged to stop whenever an error occurs . The run must 
be restarted from the beginning after the error. Each successful 
run (of length t ) will be followed by a period of printing or storing 
intermediate results which will allow the following run to be re- 
st8.rted at that point. We shall assume this period to be of length t 
The difference between this method and the preceding one is that an 
unsuccessful run is cut short in this case. Although formula (l) is 
unaltered the ones which follow are changed somewhat. This is be- 
cause the total time taken by a run (on the average) is no longer 



t^ + t^ but is rather 

1 d 



^1 2 



(t^ + t^) exp[.(t^ + t^VT] + -y- r U exp[-(u/T) ] dt 

= T ^1 - exp[-(t^ + t^VT]]^ 

The first term on the left is due to successful runs while the second 
term is due to unsuccessful runs . The formulae corresponding to the 
previous analysis are 

^= T(exp[(t^ + t^VTj -ly ' ^2'^ 

t^ - t{i - exp[-(t^ + t2)/T ]} = , (3' ) 

t^ =-^t^ , (5' ) 

F = 1 - -^ — -^ , (6-) 

t 2 T ^ ^ 
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F = 1 -V 2t /t" 



(7-) 



Here we see that the optimum length of run is increased by a factor 
ofVF, and that the wasted time is decreased by the same factor. 
Although no closed expression can be obtained corresponding to (k) 
we note that the remarks concerning t^ in the neighborhood of T 
still apply. Qualitatively the behavior of F as a function of t 
is still the same so F remains nearly unity when t « t « T. 
Thus a careful choice of t^ and an exact knowledge of T are still 
not important. 

13.2 INPUT CHECKING. Often, an overall check on the 
entire program is not possible and individual checks must be made 
on separate operations. The operations of reading and punching 
are more subject to error than other operations and are also easier 
to check in most cases. Standard techniques are also available 
for checking input and output operations so that relatively little 
burden is placed on the programmer. 

The sum check is commonly used to check input-output 
operations. It derives its name from the fact that a primitive 
form of the sum check consists of merely adding together all the 
words in the memory before the calculation begins, but after the 
program has been read into the computer. If this sum agrees with 
a previously obtained sum, one concludes that the program has been 
placed in the memory in exactly the same way as it was when the 
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sum was originally obtained. Agreement between the sums is taken 
to mean that no error was committed during the reading of the pro- 
gram. 

A sum check is particularly useful for checking the reading 
of a master tape which is to be used many times. Such tapes eventually 
wear out and a sum check provides an indication of this as well as a 
check on the operation of the computer. Library Program X 7 is an 
interlude which computes the sum of the memory and compares the 
result with a sum which is especially prepared on tape the first 
time the tape is used. A somewhat more reliable check than would 
be obtained by simply summing is obtained in this program by per- 
forming a circular right shift of the accumulated sum during the 
summing process. By this technique one may avoid the possibility 
that two errors occurring in the same digit position, in different 
words may compensate, thus yielding a correct sum. Misplacement 
of words is also detected by this technique. A typical interlude 
to sum the memory might be the one in Table 13-1. 
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00 


F -1 
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26 



Place accumulated sum in A and 
perform circular right shift. 

Add next word to sum and restore 
in zero. 

Step address to select next word. 

Test for the end of the summing 
process . 

Read previous memory sum from tape 
and compare with memory sum obtained 
this time. 

If sums do not agree stop on an FF 
order; othearwise return to the D.O.I. 

Test constant. 
Enter interlude. 
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■Table 13.I 
Sum Memory Interlude 

After computing the "memory sum" this Interlude vUl read a previously 
prepared word from the tape which It compares with the memory sum. If 
the two disagree the routine will stop on an F^ order. This order has 
heen specially chosen to indicate the failure of programmed checks. It 
Should be used whenever a programmed check is made Instead of the OF 
stop so the computer operator will be able to Identify the nature of 
the stop. 
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So far we have discussed a method for insuring the correct 
reading of a program tape or master tape. A data tape which is to be 
used many times may be treated in the same way. Tapes which are pre- 
pared by the computer constitute a special problem and may be best 
discussed in connection with output checking. There still remains the 
problem of checking the reading of data tapes which are prepared by 
hand and which need to be read only once. No entirely satisfactory 
method for checking such tapes has been developed which does not 
involve reading the tape twice. If arrangements are made for the 
master program to read a tape twice one may store the numbers the 
first time they are read and then compare them directly on the second 
reading. If memory space does not permit storing the entire set of 
numbers which are to be read one may effect an equally valid check 
by computing a sum of the numbers as they are read the first time 
and retaining only this sum to compare with a similar sum computed 

on the second reading. 

Other ways of checking data tapes involve special preparation 
of the tape. One such technique is to follow the punching of each 
number by the sum of its digits modulo 9- These sums must be com- 
puted by hand as the tape is prepared thus complicating the preparation 
of the tape and increasing the likelihood of human errors. Ae the 
tape is read the corresponding sums may be formed by the computer (a 
special input routine is required) and the reading of the tape checked. 
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A somewhat simpler technique requires the use of a 
code in the punching of characters. Punchir^ l„ code may be 
simplified If one covers the Teletype keys with caps which 
indicate the coded characters while the tape is being punched. 
A convenient code uses the 5th hole of a character as a parity 
check on the other four; for example, the first four holes might 
be used to indicate a sexadecimal character as they are at present 
and the 5th hole would he punched or not punched so as to make the 
number of holes per character odd. A special Input routine would 
be used to interpret the coded digits and check for errors. An 
Odd number of errors occurring in the reading of a single charac- 
ter would he detected by means of such a check since the parity 
would then be altered. 

In order to obtain a typewritten copy of the coded tape 
It is necessary to clamp the 5th pin on the transmitter-distributor 
unit Of the Teletype in the down position so as to suppress the effect 
Of the 5th hole. Thus all carriage returns must be Inserted manually 
as the typewritten copy is formed. Another drawback to this method 
Of Checking is that a character may be Ignored by the reader or read 
as two Identical characters without an indication of error. In 
order to check that the correct number of characters has been read 
one must Include some additional restriction on the typing of the 
data. Such a restriction might be that all numbers which are typed 
Should have an even number of digits. Ifon-slgnlf leant digits would 
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then have to he included In some numbers thus increasing the time 
spent typing as well as the time required to pass the tape through 
the reader. Still the time would be less than that required to read 
the tape twice. 

13*3 OUTPUT CHECKING . The checking of punched results 
is usually a simple operation compared to the checking of the reading 
of hand-prepared tapes . A standard technique is to arrange to read 
the output tape back into the computer at the completion of the problem. 
If the punched results are retained in the memory they may be compared 
directly at this time. If it is not possible to store the punched 
results one may compute a sum of all results as they are punched and 
then recompute this sum as the output tape is read back into the 
computer for checking purposes. 

If punched output from the computer is intended to be read 
back into the computer at a later time it is convenient to include a 
sum of the punched results on the output tape. This sum may be com- 
puted during the punching process and finally punched at the end of 
the tape. When this tape is read back again it may be checked for 
punch or reader errors by recomputing the sum as it is read and com- 
paring this sum with the sum punched at the end of the tape. 

No convenient method exists for checking the Teletypewriter 
or the cathode ray tube output. This is because no input device exists 
which allows one to read printed or photographed results into the 
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computer. Three methods exist for checking such results but all 
three involve considerable human labor which must be done without 
error if a check is to be made. These methods are. 

1. Print the results twice from a checked tape 
and compare the two printed copies by eye. 

2. Copy the printed results by hand onto Tele- 
type tape and then use the sum check method 
of the previous paragraph to check this tape, 
treating it as though it were computer output. 

3. Include some redundancy in the printed results 
so that the output may be checked by hand 
computation. One might, for example, follow 
the printing of each number by the sum of its 
digits modulo 9. This computation could then 
be made by hand using the printed results and 
could be compared. 

13.4 MEMORY CHECKIMG. As we have seen In the previous 
sections of this chapter, the most satisfactory check of a calculation 
is an overall check which checks all operations performed during the 
calculation simultaneously. When such checks are impossible, however, 
one must resort to the use of a combination of partial checks. One 
then Wishes to check the operation of the input-output equipment, 
the memory, and the arithmetic unit by means of Independent checks. 
The Checking of the high-speed memory is also usually treated separately 
from the checking of the drum memory. 
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During a long calculation it sometimes happens that words 
in the high-speed memory become altered for no apparent reason. If 
the altered word happens to be an order pair in the program the result 
will usually be a malfunction of the program which in many cases will 
not be difficult to detect. Yet there are many memory failures which 
can produce subtle changes in the results, and for this reason a memory 
check must be made. Some methods of arithmetic checking are complete 
enough to include the memory and no special check of the memory is 
needed, but if this is not the case some special method must be devised. 

Library program X 10 is a general subroutine for checking 
all parts of the memory which are not altered after the first few 
minutes of calculation. Provision is also made for removing the 
program being checked from the computer each time such a check is 
made. Complete input checking is also incorporated, in this program 
as an additional feature. This versatile subroutine may either be 
arranged so that the removal of the program from the machine is under 
the control of the operator or it may occur at intervals determined by 
the program. When removal is to occur the master program must enter 
X 10 as a closed subroutine. It causes a tape to be punched con- 
taining information concerning the subroutine link, those words in 
the memory which may suffer modification, and finally the memory sum. 
A special tape is read at this time to specify the words which may 
be modified. After punching is completed the program stops on an OF 
order. To resume the program at a later time one merely reads the 
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initial tapes as if the program were being started from the 
beginning and then arrests its progress with the white switch. 
Using a bootstrap start the tape which was punched just prior 
to removal may be read. All words will be restored to their 
condition at the time of entry to X 10 and a sum check performed, 
If the sum check agrees with that punched previously one is 
assured that words have been read correctly either at this 
reading or the previous one, that the specially punched tape 
was punched correctly, and that no word in the memory suffered 
a change during the previous run with the exception of those 
words specially designated. These latter words are usually 
checked by the arithmetic check. 

Checking of the drum memory is achieved more easily 
than checking of the high-speed memory. In most applications of 
the drum memory words are transferred to and from the drum in 
blocks of perhaps ten or more. The standard block transfer 
routine effects a check on the drum by recording a sum check 
following each block which is recorded on the drum. Upon play- 
back from the drum a similar sum is computed to be compared 
with the one which was recorded. In this way the operation of 
the drum is completely checked without special attention on the 
part of the programmer. Only when individual words are to be 
transferred to and from the drum must special arrangements be 
made for checking. 



13-1^ 



15.5 AEITHMBTIC CHECKING. In most cases a characteristic 
identity of the problem "being solved is used for arithmetic checking 
purposes . Two examples are presented here to show how arithmetic 
checking is carried out. 

a) kss\me that n values of cos -^ k with k = 0, 1^ ..., n - 1 

are to he used in a calculation and that 2. is non-integral . Checking 

n-1 2rtni 
may he accomplished by means of the identity ^^ cos -^ k = 0. As 

the cosine functions are computed they may be summed and the final 

result tested to check that it is sufficiently close to zero. 

b ) Assume that the n products ax^ are to be formed where 

a is a given number and x^ is given for i = 0, 1, ..., n - 1. These 

n-1 n-1 
products may be checked by computing Xq x^ and ^Eq ax^ so as 

n-1 n-1 
to allow a check by the identity ^Zq ax^ = a^Z^ x^ . 

In case a convenient identity cannot be found one may 
sometimes avoid duplicating all calculations by using an approximate 
calculation as a check. An integration process for example, may be 
checked by carrying out the same process using a larger step (or 
mesh size) during the checking run. The checking run will thus take 
considerably less time than the original, yet it will detect any 
malfunction which produces a gross error. 

15'^ SUMMARY. Although the operation of the Illiac is 
reliable by computer standards one is rarely justified in using Illiac 
results for scientific purposes without checking them in some way. 
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Sometimes It is possible to check results visually or by hand, but 
when this is not possible one must maJ.e use of a programmed checking 
technique. As has been seen the simplest and most effective pro- 
grammed checks are those which check the solution of the problem 
as a whole rather than the individual operations. Such overall 
checks are possible when the solution may be substituted back into 
an equation or when some fundamental redundancy exists in the results 

Even when an overall check is not possible one may obtain 
a valid check by combining checks on input-output, memory and arith- 
metic units. Although such individual checks are often difficult to 
perform they are always preferable to running the problem twice and 
thereby doubling the time taken. A summary of methods of checking 
Individual operations is given in Table I3.2. 
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Operation Checked 



Checking Methods 



1. Reading of program tape or 
reading of data tape which 
is to be used several times 

2. Reading of computer output. 



3. Reading of data tape to be 
used once. 



h. Punching of output tape 



5. Printing of results. 



6. Storage in high-speed memory. 



7. Transfer to and from the drum 
and storage on the drum. 

8. Arithmetic operations. 



Sum check such as Library 
Routine X 7 or X 10. 



Sum check to be punched on 
computer output tape. 

a. Read tape twice and arrange 
that one reading shall check 
the other. 

b. Use a parity check code on 
characters which are read. 

c. Use a modulo check on the 
digits of each number. 

a. Compare with stored results 
by reading punched results. 

b. Compute sum check as punched 
tape is read and compare 
with previously computed and 
stored or punched sum check. 

a. Print twice and compare by 
visual inspection. 

b. Prepare tape by hand from 
printed page and check as 

in i^-b. 

c. Use a modulo check on printed 
numbers and compute by hand. 

a. Use a sum check routine such 
as X 10. 

b. Combine with arithmetic checks 

Use a sum check such as the one 
included in N 7* 

a. Use an identity which relates 
the various quantities computed 

b. Use an approximate method which 
is faster than the method being 
used in the actual calculation. 



Table 13-2 
Checking Methods for Various Operations 
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CHAPTER l4 



THE PROGRAM LIBRARY 



The problem of planning and carrying out a large-scale 
computation is often almost entirely a matter of organization in 
which the detailed parts of the computation may be subcomputations 
which are common to many problems. Consequently, it is of great 
importance for the efficient use of a large-scale computer to have 
readily available to the programmer means for carrying out these 
subcomputations. For example, it would be very wasteful if each 
programmer had to prepare his own routine for converting binary 
numbers in the machine to decimal numbers. 

The Illiac library is a collection of routines which have 
been prepared to make programming easier o They have all been 
machine tested and copies punched on tape are kept available for 
useo All are prepared for use with the Decimal Order Input (See 
Chapters h and 5)" 

ll^„l LIBRARY ORGAITIZATION. Each Illiac routine bears 
a label consisting of a letter (or possibly two letters) followed 
by a number. The library has been divided into categories according 
to the types of operations performed, and the letter In a label 
denotes the category to which the routine belongs » Occasionally a 
routine, though belonging to one category, is useful only in 
connection with operations in another category (as, for example, a 
floating decimal auxiliary), and in such cases the letter of the 
second category has been added to the label. The number in the label 
is simply a serial number within the category. 

li^-1 



An accession serial number is also to be found on each 
routine; this is a Laboratory reference number and has no connection 
with the categories . 

^^•2 PROGRAM DESCRIPTIONS A ND SUMMARY ATTOTgT.g. Each library 
routine is written up in two ways. One is a complete description of 
the method used with the list of orders of the routine. The other is 
a svmmary sheet which gives enough information to enable the routine 
to be used correctly and to be compared for usefulness with another 
routine of the same type. Programmers can obtain copies of summary 
sheets by inquiring in the Teletype Room. Complete program descriptions 
are available for reference in the Teletype Room and a user who must 
make frequent reference to a particular program or who needs to modify 
it for his own use may obtain a copy from the Teletype Operator. 

1^.3 LIBRARY TAPES, All library programs are prepared for 
use with the Decimal Order Input (Library Routine X l) and each is 
punched on heavy gray Teletype tape and filed in a cabinet accessible 
to all programmers. To use a library program the programmer need 
only help himself to the tape, copy it himself or ask the Teletype • 
Operator to assist him, and replace the library copy. 

1^.^ ILLIAC LIBRARY CATEQOR TPlFt . The categories of the 
Illiac library are as follows; 
Programmed Arithmetic 

a: Floating Point 

B. Other Programmed Arithmetic 
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Code Checking 



C . Post mortem checks 

D. Dynamic code checks 

Integration 

E . Quadra tiire 

F. Ordinary differential equations 

G. Partial differential equations 

Operations on Functions 

H. Zeros and minima 

I . Interpolation 

J. Operations on polynomials and power series 

K. Approximations and statistics 

Linear Algebra 

L. Simultaneous linear equations 

M. Other operations on matrices and vectors 

Input and Output 

N. Number input 
0. Scope output 
p. Printing and pimching 

Mathematical Logic 

Q. lyfetthematical logic 

Particular Functions 

R. Roots and fractional powers 

S. Logarithmic, exponential and hyperbolic functions 

T. Trigonometrical functions 

V. Other special functions 

Organizational 

W. Counting, sorting and selecting 
X. Program preparation 

Miscellaneous 

Z. Miscellaneous complete programs 
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1^-5 ILLIAC LIBRARY Library routines which are of most 

interest or most heavily used are placed in the active library, while 

seldom-used routines are placed in the inactive library. On April 1 

1955^ the active Illiac library contained the following routines: 

LABEL SERIAL TITLE 

Floating Becimal Arithmetic Routine 
Floating Decimal with Standard Auxiliary 
Convert a Number from Floating Decimal 
Representation to a Normal Machine Form 
1.7 Precision Floating Binary Arithmetic 
with Floating Decimal Conversion 
Complex Number Arithmetic 
Floating Decimal Routine and Auxiliaries 

Post Mortem of the Decimal Order Input 

Post Mortem for Fractions and their Locations 

Post Mortem for Integers and their Locations 

Print Sets of Order Pairs and their Locations 

Address Search Routine 

Drum Sexadecimal Post Mortem 

Drum Fraction Post Mortem 

Print Sets of Floating Decimal Numbers and 

their Decimal Locations 

Combined Post Mortem 

Drum Order Pair Post Mortem 

Check Point Routine 
Control Transfer Check 

Integration by Simpson's Rule 

Solution of a System of Ordinary Differential 
Equations 

Solution of a System of Differential Equations 
by Milne's Method 

Integration of n Simultaneous Second Order 
Differential Equations with Initial Conditions 
Specified 
^Al 122 Second Order Linear Differential Equation with 

Two Point Boundary Condition 

^1 78 LaPlace's Equation - Liebmann Method 

G2 98 Poisson's Equation - Liebmann-Frankel Method 

HI 71 Inverse Interpolation 

H2 72 A Search for the Real Roots of f (x) = 



l)f-it 



Al 


63 


A2 


128 


A3 


125 


Ak 


87 


A5 


138 


A6 


15^ 


CI 


kS 


C3 


9h 


Ch 


97 


C5 


103 


C6 


lh6 


C7 


163 


C8 


l6h 


C9 


174 


CIO 


175 


Cll 


176 


Dl 


95 


d4 


70 


E2 


158 


Fl 


llif 


F2 


115 


F3 


129 



H5 
H4 


80 
81+ 


H5 
h6 


85 
86 


HFl 


111 


11 


6? 


JJ. 


159 


K2 


135 


K5 
K5 


137 

182 



L2 


113 


L3 


100 


L^ 


166 



LABEL SERIAL TITLE 

Minimization of a Function of Two Variables 
Minimization of a Function of Four Variables 
Minimization of a Function of n Variables 
Minimization of a Function of n Variables 
Treating One Variable at a Time 
Zero of a Solution of a Differential Equation 

Interpolation 

Roots of a Polynomial 

Product Moment Correlations, Means, Standard 
Deviations, Variances and Covariances 
Least Squares 
Autocorrelations 

LI 112 Solution of a Set of Simultaneous Linear Algebraic 

Equations 

Automatic Linear Equation Solver 

Complete Linear Equation Solver 

Half -precision Solution of a Set of n Simultaneous 

Linear Algebraic Eqimtions 
L5 165 Ha If -precis ion Automatic Linear Equation Solver 

Eigenvalues and Eigenvectors of a Symmetric Matrix 

Automatic Inversion of a Symmetric Wbtrix 

Automatic Eigenvalue-Eigenvector Program 

Closed Eigenvalues and Eigenvectors 

Solution of the Determinantal Equation j A - A3 | = C 

New Matrix Multiplication 

Estimation of Coramunalities by the Method of 

Maximum Likelihood 

Matrix Multiplication (Closed) 

Triangularization of a Matrix and its Determinant 

Linear Programming 

Matrix Multiplication with Floating Decimal Auxiliary 

Input One Number from Tape, Integer or Fraction 
Input a Sequence of Decimal Fractions 
Drum Fraction Input 
Drum Integer Input 

Transfer Blocks of Words from the Memory to the 
Drum or from the Drum to the Memory 
N8 181 Read One N'-umber from Tape as Integer or Fraction 

02 Ih-^ Routine to Plot Points and Axes 

03 l^^-T Display Numbers and Letters on the Cathode Ray Tube 
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MO 


li+l 


M2 


118 


M3 


117 


nh 


136 


M5 


139 


m6 


1^+5 


MIO 


153 


Mil 


160 


Ml? 


173 


MI5 


183 


MAl 


IOI+ 


Nl 


61 


N2 


88 


N5 


167 


n6 


168 


NT 


172 



LABEL SERIAL 



Ok 


161 


05 


162 


06 


177 


PI 


75 


P2 


52 


P3 


53 


Fk 


55 


P5 


37 


p6 


25 


PIO 


133 


Pll 


148 


Rl 


116 


R2 


105 


H3 


106 


RAl 


92 


RA2 


99 


SI 


^3 


S2 


64 


S3 


130 


SA2 


127 


SA3 


155 


T^ 


i4o 


T5 


157 


TAl 


126 


TA2 


156 


VI 


82 


V2 


120 


Yk 


178 


V5 


184 


XI 


18 


X3 


90 


X7 


l42 


X8 


144 


X9 


169 


XIO 


170 


XAl 


121 



TITLE 

Fast Fraction Display 4x6 Raster 

Fast Character Display 4x6 Raster 

Punch to Cathode Ray Tube Conversion Program 

Print One Number Fractional or Integer in a 
Manner Determined by a Program Parameter 
Print (A) with or without Sign to n Places as 
Determined by a Program Parameter 
Print n Digits of an Integer with or without 
a Sign 

Zero Suppression Integer Print 

Print One Number in a Parameter Set Layout 

Single Column Print 

Print Headings Interlude 

Rounded or Unrounded Print 

Square Root Routine 

Integral Root Al/t 

Fractional Power Routine 

Floating Decimal Square Root Auxiliary 

Floating Decimal Cube Root Auxiliary 

Natural Logarithm 

Exponential 

Logar-lthm 

Exponential Auxiliary for Floating Decimal 
Natural Logarithm 

New Arctan X Routine 

Sine Cosine 

Sine Auxiliary for Floating Decimal 

Arc Tangent Auxiliary to the Floating Point Routine 

Legendre Polynomials 
Tchebyscheff Polynomials 
Fourier Analysis 
Spherical Bessel Functions 

Decimal Order Input 

Constant Listing Auxiliary 

Sum Check 

Assemble Blocks 

Program Interruption Routine (Without Sum Check) 

Program Interruption Routine with Sum Check 

Floating Point Constant Listing Auxiliary 
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CHAPTER 15 
DEFINITION OF TERMS 



^ A 2-5 



The accumulator register or the contents of the accumulator 
register. 



Access Time 



The time to get a number from the memory to a register or the 
time to get a number from a register to the memory. 



Accumulator 2-5 



The register into which numbers can be put by addition; same 
as A register. 

Action Cycle 

The type of memory cycle in which the memory is connected to 
the arithemtic unit for the transfer of a number. 

Addend 2-5 

The number in the memory to be added to the accumulator. 

Address 3-1 

The right-hand ten binary digits of a twenty-digit order; 
usually the location of a number in the memory. 

Analog 

Indicates a continuous operation (as opposed to digital). 

Arithmetic Unit 2-5 

The part of the computer which is required to carry out 
arithmetic operations. 

Ag 3-9 

The double-length register including A and Q but excluding 
the first digit q^ of the Q register. 

* Numbers refer to pages . 
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Augend 2-5 



The number in the accumulator to which another number is 
to be added. 



-B- 



Binary Switch 7-4 

A portion of a program which is used to cause two separate 
sequences of orders to be followed alternately. 



Bit 2-1 



One binary unit of information or one binary "digit." 
Blocking; Order 8-6 

A transfer order inserted in a program to stop the "normal" 
sequence of orders so that some check can be made on the program. 

Blunder 8-1 

A mistake in programming, or more usually a mistake in pro- 
gramming which is most obvious even without running the program. 

Boolean 

Pertaining to the symbolic logic of George Boole, the 
mathematician. 

Boots traTJ Start 5-I3 

QolTlZ Imm' "''I'' ^"^l^ '^' ^'^ °^ '^^ original order pair 
y0026 and 40000, makes it possible for the machine to continue 
Inputting and storing under the control of orders brought into 
the machine from the tape. ^ 
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Carriage Return 9-3 

The paper tape code which will cause the typing carriage of a 
printer to go to the left-hand margin of the paper. 
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Cathode Ray Tube Output 12-1 

The cathode ray tube unit which can be used to display results 
from the computer directly by output orders. 

Check Point 8-6 

A point in a routine where some check is to be made. 

Circular Right Shift I3-7 

A right shift of AQ in which use is made of the digits shifted 
into Q. 

Clear 2-k 

The operation of changing all the digits of a register to a 
common state, usually 0. 

Code 

A group of more or less arbitrary symbols used to represent 
some other group of things. 

Coding Error 8-1 

An error in the detailed preparation of a program. 
Collate 3-27 

Digitwise logical product, same as extract. 
Comparer 9-l6 

A device for checking the identity of two tapes. 
Complement 2-2 

See one's complement and two's complement. 

Conditional Transfer 3-5 

An instruction which will cause a transfer from the pattern 
of taking orders sequentially if the sign digit of the 
accumulator is 0. 

Control 1-4 

That part of the machine which serves to "control" the arith- 
metic unit, the memory and the input-output. 
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Control Counter 8-^ 

Usually the order counter. 

Control Transfer 3-5 

An instruction which may cause a transfer from the sequential 
pattern for handling orders; see unconditional control transfer 
and conditional control transfer. 



Counter 



Usually a short routine which counts the number of iterations 
carried out in some program. 



-D- 

Diagnostic Routine 

A routine for finding and diagnosing a machine fault. 
Digital 1-2 

Pertaining to digits; discrete units. 

Digitwise-complement 2-2 

Same as one's complement; a binary number is the digitwise 
complement of another binary number if and only if each digit 
of the number disagrees with the corresponding digit of the 
other number. 

Directive 5*-^ 

An order, usually on a tape, which specifies the location at 
which a subroutine is to be stored. 

D.O.I. 5-1 

The decimal order input routine. Library Routine XI; used as 
a general purpose input routine for almost all programs. 

Double Precision 

Pertaining to numbers of approximately two kO binary digit 
number lengths; hence requiring two registers or locations 
for a single number. 
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Drum 

A rotating cyclinder with a ferromagnetic coating used as a 



memory . 



Error 



Mathematically the difference "between the correct result and 
the computed result; often used in place of blunder. 

Error Frequency 13-2 

The average number of machine errors per unit time. 

Even Order 3-1 

The left-hand order of an order pair, thus using digits 2 
through 2 "19. 

Extract Order 3-2? 

An operation which puts a 1 in the Q register wherever the 
number in the designated memory location and the number in Q 
are both 1, leaving all other digits zero; also called a 
logical multiply order because the resulting number in the 
Q register is the digitwise binary product of the Q register 
and the number in the memory; same as collate. 



-F- 



Fixed Address 5-2 

The numerical address remains unaltered when put into the 
memory; usually in connection with D.O.I, where a fixed 
address is followed by an F. 

Fixed Point 2-1 

The binary point is always in the same position of a register; 
Illiac is a fixed point machine. 
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Floating Point 



Pertaining to operations in which numbers are represented 
by a number multiplied by a power of a base; thus numbers 
may have different multiplying factors; floating point 
operations in the Illiac must be obtained by programming. 



Fraction 



A number with an absolute value less than 1. 
Function Digits 3-2 

The first eight binary digits of an order; the T and V digits 



-H- 



Hang up 



An unplanned stop of the machine due either to a machine 
fault or a coding error. 



-I- 
Input Routine 



A routine for Inputting other routines, usually aiding in 
conversion from binary- decimal to binary. 



Instruction 3-2 



The operation in the machine designated by the first eight 
binary digits of an order. 



Integer 



Pertaining to integers; although Illiac is a fixed point machine 
with the binary point between the first two digits it is possible 
to carry out operations using integer numbers less than 2+^0 ^ 



Interlude 5-11 



A routine which carries out operations and is then destroyed 
by overwriting as the problem continues, usually carried out 
during input of information. 
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Instruction Code 3-1 
The order code . 
Interpretive Routine ^-9 

A routine in which a sequence of operations (instead of a 
single operation) may he carried out on a number, the sequence 
designated by parameters carried with the number. 



Iteration 



A sequence of orders, usually to be executed more than once, 
and arranged to converge to some analytic result. 



-L- 



Leapfrog 2-17 

An engineering routine which may be used to test the machine. 

Left-hand Order 3-1 

The left-hand or even order of a pair of orders, thus using 
digits 2^ through 2-19. 

Line Feed 9-3 

The paper tape code which causes a Teletype to advance the 
paper one line; always in conjunction with carriage return. 

Link ^-^ 

The part of the routine used to bring a subroutine into 
operation and designate the point of return. 

Location 1-3 

The designation of a number location, address or position in 
the memory. 

Logical Product 

When applied to numbers of more than one binary digit, the 
digitwise logical product; the same as collation of two 
numbers and the same as extract. 
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Loop 8-4 



A sequence of orders which may be carried out more than once 
automatically. 



-M- 



Machine Error 

An error caused by a fault in Illiac. A malfunction. 
Malfunction 13-2 

A mistake in a computer function. A machine error. 
Memory 1-2 

A device which stores numbers; usually the cathode ray tube 
or Williams memory in Illiac. 



-N- 



Number Register 2-5 

The temporary location in which the addend, subtrahend, 
multiplicand and divisor are automatically placed by the 
control while the corresponding order is being executed; R^ 

-0- 

Odd Order 3-1 

The rightg^and order of an order pair consisting of the 
digits 2 through 2~39. 

One's Complement 2-2 

Same as digitwise complement; each digit of a binary number 
is changed to get its binary complement. 
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Order 



A set of 20 binary or 5 sexadecimal digits which is used to 
define a machine operation. 



Order Code 3-1 



A set of sexadecimal characters used to describe the operations 

in the Illiac. 



Order Counter 8-^ 



The counter which keeps a record of the location of the next 
order pair if no transfer of control is required. 



Order Pair 3-2 



Two orders which are stored together in one location of the 
memory . 



Order Ref^lster 3-^ 



Outi)ut 



The register R^ into which an order pair is transferred from 
the memory just before either of the orders can be obeyed; R^. 



A display by oscilloscope, paper tape or Teletype printer of 
the contents of some part of Illiac. 



Overflow Digits 2-3 

The digits which, as a result of a computation, would require 
digits to the left of 2^ in a register. 



-P- 

Page Printer 9-9 

An automatic typewriter. 
Paper Tape 

The common II/16 inch width paper tape used for handling 
numbers in and out of Illiac. 
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Parameter 4-7 

A number which may be altered from one computation to another 
but is held fixed during a single run. 

Partial Substitution Order 3-I7 

An order which allows the address digits only to be stored 
in the memory; a ^2 or k6 order. 

Plant k-h 

The act of putting some number into an order or routine; 
usually with reference to an address. 

Playback 

The number coming from the magnetic drum memory. 



Position 



Referring to one of the 2^ to 2"^^ places of the memory or 
of a machine register. 

Post Mortem 8-8 

The checking of the routine in the memory after it has been 
used, usually by an automatic routine which compares the 
contents of the memory with the original input tape. 

Program 1-4 

The plan of a calculation. 

Punch 9-1 

A unit to punch holes in II/16 inch paper tape, usually the 
punch providing output from Illiac. 



-Q. 



^ or g Register 2-5 



The register which must be used for the multiplier during a 
multiplication; also the register into which the quotient is 
placed during division; the contents of the Q register. 
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Quotient Register 2-5 

Same as Q register. 

Quadrature 

The process of evaluating a definite Integral by numerical means 



-R- 



The order register; bottom row of kO neon lights on Illiac. 

Range 

The set of numbers which may be handled In an Illiac register 
directly, from and Including -1, up to and Including 1 - 2"^^. 

Raster 

The array on a memory tube. 

Read- around 

An Index of the Interference of one spot with another In 
the Williams memory. 

Reader 9-1 

The device to take data from a punched paper tape for Insertion 
Into Illiac. 

Record 

The process of storing a number on the magnetic drum. 

Regeneration 

The process of refreshing the stored information on a 
Williams tube. 

Register 2-5 

A row of ^0 fllpflops which can hold a binary number. 
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Relative Address 5-2 

An address which is relative to a directive and hence must 
be added to the directive address to get the true memory 
location. 

Remainder 2-12 

The quantity left from the dividend after a division process 
(without regard to Illiac). 

Reperforator 9-11 

A standard Teletype machine for making paper tapes from 
electrical signals. 

Residue 2-13 

The quantity in A after a division order has been obeyed 
(in the Illiac) . 

Right-hand Order 3-1 

The order using positions 2"^^ through 2'^^; same as odd order. 
Round-off 2-12, 13 

The process of adding 2*" in multiplication and making q = l 
in division when the term is used with Illiac. 39 

Rout ine k-1 

A completed sequence of orders in coded form. 
Run 13-2 

A period of time during which a program is bein^ executed 
by the computer. 

-S- 
Scale 6-1 



The adjustment of a number to come within range of 111 



lac 
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Sexadecimal 3-2 

A number system with "base l6. 

Shift 

The process of moving a number to the right or left in a 
register. 

Sign Digit 2-2 

The first position in a register; 2 . 

Single Address Code 

A machine order code which has one address only with each 
order; Illiac has a single address code. 

Slave Tube 8-^ 

A device for displaying the entire contents of one memory 
cathode ray tube. 

Sjbore 

The memory, usually Williams memory. 

Subroutine ^-1 

A routine arranged according to a standard pattern so that 
it can be easily used as part of other routines-. 

Sum Check 13"^ 

A check on computer reliability made by adding a set of 
quantities twice and comparing the two results. 

-T- 



Tape Code 9-1 

The hole patterns in a paper tape that are used to represent 
numbers and instructions. 
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Tape Comparer 9-16 

A device to compare two paper tapes. 
T Digit 3-16 

The first sexadecimal digit of an instruction. 
Temporary Storage 

Locations in the memory used briefly during a calculation 
and not assigned to any number to be retained for results. 

Terminating Symbol 

A symbol on the tape indicating the end of a section of 
code or tape . 

Two's Complement 2-2 

The difference between 2 and the number whose two's complement 
is to be found; in the Illiac identical, because it is modulo 
2, to reversing each digit and adding 2-39. 



-U- 

Unconditional Transfer 3-5 

A transfer out of the ordinary sequential pattern of handling 
orders, regardless of the sign of A. 

Useful Time 13-2 

The average useful time per run is the probability that no 
error will occur multiplied by the length of the run. 



1 Dip^lt 3-6 

The second sexadecimal digit of an instruction, 
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Waste Order ^-3 

An order that serves no computational P^^P°^^ .^^^^!;%. 
inserted because of ease in coding or limitation of the 
machine because of odd and even pairing of orders. 

Williams Tube 

A cathode ray tube used for storing binary digits . 

Word 

Forty binary digits In a register or a single location^ 
in the memory. 

Working Space 

See temporary storage. 
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CHAPTKR 16 
THE DRUM STORAGE UNIT 

16.1 GENERAL DESCRIPTION . Auxiliary storage of 12,800 
words is provided for the Illiac by the magnetic drum storage unit. 
Some knowledge of the physical characteristics and mode of operation 
of the drum storage unit is necessary for its most efficient use. 

The magnetic drum is, physically, a rotating cylinder whose 
surface is coated with magnetic material. Lengthwise, the cylinder 
is subdivided into 200 tracks; associated with each track is a re- 
cording and reading head capable of recording digital inforrriation 
on the magnetic surface and capable of reading data previously 
recorded (Page 16~2). The periphery of the cylinder is subdivided 
into 64 sectors; as a sector passes under one head, the digits of 
one word are recorded or read. Thus, the selection of any one of the 
12,800 words requires a spatial selection of a track and a temporal 
selection of a sector; that is, access to one word requires selection 
of one of 200 heads and the selection of one of 64 time intervals 
during which the digits of the desired word pass under the selected 
head. 

16.2 COMMON DRUM INSTRUCTIONS . The record instruction 
86 11 OOp places the contents of A in drum location p and performs 
eleven left shifts of AQ. The playback (or read) instruction 

85 11 OOp performs eleven left shifts of Q and places the contents of 
drum location p in A. Either type of drum instruction must occupy a 
complete 40 digit word. 
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FIG. 16.1 ILLIAC MAGNETIC DRUM 



The address p of either drum instruction is specified by 
the rightmost 14 binary digits (2~ to 2~^^) of the 40 digit 
instruction and hence is obeyed modulo 2 = 16,384. It is restricted 

to the range < p < 12,800 for playback and 2560 < p < 12,800 for 

ol4 
record, with the interpretation in both cases being made modulo 2 . 

Of the 14 binary digits specifying p, the least significant six 

(2"^^ to 2 '^) determine in which of sixty- four sectors the desired 

word lies; the remaining eight (2" to 2"^^) specify the track. 

Except in special circumstances, the programmer need not be concerned 

with the breakdown of the address p into track and sector addresses, 

since the track switching is automatic and imposes no special timing 

restrictions on drum use. 

The section of the drum < p < 2560 is used for storage of 
the drum bootstrap and commonly used routines; no recording by drum 
instructions is possible in this range. If either the record or 
playback instruction has an address in the range 12,800 < p < 16,384 
(mod. 16,384), the instruction will cause the computer to stop. 

For most efficient use of the drum storage unit, the desired 
sector of the drum should be approaching the reading and recording head 
when a drum instruction is to be executed. Subroutines for transfer of 
blocks of words between the drum and electrostatic memories require time 
for execution of other Illiac instructions between successive drum 
accesses; to facilitate such use, sectors whose Illiac addresses differ 
by one are physically five sectors apart on the drum surface (Fig. l6.l). 
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The circuits for switcfdng from one track to another are 
so designed that all „ords on the sa.e sector are equally accessible. 
For this to be possible, a delay of one „ord ti.e precedes any drum 

consultation. 

As an example, suppose that the dir;its of word have been 
read from irur. track 0, sector 0. As the drum rotates (Fig. 16.1), 
sectors 13, 26, and 39 pass under the readin, heads. This interval 
Of three word times is available for calculation and for address 
modification. The dn.m instruction requesting; consultation of 
sector 1 Should be in the lUiac order register before sector 52 begins 
to pass under the reading heads. The enforced delay of one word time 
occurs as sector 52 passes under the head, whereupon the digits of any 
word in sector 1 can be read into the Illiac accumulator as sector 1 
passes under the heads. 

16.3 CALCUL.^lu, OF ACCKSS TTKK. The timing data necessary 
for efficient use of the drum storage unit can be deduced from the 
following facts: 

1. The time required for one revolution of the drum is 16.9 milliseconds, 

2. Since the digits of 64 words are stored serially around the periphery 
of each track, the time required for the digits of one word (one 
sector interval) to pass a reading or recording head is 0.264 
milliseconds. 

3. The d^m instruction must be in the Illiac order register at least 
one sector interval (0.264 milliseconds) before the digits of the 
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word can be transferred between the drum and the Illiac accumulator. 

Thus, the minimum time required for execution of a drum instruction 

is two sector intervals or 0.53 milliseconds. 
4. words corresponding to successive addresses in drum instructions 

are spaced five sectors apart on the drum surface. 

An isolated drum instruction with no definite timing relation- 
ship to other drum instructions may be regarded as occurring when the 
drum is in a random position. In this case, the time taken by the drum 
instruction will be equally likely to lie anywhere in the range between 
the minimum of 0.53 milliseconds and the maximum of one revolution time 
plus one sector interval or 17.43 milliseconds. The average time would 
thus be 9.0 milliseconds. 

A program does not usually have isolated drum instructions, 
but instead it will contain a sequence of drum instructions separated 
by other Illiac instructions. Under these circumstances, three timing 
calculations are of importance. 

1. The time t between the completion of one drum instruction and the 
completion of the next drum instruction. 

2. The time T required by all the Illiac instructions which come 
between the two drum instructions. 

3. The maximum time T available for execution of Illiac instructions 

max 

if the second drum instruction is to be obeyed as quickly as 
possible after the first. 

The commonest case is that in which successive drum instructions 

have successive addresses. The time T^^^ is then equal to three sector 

max 
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intervals or 0.79 milliseconds. If T < T then the time t is 

m£LX 

equal to five sector intervals, or 1.32 milliseconds. If the time T 

required for execution of Illiac instructions exceeds T , at least 

max 

one additional drum revolution will be required for each access, 
giving t = 18.2 milliseconds as the time between completion of drum 
instructions. 

A general calculation of t and T for arbitrary addresses 

max 

of two successively obeyed drum instructions can be made in the following 

way. If the address of one drum instruction is p, , and that of the next 

is p^, one first calculates ^(Pg - P-j_) - 2 - 3.8T and forms its residue 

modulo 64. (All times are in milliseconds.) Call this quantity R. In 

other words, R is formed by successively adding or subtracting 64 from 

the above expression so as to make R lie in the range < R < 64. The 

time t between completion of the two drum instructions is then given by 

the formula t = 0.53 •*• T f .264R. The time T is chosen as the 

max 

shortest time T to make R = 0. Thus one obtains 

'^max " -264 {Residue of [^{p^ « p^) - 2 ] (modulo 64 J} 

16.4 CALCULATION OF TIME OF EXECUTION OF ILLIAC INSTRUCTIONS . 
If access time to the drum storage is to be minimized, it is important 
that accurate calculations be made of the time of execution of Illiac 
instructions obeyed between successive drum instructions. The most 
accurate way of performing this calculation is to determine the number 
of clock periods of the Illiac electrostatic storage unit. The clock 
periods required for instructions commonly used in drum transfer routines 
are given in Table 16. 4. 
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Readout of each order pair 

Order types K, S, F, and L 

Order types 2, 4, 5, J, and 3 v^en 

control is transferred. 
Order types 85, 86, and 3 when control 

is not transferred. 







Table 16.4 



Two clock periods must be allotted for each order pair, 
including the drum instruction itself and regardless of whether 
one or two orders are obeyed. The remaining clock periods listed 
in Table 16.4 are in addition to those allotted for order pair 
readout. The clock period of the Illiac storage unit is 
0.0185 ± 0.0002 milliseconds; the maximum time T required for 
executing instructions between successive drum orders is found by 
multiplying the total number of clock periods required for the 
instructions by 0.0187. As an example, the time required for the 
record loop in library routine Yl is calculated below; 
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Location 
10 



14 



Order 

40 11 L 
L5 (n) F 



11 


86 11 


F 




00 


F 


12 








F5 11 


L 




40 11 


L 


13 








F5 10 


L 




42 10 


L 



FO 35 L 
32 10 L 



2 readout order pair 
4 order type L 



2 readout order pair 
order type 86 



2 readout order pair 
4 order type F 
2 order type 4 



2 readout order pair 
4 order type F 
2 order type 4 



2 readout order pair 

4 order type F 

2 control transfer 

32 Total 



The time T required for the record loop is therefore 32 x 0.018? =^ 0.598 
milliseconds, which is less than the time of three sector intervals 
(0,792 milliseconds). Successive words will therefore be recorded at 
the rate of 1.32 milliseconds per word, after the first word has 
been recox'ded. 

16.5 arHER TI?i!;S of drum INSI'RUCTIONS . a more general 
type of drum instruction may be written in the form 8 V, 11 T V^ p. 
If T is 0, 1, 8, or 9 the drmn instiniction will be obeyed as described 
in the first section. If T is any other digit type the right hand 
half of the instruction will be obeyed as a separate Illiac instruction 
after the drum instruction has been obeyed. In such a case the drum 
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, , A,.^r^ A 1Q2 = 2^-^ while the drum instruction 

address is interpreted modulo 8,1V^ '^ wnxx« 

is being obeyed and far this reason no hangup is possible. The 

rightM^st binary digit of V, is uniquely determned by the drum address 

p but Vt is othervd.se arbitrary. 

Possible choices for V^ include 5, 6, 7, J and L. The 

playback instructions have V^ equal to 5 or J and both produce the 
same result. It is customary to take V^ = 5 for pUyback. The 
recoil instructions have V^ equal to 6, 7, and L, and they record N(A). 
.ero, and l/2 respectively as my be deduced from the characteristics 

of the variant V, . 

If V is takai equal to 4 a playback is executed but since A 
is not cleared before carrying out the playback the result will also 
depend on the initial contents of A. The result of this instruction is 
to shift A and Q left eleven places but during the shift to insert 

^eros into right hand eleven digits of A rather than digits from Q. 

A digitv^se inclusive or (or logical sum) of this shifted quantity in 

A is formed with the word taken from the drum. 

A number of other instructions exist which affect the drum 

but which also have little utility. They include the cases: 

a. Drum instructions with the left hand address different from 11. 

b. Drum instiMctions with the left hand type digit equal to 9 

rather than 8. 

c. Drum instructions lying entirely within the right half of a 
word. (These instructions will pix^duce the conventional effect 
only if the address equals eleven mod 64.) 
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Since applications of these types of instructions to 
actual programs are so limited we give no description of their 
properties here. 

1^-6 MODIFICATION OF ILLIAC . The initial instruction 
pair 8002840000 which has been used to start the tape bootstrap was 
rewired so that the instruction pair 8500S40000 is used instead. 
This change will have no effect in the operation of the computer 
except when a bootstrap start is used, and even then the computer will 
usually respond in the same way now as it did with the original 
instruction pair. The differences may be noted as follows: 

1. No D.O.I, is now required on the front of your tape since if the 
D.O.I, is omitted it is taken from the drum automatically and 
placed in locations 3F7 to 3LL. The additional locations 3F5, 
3F6, 000, and 001 are used in this process. 

2. If a tape bootstrap such as that appearing on the D.O.I, tape is 
at the beginning of your tape, the computer will handle it in 
the same way as it has in the past except that locations 3L6 
through 3LL will be used before the bootstrap is read. 

3. If your program requires the reading of a jump instruction by means 
of a bootstrap start, it must be altered. A jump instruction 
written as 24O64OOOOO must be replaced by 00K24100N, while if the 
instruction 263F700000 is to be read on a bootstrap start this 
instruction may be omitted entirely. In any such case, the 
locations 3F5, 3F6, 000, 001 will be used and the D.O.I, will be 
replaced in the Williams memory. 
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The program which provides the facilities described above 
consists of a set of routines which are stored on the drum. The 
initial instruction pair 8500S40000 plays back a drum bootstrap which 
occupies locations 000, 001, 400, 401 on the drum and 000, 001 in the 
high speed memory. It, in turn, plays back routine 1 from 3L6-3L9 
on the drum to 3L6-3L9 in the high speed memory. Routine 1 plays 
back routixie 2 from 3LK-3LL on the drum to 3U-3LL in the high speed 
memory. Routine 2 records the contents of the high speed memory 
locations 002-OlL at locations 31F2-31LL on the drum to preserve them 
and returns control to routine 1. Routine 1 plays back routines 3 
and 4 from locations 002-013 on the drum to 002-013 in the high speed 
memoiy anql jumps to routine 3. Routine 3 reads one se^decimal 
character from the tape to deteraine which routine should be played 
back next and prepares routine 4 accordingly. If the character is 
read from the tape, routine 4 will play back part of the D.O.I, and 
routine 5 and will jump to routine 5- Routine 5 replaces the words 
in 31F2-31LL in locations 002-01L in the high speed memory and then 
plays back the rest of the D.O.I. Entry is made to the D.O.I, in 
such a way that the reading of the initial function digit of the 
directive is suppressed the first time. This is necessary since 
routine 3 has already read the initial 0. Routine 5 occupies locations 
3F5, 3F6, 3L9-001. All except locations 3F5, 3F6, 000, and 001 are 
covered by the D.O.I. 
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If routine 3 reads the character 8 from the tape a tape 
bootstrap is indicated and routine 4 is prepared so that it plays 
back routine 6 into locations 3L8-3LL and jumps to routine 6. All 
words of routine 6 are contained in routine 5, but entry is made at 
a different point. Routine 6 replaces the words in 31F2-31LL in 
locations 002-OlL and causes the next nine characters to be read 
from the tape and stored in 000 after faking the initial 8 which was 
already read. Routine 6 then .jumps to 000 and the usual tape boot- 
strap sequence has been started. 

Other initial characters have been assigned as follows: 

(l) Leapfrog 

(3) Flying Leapfrog 

(F) Drum Post Mortem Routines 

(L) Store Post Mortem Routines 

16.7 USE OF THE DRUM. The drum will most commonly be 
used by means of library subroutines, men this is done no special 
knowledge of drum instructions is required, and all problems of 
reducing access time to the "minimum" of 1.32 milliseconds are handled 
automatically by the subroutine. 

Y-1 is a typical subroutine for transferring information 
back and forth from the drum. It may be used to record a block of r 
words on successive locations of the drum starting at location q 
provided these words are present in successive locations in the 
Williams memory starting at location n. One must enter the subroutine 
by means of the instructions 
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JO n 
P 50 p 

p-fl 26 — 
00 q 

p+2 00 r 
in order to cause the recorxiing to occur. Control will be returned 
to the right hand side of p+2 after the recording. 

The sajne subroutine may be used to play back r words from 
successive locations of the drum starting at location q provided these 
words were previously recorded by the subroutine. They will be 
stored in successive locations of the Williams memory starting at any 
arbitrary location n. In order to perform this playback and store 
operation one enters the subroutine by means of the instn^ctions 

50 n 

P 50 p 

p+1 26 ■— 

00 q 

p+2 00 r. 
In addition to recording and playing back a block of words 
in minimiun time, this subroutine comp^xtes a sum check for the block 
of words, which is stored on location q + r of the drum. Upon play- 
back this sum check is recomputed and compared with the recorded sum 
check. If any word has been recorded or played back incorrectly, the 
two sum checks will not agree and the computer will stop on the 

instruction FFOIO. 

In some complex problems one may wish to play back words 
,^ich do not occupy successive locations on the drum, and if this is 
the case library subroutines are not practicable. For e^mple, it 
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my be desirable to record a „at^ix in such a way that either a 
row or a column can be played back with minimum access time. If 
the matrix has no more than 65 columns, one may achieve minimum 
access for either rows or columns by recording the elements of each 
row in successive locations but recording successive rows 65 locations 
apart. Thus, we see that T^^ (as described in 16.3) is 1.32 milliseconds 
for both rows and colu^^s and corresponds to the usual "minimum" access 
time. If the number of columns N Ues in the range 65 < N < 129 one 
must record the rows 129 locations apart, etc. 

If a technique of this sort is to be used with the drum, it is 
up to the programmer to write his own record and playback loops with 
due attention to access time and to sum checking. 
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-A- 



A register, 2-5 
Acciamulator, 2-5 

shifting without shifting Q, 7-lo 

Accxaracy, 6-T 

Add from Q, 2-8, 5-26 

Addend, 2-5 

Addition, 2-5^ 2-7* 3-29 

Address, 1-5 > 3-1 

digits, 3-8 

fixed, 5-1 

relative, 5-1 

AQ, 5-9 

Arithmetic checking, 15-15 
Arithmetic unit, 1-2, 2-5 
Augend, 2-5 
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Binary 

chopping, 7-15 

-coded decimal, 9-5 

digit, 2-1 

point, 2-1 

point shifted, 6-1 

switch, 7-^ 
Blocking order, 8-6 

example of use, 11-57 
Blunders, 8-1 

list of typical, 8-5 

discovery of, 8-2 
Bootstrap, 5-12 

start, 5-15 
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CameiTa operation, 12-^4- 
Cathode ray t\ibe display, 12-1 

orders for controlling, 12-5^ ^ 
programming for, 12-5 
example of result, 12-9 



Check, 

programmed, 15-1 

sum, 13-6 
Checking, 8-1 

arithmetic, 15-15 

drum, 15-1^ 

input, 15-6 

memory, 15-12 

methods, table of, 15-17 

of a particular program, 11 -56 
Check point, 8-6 
Coding tricks, 7-1 
Combination Teletype, 9-11 
Comparing of tapes, 9-16 
Complement, 2-2 
Conditional transfer, 5-5 j» I6 
Constants, orders and addresses as, 7-7 
Control coimter, Q-k 
Control transfer, 5-5 

reversed, 7-5 
Control -unit, l-^i- 
Convergence criteria, 7-12 
Correction tape, 8-5 
Cyclotron trajectories, 11-1 
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Data tape, I5-9 

Decimal order input, 5-1 

Definition of terms, 15-1 

Differential equations, 11-17 

auxiliary subroutine, 11-20 

Digital, 1-2 

Digitvrise complement, 2-2 

Directive, 5-^ 

Division, 2-12, 5-19 

"double precision", 2-l8 
hangup, 2 -16 
remainder, 2-15, 7-1^ 
special cases, 2-l6 

Drum, 

checking, 15-15 
playback, 5-21, 15-1^ 
record, 5-22, 15-1^ 
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Error frequency, 15-2 
Extract, 5-2? 
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FF order, 5-28, 15-8 
Fifth hole, 9-^ 
Fixed point, 2-1 
Format characters, 9-5 
Fractions, input of, 5-° 
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Glossary, 15-1 

Hjand punch, 9-lT 

-I- 



Increment add, 2-8, 5-28 

from Q, 2-8, 5-25 
Information sources, 8-i^■ 

Input, 9-1 ^ 
Input checking, 15-6 
Input-output orders, 5-'^!^ "^0 
Institute for Advanced Study, l-i 
Instruction, 5-1 
Integer operations, 2-20 
Interchange of numbers, T-H 
Interlude, 5-11 
Interpretive routine, ^"3 
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Keyboarcl perforator, 9-9 
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Leapfrog, 2-17 

Left shift, 2-i^, 2-9, 5-12 > 7-18 

Letter printing, 9-4 

Library, li<-l 

Link, h'k. 

Location, 1-5 

Logical product, 5-27 

Loop, 8-4 
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Machine error, 13-2 

Magnetic drum memory, 16-1 

Malfunction, 15-2 

Marking, 7-15 

Master tape, 15-9 

Mean time between errors, 15-2 

Memory, 1-2 

checking, 15-12 
to Q, 2-17, 5-18 

Modifications to program, 8-4 

Multiplications, 2-11, 5-20 
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Negative numbers, 2-2 
Number register, 2-5 
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One-address code, 5-1 



Order, l-k 

as a constant, 7-7 
blocking, 8-6 
code, 5-1 

execution time, 10-2 
pair, 3-2 
types, 3-6* 7 
variants, 5-8 

Ordvac, 1-1 

Output, 9-1 

Overflow, 2-5 
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Page printer, 9-9 
Parameters, k-'J 

preset, ^-11 ;, 5-7 

program, k-'J 

use of, 11-16 
Photo-electric tape reader, 9-17 
Playback from drum, 15-1^ 
Poisson distribution, 15-5 
Polynomial evaluation, J-lG 
Post mortem routine, 8-8 
Product summation, 7-1 
Program, 1-^ 

complete, 11-1 

modification, 8-^ 

tape, 15-9 
Programmed check, 15-1 
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Q register, 2-5 
Quotient register, 2-5 

for interchanges, 7-11 



-R- 



R^, 2-5 



Record on drtim^ 13-lij- 
Register, 2-5 
Remainder, 2-12 

in integer division, 7-14 
Reperf orator, 9-11 
Resetting of cycles, 7-8 
Residue, 2-13 

Right shift, 2-4, 2-9, 3-I3 
Routine, 1-4 

address search, 8-11 

check point, 8-8, 12 

closed, 4-2 

entry to, 4-3 

exit from, 4-3 

interpretive, 4-9 

post mortem, 8-8, 9 

sequence checking, 8-8, 11 
Run, 13-2 
Running time, 10-1 

simple formula, 10-4 



-S- 



Scaling, 6-1 

adjustable, 6-10 

continuous, 6-10 

of a complete problem, 11-10 
Sexadecimal, 3-2 
Sign digit, 2-2 
Starting a program, 5-6 
Stop orders, 3-6 
Store, 2-17, 3-17 
Subtraction, 2-6, 2-7, 3-29 
Sum check, I3-6 



-T- 



Tape, 



data, 13-9 
master, 13-9 
program, I3-9 



Tape code, 9-1 

caniplete, 9-5 
orders for, 9-6 
sexadecimal, 9-1 
Tape preparation, 9-1 
Teletype equipment, 9-7 
Testing for and -1, 7-6 

for 1/2, 7-11 
Transfer of control, l-5> 5-5 



-U- 



Unconditional transfer, 3-5, 5-1^ 



UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



ILLIAC ORDER CODE SUMMARY 

Input output and drum orders are not included. 
N = Storage location specified in order r 



Add, Subtract or Copy into A 

ir n III i T I I r T : i ri j ■ 



7^ microseconds 



I vUvd 






Inst Replace A by 

I I III I 1 i w i irt ii II I T I 

lA A 4- N 
LO A - N 

l£ A + |ir| 

L2 A * |n| 

A+ (N+2"?^) 

A * (N+2"^^) 

*39| 



FO 

f6 



A 4- |K+2 

A - iN+a"*^^ 



^ Ah- Q 

SO A - Q 

S6 A + |q| 

S2 A - |Q| 

Yih A + (Q+2"^^) 

KO A * (Q+2'*^^) 

kS a + |Q+2*^^| 

K2 A - |Q+2""^^' 



Inst Replace A by 



L5 
LI 



LT 



FT 
F5 

S5 

SI 

ST 
S3 







10 (Q+S"*^^) 
KT |Q-f2"^^ 



Inat Replace A by 

LJ 1/2 + N 

L9 1/2 * N 

LL 1/2 -f |n| 

LS 1/2 * |n| 

FJ 1/2 + (N+2"^^) 

F9 1/2 •* (N+2'^^) 

FL 1/2 + |N+2''^^| 

FS 1/2 - iN+2"^^| 

SJ 1/2 + Q 

S9 1/2 * Q 

SL 1/2 + |q| 

ss 1/2 w |q| 

KJ 1/2 + (Q+2"5^) 

K9 1/2 - (Q+2*"^^) 

KL 1/2 4- 1(^2*^^! 

KB 1/2 • |Q+2"'^^| I 






Oi ^ 






J I ^l :'(x 



•^/l,'c-t ■■■I--— ^.<._t. ...--' '- ■• • ^ . * -^ "^^ 






y '■.^' 



■y*^ 



Transfer from A 55 microseconds 



-I a 



Inst Effect 



Inst Effect 



i^O 



N =: A 



k6 l.h.a. JJ = A 
k2 r.h.a. N = A 



in 

h3 



N = A = 



Inst Effect 



1 


,h 


^a« 


N 


__ 


0, 






A 


= 







r 


•h 


•a* 


N 


= 


0, 






A 


= 








^9 N - A - 1/2 

kh l.h»a. N = 0, 
A = 1/2 

ij-S r*li-a- N = 0, 
A = 1/2 



Copy into Q 55 microseconds 

Instruction T^ffect 

50' (52,5i4-,56) Copy K into Q 
51 (55,55i5T) 
59 (5S,5Jr5L) 

JO (J2,jit-,J6) 
jl (J55J5,JT) 
J9 (jS,JJ,Jli) 



" " " " and clear A 

" ♦' «* '' and put 1/2 in A 

Logical product N and Q to (^ 

ft It rt »f »' " '^ and clear A 

♦f " »» " " " " and put 1/2 in A 



Multiply setting % = 

HOLD 
Replace AQ "by 



Inst 



1^ 
70 

T6 
72 



Q.N .+ 2 -^^k 

-Q.N + 2 -^^A 

Q|n1+ 2"'^^A 

-q1n1+ 2*^^A 



Inst 
T5 

7i 

77 
75 



700 microseconds 

UNROUNDED 
Replace AQ by. 



Q.N. 
-Q.N. 

q|n1 
-qInI 



ROUNDED 
Inst Replace AQ by 

-4o 



7J 
79 
7L 
7S 



Q.N. + 2 
-Q.N. + 2 
q1n| + 2 

-qInI + 2 



-1^0 
-1^0 



- 2 - 



Divide leaving residue in A 800 microseconds 

iJ^Qt Replace Q by 
66 (AQ)/Tf 
■ 67 2-^Vn if Q > 0; 2"^9 (i^Q)/jj if Q < 
^ (6S) [1/2 ^ 2-^Vn if Q > 0; [1/2 -. 2-^^(uq)] /n if Q < 

Shift A9, Address = Number of shifts (mod 6k) I 
Ion microseconds 



Ittstructions Effect 

00 

(02.04.06) Shift left 

01 Caear A 

(03.05.07) and shift left 

09 Put 1/2 in A 

(OS,OJ,OL) and shift left 



Instruction 

10 
(12,U,16) 

11 
(15,15,17) 

(1S,1J,1L) 



Effect 



Shift right 

Clear A 
«tnd shift right 

Put 1/2 in A 
and shift right 



^Transf er Control UWCONDITIOKAL TRANSFER 55 microseconds 



last Effect 
26 Jump to l.h. 

22 Jump to r.h. 

24 Stopi to l,h. 

20 Stop J to r.h- 

CONDITIONAL TRANSFER 



Ijnst 



Effect 



27 Jump to l.h., 

(37) A = 

23 Jump to r.h., 

(33) A = 

25 Stopi to l.h., 

(35) A = 

21 Stopf to r-h., 

(31) A = 



Inst Effect 

2L J\j[Brp to l.h.. 
(3L) A = 1/2 

2S Junrp to r.h., 
(3d,) A = 1/2 

2 J Stop; to l^h., 
(3J) A = 1/2 

29 Stop J to r.h., 
(39) A = 1/2 



55 microseconds (15 if no transfer) 



last 



Effect 



% 
32 

34 
30 



Junrp to l.h. if A > O 
Jump to. r ,h * if A > 
Stopj to l.h. if A > 
Stop; to r.h. if A > 

- 3 - 



Final Stop 
OF 

S. bill, 2/26/54' 
rt/ 9/26/56 



<^C>^ \J[h)' >VCh) 



^ •^jJ[f^J^^-^-~'-jgf 









** " \ ^.^ ^.a.-A<^-V *^^-'^ 



9 Z -^ 
^9 



/- 






^^ ^ ^yi4^.-^/» -/-•'^^A-- 



<> -< tt 



i2c'0 






